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ABSTRACT 


Techniques  for  sensing  the  relative  threat  of  flow  instabilities 
in  jet  engines  have  been  a  long  sought-after  goal.  This  is  the  final 
report  of  a  study  aimed  at  determining  the  feasibility  of  certain  types 
of  fluid  jet  amplifiers  used  as  fast  pressure  and  temperature  sensors. 
The  background  of  the  project,  the  environment  in  which  the  fluid 
sensors  must  operate,  and  initial  results  have  been  presented  in  the 
first  annual  report. 

The  more  theoretical  aspects  of  the  present  report  include  a 
study  of  the  statics  and  dynamics  of  the  interaction  of  a  fluid  jet  and 
a  receiver  port,  the  dynamics  of  the  controlled  deflections  of  jets  in 
fluid  amplifiers,  and  the  dynamic  analysis  of  fluid  systems.  The  more, 
applied  aspects  include  the  design  of  a  fluid-relay  stall-cell  sensor,  a 
fluid  diode  possibly  useful  for  sensing  blade  wakes  or  stall  cells,  and  a 
unique  reference -frequency  pneumatic  oscillator. 
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1.  INTRODUCTION 


Techniques  for  sensing  the  relative  threat  of  flow  instabilities 
in  jet  engines  have  been  a  long  sought-after  goal.  The  present  study 
has  been  directed  toward  the  possible  use  of  fluid  jet  amplifiers  with 
no  moving  parts  to  sense  stall  cells  and  blade -wake  widths,  and  to  ex¬ 
tend  the  basic  knowledge  of  these  complex  6 -port  amplifiers  which  also 
have  many  other  military  and  non-military  uses.  In  addition,  a  pneumatic 
oscillator  was  developed,  using  some  moving  but  non-sliding  parts,  which 
has  a  frequency  largely  insensitive  to  pressure  and  temperature  changes. 

The  background  of  the  present  project,  the  environment  in  which 
fluid  sensors  must  operate,  and  initial  results  under  the  contract  are 
presented  in  the  first  annual  report. 

The  Engineering  Projects  Laboratory  depends  on  Master's  and 
Doctor's  degree  candidates  for  the  bu.lk  of  its  effort.  It  is  impossible  to 
accurately  predict,  let  alone  control,  the  durations  of  several  adequate 
thesis  investigations,  particularly  at  the  doctorate  level.  Consequently, 
only  two  Master's  theses  have  been  completed  under  the  contract  (a  third 
was  abortive),  while  two  most  pertinent  doctoral  theses  and  one  Master's 
thesis  are  unfinished,  a  fact  reflected  in  this  report.  Estimated  completion 
dates  for  all  three  are  in  August  of  this  year.  All  theses  are  or  will  be 
available  at  cost  from  the  Engineering  Projects  Laboratory,  Document 
Room  3-154,  M.I.T.,  Cambridge  39,  Massachusetts. 

The  first  few  chapters  emphasize  the  more  theoretical  aspects  of 
the  work,  and  the  last  few  the  more  practical.  Chapter  2  describes  an 
extensive  study  by  K.  N.  Reid  into  the  statics  and  dynamics  of  the  inter¬ 
action  of  a  fluid  jet  and  a  receiver  port.  A  partially -theoretical  partially  - 
experimental  approach  toward  static  diffuser  optimization  with  a  non-uniform 
inlet -velocity  profile  is  presented,  although  the  experimental  work  itself  ib 
as  yet  mostly  undone,  The  dynamics  of  fluid  lines  with  a  jet  incident  at  one 
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end  was  investigated  in  detail  both  analytically  and  experimentally,  with 
good  agreement.  It  is  shown  that  most  oscillations  can  be  quantitatively 
predicted  from  static  measurements  of  the  active  upstream  pressure-flow 
characteristics  at  the  jet-end  of  the  line,  plus  theoretical  considerations 
for  the  passive  part  of  the  line  and  associated  volumes,  restrictions,  etc. 
Viscous  wall  shear  is  shown  to  substantially  dampen  most  oscillations. 

(L.  E.  Johnston,  a  doctoral  fellowship  student  unsupported  by  this  contract, 
is  studying  the  problem  of  two  parallel  lines  coupled  by  a  jet  incident  at  one 
end. ) 

The  dynamics  of  the  controlled  deflection  of  jets  in  fluid  amplifiers 
and  relays  is  the  subject  of  Chapter  3,  prepared  by  R.  J.  Gurski.  Contin¬ 
uing  from  the  presentation  in  the  first  annual  report,  a  mathematical  model 
for  the  dynamics  of  the  jet  under  control  is  presented.  The  dynamic  control 
impedances  as  well  as  the  dynamic  response  of  a  system  is  the  goal  of  this 
as  -yet -incomplete  effort.  Considerable  use  of  digital  and  analog  computers 
is  involved. 

The  dynamic  analysis  of  systems  is  discussed  briefly  in  Chapter  4 
including  references  to  Engineering  Projects  Laboratory  publications  on 
fluid  lines  and  a  general  approach  toward  predicting  the  stability  of  fluid 
systems.  This  latter  development,  sponsored  by  the  present  contract, 
is  presented  in  detail  at  the  Harry  Diamond  Laboratories'  Second  Fluid 
Amplification  Symposium.  Also  in  Chapter  4  the  effect  of  bends  on  the 
dynamic  response  of  lines  is  discussed. 

In  Chapter  5,  S.  D.  Graber  presents  his  work  on  bistable  fluid- 
jet  relays  motivated  by  the  desirability  of  sensing  stall  cells  in  axial-flow 
compressors.  Background  work  was  presented  by  R.  J.  Gurski  in  le 
first  annual  report.  Graber  concentrates  on  the  static  stability  of  the  jet. 
Sensitivity  of  the  device  depends  on  a  small  stability  margin  for  the  deflected 
states  of  the  jet;  rapid  dynamic  response  depends  on  a  large  instability  for 
the  center  position  of  the  jet.  These  conditions  *un  counter  to  one  another 
but  an  optimum  configuration  exists.  Dynamic  experimental  results  reveal 
a  frequency  limitation. 

Chapter  6  by  A.  H.  Greenleaf,  presents  the  development  of  a  fluid 
oscillator  which  is  intended  to  be  efficient  and  nave  a  frequency  insensitive 
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to  large  operating  temperature  and  pressure  changes.  Mass -spring 
oscillators  of  several  types  were  considered,  and  a  circularly-oscillating 
wire -supported  disk  adopted.  Pressure  insensitivity  results  from  a  tan¬ 
gential  driving  force,  perpendicular  to  the  radial  spring  force.  Consequent 
dynamic  stability  proDlems  are  analyzed  and  then  configurationally  avoided. 
Temperature  sensitivity,  resulting  from  changes  in  Young's  modulus  for 
the  wire  springs,  is  compensated  by  an  axial  sub-buckling  force  induced  by 
a  temperature -sensitive  sealed  allows.  Although  the  final  breadboard 
model  worked  with  respect  to  its  significant  features,  further  minor  develop¬ 
ment  is  indicated  for  a  truly  viseful  oscillator. 

In  the  first  annual  report  a  fluid -diode -blade -wake  sensor  was 
proposed,  and  linear  analyses  and  initial  experimental  results  discussed. 

In  Chapter  7  of  the  present  report  the  conclusion  of  this  work  is  presented. 

A  nonlinear  dynamic  analysis  is  given.  Of  the  three  types  of  diodes  studied 
by.  W.  B.  Bails  only  the  vortex  diode  has  marginally  adequate  experimental 
static  characteristics,  but  unfortunately  this  diode  has  dynamic  characteristics 
which  fall  short  of  the  high-frequency  goal.  The  Tesla  diode,  patented  some 
40  years  ago,  was  given  brief  static  testing  anc.  found  to  be  superior  to  all 
the  others.  All  the  devices,  however,  are  shown  to  be  ineffective  at  low 
Reynolds  numbers,  greatly  limiting  small-scale  applications  such  as  the 
desired  sensor. 
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2.  JET  -RECIEVER  STUDIES 


2.1.  Background 

Various  aspects  of  the  over-all  problem  of  jet -receiver -diffuser 
steady-state  and  dynamic  performance  has  been  discussed  by  the  writer 
in  previous  reports  on  this  and  another  contract,  Ref.  1  and  2,  as  well 
as  in  the  writer's  M.I.T.  ScD.  Thesis  Proposal  (Ref.  3).  The  discussions 
of  static  jerformance  dealt  primarily  with  problem  definition  and  the 
adequacy  and  limitations  of  the  available  literature  in  providing  a  basis 
for  static  performance  prediction.  It  was  shown  that  the  over -all  problem 
of  static  performance  prediction  can  be  separated  conceptually  into  several 
related  fundamental  hydrodynamic  problems.  For  convenience,  these 
problems  were  classified  as  follows:  flow  through  "short  tube"  orifices  and 
nozzles,  flow  of  "controlled"  submerged  jets,  diffusion  of  nonuniform 
turbulent  streams  in  constant -area  passages,  and  conversion  of  kinetic 
to  potential  energy  in  "slow  expansion"  and  "sudden  expansion"  diffusers. 

In  the  previous  annual  report  on  the  present  contract^  Rc-f.  1), 
the  problem  of  jet-receiver  dynamic  interactions  was  discussed.  The 
potential  "causes  and  effects"  of  receiver  pressure  oscillations  were  briefly 
dJscu8sed.  A  simple  experiment  was  outlined  which  was  designed  to  study 
the  dynamic  behavior  of  the  interaction  region.  The  experiment  was  designed 
so  that  the  reflection  coefficient  of  the  source  terminal  (i.e. ,  the  interaction 
region)  could  be  measured.  No  significant  experimental  results  were  avail¬ 
able  at  the  time  of  writing. 

The  goal  of  the  over -all  program  of  study  of  jet-receiver  per¬ 
formance  is  to  provide  a  critical  portion  of  the  information  necessary  for 
rational  design  of  jet  modulator  systems.  Included  in  this  report  (Secs.  2.  2 
and  2.3)  are  analytical  and  experimental  results  obtained  to  date  on  certain 
of  the  jet-receiver  static  and  dynamic  problems  thought  to  be  particularly 
critical  in  terms  of  accomplishment  of  the  desired  goal. 
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2.2.  Receiver-Diffuser  Steady  -State  Performance 

2.2.1.  General  Characteristics  of  Diffusion  Processes 

The  underlying  problem  of  receiver -diffuser  design  is  the 
establishment  of  an  optimum  internal  geometry  based  on  a  desired 
criterion  of  merit.  For  example,  one  interesting  and  important 
"optimal"  problem  involves  maximizing  the  ratio  of  the  power  recovered 
in  the  receiver -diffuser  to  the  power  available  at  the  nozzle  supply  con¬ 
ditions.  Regardless  of  the  type  of  optimal  problem,  the  function  of  the 
receiver -diffuser  is  basically  the  same:  that  of  efficiently  transforming 
the  kinetic  energy  of  a  high  velocity  submerged  jet  into  controllable 
potential  energy. 

It  is  apparent  on  consideration  of  the  fundamental  character  of 
a  free  submerged  jet  that  the  receiver-diffuser  inlet  velocity  profile  will 
be  nonuniform  in  character.  Perhaps  the  most  important  single  problem 
affecting  steady-state  performance  prediction  is  a  complete  characterization 
of  the  receiver  inlet  velocity  profile.  Certainly  a  complete  description  of 
this  velocity  profile  must  take  into  account  the  critical  geometrical  and 
operational  parameters.  One  interesting  question  is:  how  far  back  (toward 
the  nozzle)  into  the  otherwise  free  submerged  jet  is  the  disturbance  intro¬ 
duced  by  the  receiver  actually  propagated?  Answers  to  these  questions  are 
goals  of  future  experimental  studies. 

Having  described  the  characteristics  of  the  flow  in  the  "free" 
region  between  the  nozzle  exit  and  the  receiver  inlet,  the  next  problem 
to  be  encountered  is  a  description  of  the  flow  within  the  receiver-diffuser. 
Diffusion  of  a  stream  in  internal  flow,  in  the  sense  of  increasing  static 
pressure,  is  accomplished  by  means  of  velocity  profile  "smoothing"  and/or 
simple  area  change.  It  is  well  known  that  diffusion  takes  place  in  a  diverging 
passage,  the  most  common  type  of  diffuser,  due  to  the  area  change  itself. 

It  is  often  forgotten  that  changes  in  the  velocity  profile  also  affect  the  dif¬ 
fusion  process,  i.e. ,  the  ability  of  a  passage  to  transform  energy  from 
kinetic  to  potential.  Equally  important,  but  less  well  known,  is  the  fact 
that  diffusion  may  take  place  in  a  constant -area  passage  due  to  the  effect 
of  velocity  profile  smoothing. 


In  order  to  understand  the  process  of  constant -area  diffusion, 
it  is  necessary  to  introduce  the  concept  of  uniformity  of  a  fluid  stream. 
Uniformity  refers  to  the  closeness  to  which  the  transverse  profile  of 
the  longitudinal  velocity  approaches  a  "rectangular’1  profile  (see  Fig.  1). 
A  convenient  quantitative  measure  of  the  uniformity  of  a  velocity  profile 
is  the  momentum  coefficient  defined  as  follows: 


true  momentum  of  the  stream 
a  =  - ■ — - ... - - — - — . —  — 

momentum  based  on  the  continuity  average  velocity 


Thus 
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Considering  now  the  constant-area  passage  shown  in  Fig.  2,  the 
law  of  conservation  of  momentum  requires  that 
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where 

Pj  =  static  pressure  at  station  1 

P^  =  static  pressure  at  station  2 

Vj  =  =  continuity  average  velocity 

a.j  and  a.,  =  momentum  coefficients 

r  =  wall  3hear  stress 
w 

P  =  fluid  mass  density 

Thus  the  difference  in  static  pressure  at  two  stations  is  a  function  of  the 
difference  in  the  velocity  profile  shapes  (given  by  a,  and  a <*)  at  the  two 
sections  and  the  wall  friction  forces  which  act  between  the  two  stations. 

It  is  well  known  that  the  velocity  profile  of  the  flow  in  a  constant -area 
passage  tends  to  stabilize  along  the  passage,  i.e.,  asymptotically  approach 
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Fig.  1.  Concept  of  Velocity  Profile  Uniformity 


Fig.  2.  Velocity  Profile  Smoothing  in  a  Constant -Area  Passage. 
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a  "fully-developed"  profile.  Wi^h  laminar  flow,  the  profile  always 
approaches  the  well-known  parabolic  profile  at  a  distance  sufficiently 
far  down  the  passage.  Likewise,  in  turbulent  flow,  the  profile  tends 
asymptotically  to  the  familiar  n~th  power  law  profile.  It  is  interesting 
that  the  value  of  a  for  a  1  /7 -th  power  profile  is  approximately  1.02. 

If  the  profile  at  a  st’tion  in  the  passage  is  highly  nonuniform 
(a  is  much  grtater  than  the  corresponding  fully  developed  value),  then 
Eq.  (2)  shows  that  an  increase  in  profile  uniformity  may  cause  an  in¬ 
crease  in  static  pressure  which  outweighs  the  decrease  due  to  wall 
shear.  Certainly  if  the  passage  is  long  enough,  wall  friction  will 
eventually  dominate.  It  should  be  noted  that  constant -area  diffusion 
is  accompanied  by  inherent  energy  loss,  as  in  any  diffusion  process. 

An  empirical  approach  to  the  study  of  constant -area  turbulent  diffusion 
is  discussed  in  Sec.  2.2.3. 

In  a  passage  whose  area  increases  in  the  direction  of  flow, 
the  "profile  smoothing"  form  of  diffusion  is  superimposed  on  the  area 
change  form.  Area -change  diffusion,  in  the  sense  of  increasing  static 
pressure,  ceases  when  the  area  change  effect  is  balanced  by  the  wall 
shear  effect.  Energy  losses  in  area-change  diffusers  primarily  result 
from  velocity  profile  smoothing,  wall  shear,  and  separation.  Experience 
indicates  that  a  nonuniform  velocity  profile  at  the  entrance  to  a  "slow 
expansion"  area -change  diffuser  leads  to  a  significantly  larger  over-all 
energy  loss  than  would  occur  with  a  more  uniform  inlet  profile  (see 
Sec.  2.2.4). 

On  the  basis  of  this  qualitative  understanding  of  diffusion 
processes,  past  experience  with  jet  pumps  and  ejectors,  and  con¬ 
siderations  of  practical  geometry  from  the  standpoint  of  construction, 
the  following  hypotheses  are  made  (Ref.  2): 

1.  There  exists  a  "best"  length  of  constant -area  passage 
together  with  a  "best"  geometry  of  the  increasing- 
area  passage  such  that  the  efficiency  of  the  diffusion 
process  is  maximized  for  a  given  set  of  conditions. 
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2.  There  exis*  a  "bes^'length  of  constant -area 
passage  together  with  a  ‘'best"  geometry  of  the 
increasing -area  passage  such  that  the  desired 
"shape"  of  the  static  pressure -flow  characteristic 
is  obtained  at  the  maximum  possible  efficiency. 

These  hypotheses  must,  of  course,  be  verified  by  suitable 
analysis  and  experiment. 

A  combined  analytical  and  experimental  program  is  underway 
which  is  aimed  at  defining  the  term  "best"  in  a  quantitative  way.  Pre¬ 
liminary  results  of  this  study  are  given  in  Secs.  2.2.3  and  2.2.4. 

2.2.2.  Constant -Area  Turbulent  Diffusion  in  Jet  Modulators 

In  a  sense,  the  problem  of  constant-area  turbulent  diffusion 
belongs  to  a  class  of  problems  known  as  "entry  length"  problems.  A 
brief  discussion  of  the  character  of  these  problems  is  given  in  Ref.  2. 
The  velocity  profile  development  of  a  turbulent  stream  in  the  constant- 
area  diffuser  of  a  jet  modulator  differs  in  the  following  respe  ..s  from 
that  encountered  in  the  usual  entry  length  problems: 

1.  The  velocity  profile  at  the  inlet  section  of  the  passage 

is  very  nonuniform  and.  in  fact,  may  be  nonsymmetrical 
as  well, 

2.  The  flow  is  already  "fully  developed"  turbulent  shear 
flow'  or  at  least  has  a  relatively  high  turbulence  intensity 
when  it  enters  the  passage, 

3.  Development  of  the  velocity  profile  in  the  direction  of 
flow  is  primarily  due  to  the  transfer  of  momentum 
between  fluid  layers  under  the  action  of  turbulent 
shear  stresses.  The  wall  shear  stress  is  small  in 
comparison  with  the  turbulent  tzresses,  at  least  for 
several  diameters  along  the  passage, 

4.  The  development  of  the  profile  is  accompanied  by  an 
initial  rise  in  static  pressure;  ultimately  the  point  is 
reached  where  wall  shear  becomes  important  and  the 
static  pressure  begins  to  decrease. 
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5.  Transition  from  an  unusually  nonuniform  profile  to 

one  having  the  3ame  degree  of  uniformity  as  the  corres¬ 
ponding  fully  developed  turbulent  pipe  flow,  takes  place 
in  only  a  few  diameters  rather  than  the  customary  20  -  50 
diameters  required  in  the  usual  entry  length  problem. 

Figure  3  serves  to  illustrate  schematically  what  might  be 
expected  to  occur  in  constant -area  diffusion.  In  Cases  a  and  b,  it 
is  assumed  that  the  flow  in  the  passage  is  turbulent  and  that  the  velocity 
profile  is  identical  to  the  profile  in  the  corresponding  section  of  an  un¬ 
obstructed  free  submerged  jet.  These  two  cases  represent  one  end  of 
the  spectrum  of  interest,  i.e.  ,  the  situation  which  should  exist  when  no 
throttling  takes  place  downstream  of  the  diffuser  (receiver  input  impedance 
is  essentially  zero).  When  considerable  throttling  takes  place  (receiver 
input  impedance  becomes  very  large),  the  flow  in  the  passage  may  he 
laminar  rather  than  turbulent.  Certainly  as  the  receiver  input  impedance 
approaches  infinity,  the  f’ow  in  the  passage  must  proceed  from  turbulent 
to  laminar.  Cases  c  and  d  of  Fig.  3  illustrate  cases  of  high,  but  not 
infinite,  receiver  input  impedance  (i.e.,  finite  through  flow).  In  all 
cases,  the  corresponding  fully  developed  profile  is  approached.  Very 
little  beyond  qualitative  reasoning  exists  to  substantiate  the  schematic 
representations  of  Fig.  3.  However,  the  results  of  a  single  experiment 
reported  in  a  1950  PhD.  thesis  by  H.  L.  Grirmrutt  (Ref.  4),  provides 
some  measure  of  verification  for  Case  a.  Grimmett  measured  the  velocity 
profile  in  a  circular  constant -area  duct  displaced  15  diameters  from  an 
axially -aligned  nozzle.  It  should  be  noted  that  there  was  no  throttling  of 
the  flow  leaving  the  downstream  end  of  the  duct,  and  consequently,  the 
static  pressure  at  the  exit  plane  was  atmospheric.  These  data  are  re¬ 
produced  in  Fig.  4.  Since  the  well  streamlined,  non -throttled  duct  places 
virtually  no  obstruction  in  the  field  of  the  submerged  jet,  the  duct  inlet 
velocity  profile  is  indistinguishable  from  the  corresponding  portion  of  the 
free  jet  profile.  It  is  evident  that  the  transition  from  the  higliy  nonuniform 
inlet  profile  to  a  seemingly  fully  developed  profile  takes  place  in  4-5  diameter 
It  is  interesting  that  the  point  at  which  the  vvall  static  pressure  first  peaks 
occurs  somewhat  before  the  apparent  point  at  which  the  profile  is  fully 
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Fig.  4.  Velocity  Profile  Development  in  a  Circular  Constant-Are 
Displaced  15  Diameters  from  an  Axially-Aligned  Nozzle. 
(From  data  tabulated  in  Ref.  4) 


developed.  Furthermore,  the  static  pressure  data  show  the  relative 
urt'mportance  of  the  total  integrated  wall  shear  stress  in  the  first  4  or 
5  diameters  of  the  duct.  The  Important  conclusion  is  the  fact  that  the 
velocity  profile  smoothing,  under  the  action  of  turbulent  shear  stresses, 
takes  place  much  more  rapidly  than  the  profile  "shaping"  as  a  result  of 
wall  boundary  layer  growth.  These  data,  although  admittedly  much  too 
limited  in  scope,  provide  the  basis  for  the  previous  hypotheses  concerning 
the  length  of  duct  necessary  to  achieve  a  profile  of  high  uniformity. 

2.2.3.  Empirical  Approach  to  the  Study  of  Constant -Area 
Turbulent  Diffusion 

It  is  entirely  possible  that  jet  modulator  static  performance 
can  be  predicted  adequately  by  simplified  means  without  ever  really 
being  concerned  with  the  details  of  t  e  diffusion  processes.  However, 
even  if  this  ideal  could  be  achieved,  it  is  unlikely  that  there  would  exist 
a  truly  rational  basis  for  optimization  studies.  Consequently,  one  seems 
justified  in  considering,  insofar  as  is  feasible,  the  details  of  the  flow 
processes . 

In  establishing  an  analytical  basis  for  optimizing  the  over -all 
diffusion  process,  it  is  necessary  to  have  a  knowledge  of  the  velocity 
distribution  in  the  diffuser  passages.  On  seeking  a  solution  to  the  constant 
area  turbulent  diffusion  process,  one  is  faced  with  the  inadequacy  of  our 
understanding  of  even  the  most  simple  turbulent  flow  processes. 

Truly  no  "exact"  analytical  solution  to  the  basic  equations  of 
motion  has  ever  been  found  for  a  turbulent  flow  problem.  Nevertheless, 
the  well-known  semi-empirical  hypotheses  established  oy  Prandtl, 
Reichardt,  and  Taylor  have  been  applied  with  a  reasonable  degree  of 
success  in  a  number  of  important  engineering  problems  and  present 
some  hope  for  the  present  problem.  Of  course,  a  solution  to  the 
turbulent  diffusion  problem  would  not  be  complete  by  itself.  The  truly 
useful  solution  must  be  applicable  regardless  of  flow  regime  (i.  e. , 
turbulent,  transition,  or  laminar}.  It  seems  unlikely  that  a  single 
solution  can  be  found  that  will  be  universally  applicable.  What  is  more 
probable  is  that  two  or  more  solutions  can  be  "matched"  judiciously! 
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The  turbulent  diffusion  problem  can  now  be  stated  as  follows: 

Given  the  initial  distribution  of  velocity  and  static  pressure  at  a  certain 
passage  section  (i.  e.,  at  the  inlet  section  of  the  passage),  what  is  velocity 
profile  at  another  section  downstream?  Certainly  a  hint  for  a  solution  is 
offered  by  experimental  evidence  (Refs.  4  and  5)  that  the  turbulent  mixing 
process  is  far  more  rapid  than  the  boundary  layer  growth.  The  success  of 
Reicharcl.'s  inductive  theory  of  turbulence  when  applied  to  submerged  jets 
and  other  problems  which  have  a  principal  flow  direction,  suggests  that 
his  theory  may  be  applicable  to  the  present  problem.  This  theory  has 
been  applied  with  apparent  success  to  problems  of  free  jets,  jets  dis¬ 
charging  parallel  to  a  wall,  and  jets  discharging  at  the  center  of  a  duct, 

-f- 

by  a  group  of  investigators  at  the  Engineering  Experiment  Station  of  the 
University  of  Illinois  (Refs.  6  and  7).  The  technique  developed  to  handle 
the  case  of  a  jet  issuing  at  the  center  of  the  duct,  as  in  a  jet  pump,  is 
empirical  n  nature  in  that  two  constants  must  be  determined  experimentally. 
Such  is  typical  for  turbulent  flow  problems.  An  empirical  expression  is 
presented  in  Ref.  7  for  the  prediction  of  the  transverse  total  momentum  dis¬ 
tribution  at  the  various  sections  of  the  duct.  Recently,  R.  Moissis  and 
P.  Griffith  of  M.I.  T.  successfully  utilized  the  technique  developed  by  the 
previous  investigators  for  the  "ducted  jet",  to  predict  the  total  momentum 
distribution  in  the  wake  of  an  air  bubble  in  developing  two-phase  slug  flow. 
(Ref.  8).  The  same  general  approach  will  be  followed  in  the  analysis  below. 

Consider  now  the  steady,  incompressible  axially -symmetric 
turbulent  flow  in  a  constant -area  circular  passage.  (The  results  of  an 
analysis  for  a  two-dimensional  passage  will  follow  the  present  development.) 
Boundary  layer  development  at  the  walls  is  neglected.  The  theoretical 
reasoning  of  Reichardt's  approach  is  discussed  in  Refs.  9  and  10.  The  theory 
is  baaed  on  the  momentum  equation  for  the  time -mean  averages  of  the  velocity 
in  the  principal  direction  of  flow.  Thus 


+  The  PhD.  thesis  by  H.  L.  Grimmett  mentioned  in  Sec.  2.2.2.  is  part  of 
this  larger  study. 
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where 

p  =  fluid  mass  density 

u  =  instantaneous  velocity  in  the  x-direction 

u  =  time -mean-average  velocity  in  the  x-direction 

v  =  instantaneous  velocity  in  the  r -direction 

v  -  time -mean-average  velocity  in  the  r~direction 

x  =  distance  along  axis  of  cylindrical  passage 

r  =  radial  position  in  cylindrical  passage  measured 
from  axis  of  passage 

P  =  static  pressure;  P  =  mean  static  pressure 


The  second  basic  equation  in  Reichardt's  theory  is  of  an  empirical  nature 
and  has  the  form 


-  A  3(u  ) 
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The  term  uv  represents  the  turbulent  shearing  Rtrees.  The  parameter  A 
is  called  the  "momentum  transfer  length.  "  It  has  the  dimensions  of  length 
and  is  assumed  to  be  a  function  of  x  only.  Equation  (4)  thus  relates  the 
turbulent  shear  stress  to  an  easily  measured  mean  flow  quantity  by  means 
of  a  "hypothetical"  transfer  length.  In  words,  Eq.  (4)  states  that  the  flux 
of  x-component  of  momentum  which  is  transferred  in  a  transverse  direction 
is  proportional  to  the  transverse  gradient  of  momentum. 

Combining  Eqs.  (3)  and  {4)  and  rearranging  yields 


l£.  +  _3_  (pu2)  .-1  A  i- 
9x  9x  r  9  r 


0 


(5) 


Assuming  the  pressure  gradient  in  the  transverse  direction  (i.e.,  9p/9r)  is 
of  small  order  compared  to  the  longitudinal  gradient  9p/9x,  dp/dr  can  be  added 
to  both  sides  of  Eq.  (5).without  destroying  its  validity.  Thus 
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3  .  -  ,  2.  1  A  3 

-  (p  +p  u  ) - A  - 

8x  r  dr 


3  /-  2, 

r  -  (p  +  pu  ) 

3r 


Defining  m  =  p  +  p  u  gives 


J-A-i 


i=J  =  o 

3r 


In  order  to  simplify  the  boundary  conditions,  it  is  convenient  to  make 
the  following  definitions:  r+  =  2r/D,  where  D  =  passage  diameter 
and 

M  =  (m  -  mJ/mQ 


where 


m^  =  (p+pu  )  at  x— <*-oo;  taken  as  the  fully  developed 

turbulent  flow  total  momentum 
_  - £ 

mQ  =  (p  +  pu  )  at  x  =  0;  x  =  Ois  the  total  momentum  at  the 
entrance  of  the  passage  and  is  assumed  to  be  known. 


Using  these  definitions,  je.q.  (7)  can  be  written  in  dimensionless  form  as 


follows : 


+  3M  n 
r  -  -  0 


A  3 

D2  3r  + 


This  is  a  linear,  second-order  partial  differential  equation  in  the  variable  M. 
Four  independent  boundary  conditions  can  be  specified: 

1.  {3M/3r  +  )r  +  _  q  =  0,  from  symmetry. 

2.  (3M/3r+)r+_1  =  0,  since  wall  friction  is  neglected. 


3.  (M)x  ^  =  i.e. ,  fully  developed  value  of  M  is  a 

useful  reference.  It  can  be  calculated  from  an  accepted 
power  law  profile. 

4.  <M)X=0  =(m0-m«ymo  '  1  '  mco/mo  =  f(r+,; 
this  value  of  M  is  assumed  to  be  known. 
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Solution  of  Eq.  (8)  by  separation  of  variables  and  subsequent 
substitution  of  the  first  two  boundary  conditions  above  yields  the  product 
solution. 


where 


A  J  (a  r  )  exp 
non  ^ 


f  4  A 


O' 


D" 


J  =  Bessel  function  of  zero-th  order 
o 

=  Bessel  function  of  first  order 

a  =  n-th  root  of  J, 
n  i 

A  =  a  constant 
n 


(9) 


The  function  A  =  A  (x)  which  appears  in  the  exponential  (i.  e.  ,  the 
momentum  transfer  length)  must  be  determined  empirically.  Given  a  set 
of  experimental  values  of  M,  the  function  A  (x)  can  be  defined  in  such  a 
manner  that  experimental  and  theoretical  values  of  M  are  identical.  Re¬ 
calling  the  definition  of  M,  it  is  obvious  that  measurements  of  velocity 
distribution  and  static  pressure  at  various  stations  along  the  passage  will 
yield  the  desired  experimental  values  of  M.  Equipment  has  been  designed 
and  is  now  being  constructed  for  the  purpose  of  making  measurements  ox 
this  type.  A  pitot  probe  will  be  used  for  these  measurements.  Additional 
measurements  will  be  made,  using  X-array  hot-wire  probes,  of  the  actual 
turbulent  shear  stress  distribution  at  various  passage  stations.  These 
latter  measurements  will  provide  the  additional  information  necessary  for 
direct  integration  of  the  original  equation  of  motion  (i.  e. ,  Eq.  (3)).  The 
validity  of  the  simplifying  assumptions  made  in  establishing  Eq.  (3)  and, 
indeed,  the  whole  of  the  Reichardt  theory,  can  then  be  determined. 

It  would  be  highly  desirable  to  establish  an  explicit  relation  for 
A  =  A  (x).  The  assumption  has  already  been  made  that  A  is  independent  of 
This  turns  out  to  be  a  good  assumption  when  the  Reichardt  hypothesis  is  used 
for  free  jet  analyses,  in  contrast  to  free  jets,  A  cannot  grow  indefinitely 
(as  x  increases)  in  a  constant  diameter  passage.  In  the  absence  of  wall 
friction,  it  seems  reasonable  that  the  presence  of  a  wall  should  have  no 
other  effect  on  A  than  to  place  an  upper  limit  on  its  growth.  A  purely 

-  18  - 


i 


■jggpgwp 


arbitrary  function,  supsjest.ed  in  Ref.  7.  that  reduces  to  the  free  jet 
relation  for  small  x,  but  which  imposes  on  A  an  upper  limit  which  is  a 
function  of  the  passage  diameter  is 


2 

c  x 


(10) 


where  c  and  b  are  arbitrary  constants.  By  having  such  an  arbitrary 
function,  it  would  seem  that  one  has  a  lot  of  latitude  in  correlating  ex¬ 
perimental  data.  It  is  to  be  expected  that  the  form  of  Eq.  (10)  will  be 
reasonable  for  a  passage  with  little  downstream  throttling.  In  this  case, 
the  value  of  A  at  the  entrance  of  the  passage  should  be  very  nearly  identical 
with  that  of  the  corresponding  free  jet.  Whether  Eq.  (10)  is  of  a  reasonable 
form  when  there  is  considerable  downstream  throttling  -emains  to  be 
established. 

Assuming  the  utility  of  Eq.  (10),  it  is  possible  to  carry  the 
analysis  further.  Combining  Eqs.  1 9)  and  (1C)  yields 


=  2 


A  J  (a  r  )  exp^  - 
n  o'  n  M 


f  ,  2  2  _ 

i  2c  an  F bx 


b2  LD 


In  (1  +  ±x_  ( 


The  remaining  two  boundary  conditions  can  now  be  utilized.  These  are 


and 


M  — 0  as  x — *-co 

M  =  f  (r+  ). 
o 


The  first  condition  is  satisfied  by  the  exponential  term  of  Eq.  (11)  which 
tends  toward  zero  as  x  increases  without  limit.  The  second  condition  yields 


oo 

f(r+)  =  ^  J  (a  r+) 

'  '  ,  no  '  n 

n=  1 


(12) 
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By  virtue  of  the  orthogonality  of  the  Bessel  series 


A 

n 


2 


J  ' 
o 


1 

f(rT)  Jq  Unr  +  )  r+  dr  + 

*0 


(13) 


provided  the  series  is  uniformly  convergent.  (Ref.  11)  Now  any  arbitrary 
initial  profile  of  the  total  momentum  function,  Mq,  can  be  imposed  on  A^. 
Except  for  the  case  of  little  downstream  throttling,  a  rational  means  for 
evaluating  Mq  under  all  downstream  load  conditions  does  not  exist  presently. 
The  earlier  mentioned  experimental  program  is  aimed  at  providing  this 
?  '•'formation. 


There  is  as  yet  no  known  way  to  predict  the  values  of  the  empirical 
constants,  b  and  c,  in  F.q.  (11).  Resource  must  be  made  to  judicious  ex¬ 
periment.  However,  point  by  point  application  of  Eq.  (11)  is  not  a  satis¬ 
factory  means  of  evaluatin'  b  and  c  from  experimental  data.  It  is  possible, 
however,  to  greatly  simplify  the  evaluation  process  by  defining  a  one- 
dimen3ional  function  which  is  integrable  in  r+ .  The  integral  of  M  itself 

is  not  a  suitable  function  since,  by  conservation  of  momentum,  it  is  always 

2 

zero.  However,  the  integral  of  M  is  not  conserved,  and  may  be  used  to 

define  a  one -dimensional  function,  I,  as  follows: 

•> 

l 

I'--  (  M2  x  2Trr+  dr+  (14) 

0  J 


The  function  I  appropriately  may  be  termed  the  ’’mixing  index." 

For  any  set  of  values  of  b  and  c,  an  equation  resulting  from  a 
combination  of  Eqs.  (11),  (13),  and  (i4)  yields  an  analytical  expression  for 
the  variation  of  I  as  a  function  >_■;  ;t  only.  In  addition,  the  value  of  I  corres¬ 
ponding  to  any  position  x  may  be  found  by  performing  graphically  the  inte¬ 
gration  of  Eq.  \14).  In  the  latter  case,  the  values  of  M  are  obtained  from 
experimental  data.  By  trial-and-error  the  values  of  b  and  c  can  be  chosen 
such  that  the  "analytical"  and  "measured"  values  of  J.  show'  good  agreement. 
The  success  of  this  procedure  awaits  the  results  of  the  experimental  program} 
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Consider  now  the  case  of  constant -area  turbulent  diffusion  in  a 
two-dimensional  passage.  The  momentum  equation  for  the  time-mean- 
a%rerages  of  the  velocity  in  the  principal  direction  of  flow'  becomes  for  this 
case 

Z  - 

3u  3uv  9p 

p - +  p  -  =  -  — — 

3x  3y  3x 

where  x  is  the  coordinate  in  the  principal  flow  direction  and  y  is  the 
coordinate  in  the  tr*  asverse  direction  measured  from  the  axis  of  flow. 

Again,  using  Reichardt's  hypothesis  for  the  turbulent  shear 
stress  leads  to  the  following  partial  differential  equation 

3M  _  4  A  32M 

3x  2  a  +2 

h  By 

where  his  the  minimum  transverse  dimension  of  the  passage  and  y  =  y/h. 
The  boundary  conditions  are  of  the  same  form  as  for  the  circular  passage. 
Again  making  use  of  the  separation  of  variables  technique  for  solving  partial 
differential  equations,  the  following  product  solution  for  M  is  obtained 


From  this  point  on,  the  same  considerations  hold  as  in  the  circular 
passage  case. 

2.2.4.  Area -Change  Diffuser  Characteristics 

The  theoretical  prediction  of  flow  conditions  in  area -change 
diffusers  has  long  been  an  unsolved  problem  in  fluid  mechanics.  Re¬ 
searchers  for  over  50  years  have  relied  almost  totally  on  experimental 
techniques  in  the  design  of  diffusers.  Two  basic  problems  have  been 
stud'ed  extensively  in  connection  with  the  design  of  straight -walled  diffusers. 

1.  Prediction  of  optimum  geometry  for  most  efficient 
transformation  of  kinetic  to  potential  energy,  and 
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flow  regimes  encountered  with  various  degrees  of 
separation. 

Many  basic  problems  remain  unanswered,  however.  Nevertheless,  con¬ 
siderable  experimental  cat'1  are  available  on  the  "efficiency"  of  circular, 
square,  and  rectangular  diffusers  having  various  area  ratios,  divergence 
angles,  inlet  boundary  layer  conditions,  and  Reynolds  numbers.  A  brief 
bibliographical  review  of  the  literature  is  given  by  the  writer  in  Ref.  2. 

At  the  outset  of  the  receiver-diffuser  study,  it  was  believed  that 
considerable  data  existed  which  could  be  used  directly  in  predicting  per¬ 
formance  of  jet-modulator  area-change  diffusers.  The  results  of  a  fairly 
comprehensive  literature  survey  showed  that  considerable  information  is 
lacking.  It  is  diffucult  to  make  direct  comparisons  of  the  available  results 
of  isolated  studies  since  the  experimental  test  conditions  are  not  always 
well  <^fined  or  even  controlled.  In  addition,  there  is  no  common  agree¬ 
ment  in  the  literature  on  the  most  convenient  definition  of  efficiency.  Con¬ 
sequently,  experimental  results  giving  "efficiency"  as  a  function  of  diffuser 
geometry  and  various  operating  conditions,  is  often  extremely  co”  fusing. 
Even  more  important  is  the  fact  that,  almost  without  exception,  available 

diffuser  performance  data  are  restricted  to  relatively  high  diameter 

4 

Reynolds  numbers  (say  10  and  up).  It  is  likely  that  in  most  future 

4 

applications  of  jet  modulators,  Reynolds  numbers  between  1  and  10 
will  be  encountered. 

It  is,  of  course,  fortuitous  that  so  much  experimental  data  exist 
for  the  high  Reynolds  number  range.  The  important  trends  and  techniques 
established  by  the  previous  studies  will  be  invaluable  in  the  efficient  planning 
of  future  experiments.  Indeed,  ..he  hypotheses  concerning  optimum  diffi’ser 
geometry  discussed  in  Sec.  2.  2.  2.  were  based  almost  totally  on  available 
diffuser  performance  data.  Since  the  available  data  are  so  important  to 
the  over-all  goal  of  the  static  performance  study,  some  of  thebest  avail- 
ble  pertinent  diffuser  data  are  reviewed  below.  Only  the  results  of  studies 
with  conical  diffusers  will  be  discussed  below.  The  interested  reader  is 
referred  to  Refs.  12,  13,  and  14  for  extensive  discussions  of  diffusers  of 
all  shapes  and  forms. 
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Oi  the  various  efficiencies  used  in  the  literature,  the  most 
common  and  probably  the  most  useful  in  the  majority  of  situations  is 
the  "static  pressure  efficiency",  or  as  sometimes  termed,  the  "energy- 
conversion  efficiency".  This  efficiency  is  defined  as  the  ratio  of  the 
rate  at  which  energy  is  transformed  to  the  rate  at  which  energy  is  supplied 
for  transformation.  Consider  a  passage  which  expands  from  an  initial 
area  to  a  final  area  A^.  Let  P  be  the  static  pressure,  p  the  fluid 
mass  density,  u  the  component  of  velocity  parallel  to  the  axis  of  the 
passage,  and  v/  the  absolute  velocity  at  any  point.  Thus  the  actual  energy- 
conversion  or  static  pressure  efficiency  is 


^  P  udA  -  ^A}  P  udA 

f  1  2  r  i 

\  —  uw  udA  -  \  —  p 


\  pw2  udA 

A1  2 


J  ? 


p  w  udA 


If  the  flow  is  purely  axial  at  sections  1  and  2,  w  =  u  and  p  is  a  constant 
over  the  section  (i.  e.,  no  streamline  curvature).  Then 


(p2  -  pp  a,  Aj 


where 


f  1  3  Cl 
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Equation  (16)  can  be  simplified  by  defining  a  kinetic  energy  factor  3 


as  follows : 


u  dA 


1  -3  A 

—  p  u  A 

2 
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The  kinetic  energy  factor  is  analogous  to  (but  not  equal  to)  the  momentum 
factor  a,  defined  in  Section  2.2.  3,  in  that  it  accounts  for  the  nonuniformity 
of  the  velocity  profile.  Utilizing  this  definition  and  recognizing  that 
continuity  requires  u^A^  =  the  following  expression  is  obtained: 


(19) 


Many  investigators  have  used  a  one -dimensional  approximation  to  Eq.  (19) 
as  follows: 


(20) 


Equation  (19)  reduces  to  Eq.  (20)  if  it  is  assumed  that  the  velocity  profiles 
at  the  entrance  and  exit  sections  are  uniform  (i.  e„ ,  Pj  =  &2  = 

The  efficiency  depends  on  many  factors,  the  moat  important  of 
which  is  the  rate  of  expansion  of  the  flow.  This  is  not  surprising  since  it 
is  this  factor  which  determines  the  magnitude  of  the  opposing  pressure 
gradient.  Some  results  of  experiments  on  the  flow  between  straight 
diverging  walls  (conical  section  only)  by  Gibson  (Ref.  15)  and  Peters  (Ref.  16) 
are  summarized  in  Fig.  5.  Efficiency,  r^,  is  plott<  d  against  20,  the 
included  angle  between  the  diverging  walls.  The  efficiency  is  low  at  very 
small  angles,  rises  to  a  peak  in  ♦’he  vicinity  of  5  -  8  degrees,  rapidly  falls 
a3  20  increases  until  a  low  point  is  reached  in  the  vicinity  of  60  degrees,  and 
then  slowly  approaches  the  theoretical  value  given  by  Borda  at  20  =  180°. 

The  Borda  limiting  condition  (i.  e.  ,  a  sudden  expansion  diffuser)  is  derived 
from  simple  continuity  and  momentum  considerations;  the  following  ex¬ 
pression  results: 
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^  J  sudden  A9 

j  expansion  1  +  - — 

A! 

Perhaps  the  most  important  feature  of  the  results  of  Fig.  5  is  that  the 
efficiency  reaches  a  peak  value,  which  is  quite  high,  in  the  vicinity  of 
5-8  degrees. 

To  understand  the  reason  for  the  peak  in  the  efficiency  curve, 
consider  a  family  of  diffusers  having  identical  area  ratios  (A^/A^  =  const.). 
If  the  included  angle  is  small  (say  <  5°),  the  wall  friction  losses  are  large 
since  the  passage  must  be  very  long.  However,  the  "separation"  losses  are 
extremely  small.  Conversely,  if  the  included  angle  is  large  (sav  >  15°), 
the  passage  length,  and  therefore  the  wall  shear  losses,  are  greatly  mini¬ 
mized.  The  separation  losses,  though,  are  greatly  increased  as  the  in¬ 
cluded  angle  is  increased.  The  peak  efficiency  condition  then  represents  an 
optimum  trade-off  between  -wall  shear  and  separation  losses. 

A  factor  which  is  often  obscured  in  the  presentation  of  diffuser 
performance  data  is  the  effect  of  inlet  and  exit  velocity  profile  nonuniformity. 
Although  efficiency  calculations  based  on  Eq.  (19)  take  into  account  profile 
nonuniformity  directly,  most  investigators  have  not  bothered  to  make  the 
additional  (and  tedious)  measurements  necessary  to  establish  (3^  and  (3^. 

The  approach  used  by  most  investigators  (if  they  bothered  at  all)  has  been 
to  study  the  effect  of  various  inlet  velocity  profiles  produced  by  using  dif¬ 
ferent  lengths  of  constant  diameter  inlet  pipe.  For  a  very  short  inlet  pipe, 
the  value  of  p^  is  approximately  equal  to  1.  As  the  inlet  pipe  length  is 
increased  to  say  30  diameters,  the  value  of  p^  is  very  nearly  that  of  a  fully 
developed  turbulent  flow  profile  (p^  -  1.05,  dependent  on  Reynolds  nurr1.' 

The  percentage  difference  between  the  efficiencies  calculated  from 
Eqs.  (19)  and  (20)  may  be  as  small  as  a  few  per  cent  or  as  large  as  10  -  15 
per  cent.  Even  if  the  diffuser  entrance  profile  is  approximately  uniform 
(i.e.,  Pj  =  1.0),  the  exit  profile  normally  is  quite  nonuniform  (P2  >1.0). 
Figure  6  is  illustrative  of  typical  flow  conditions  associated  with  a  conical 
diffuser  having  "near  optimum"  geometry. 
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Due  to  the  fact  that  the  diffuser  outlet  velocity  profile  is 
typically  quite  nonuniform,  the  addition  of  a  constant  diameter  mixing 
section  downstream  of  the  exit  section  may  have  an  important  effect  on 
the  diffuser  efficiency.  A  diffuser  with  constant  diameter  inlet  and  outlet 
sections  is  shown  in  Fig.  7a.  If  the  velocity  profile  at  the  diffuser  exit 
is  more  nonuniform  than  a  fully  developed  profile  (and  in  most  cases  it  is), 
the  profile  becomes  progressively  smoother  in  the  discharge  section  which 
follows.  As  a  result,  the  static  pressure  continues  to  rise  in  the  discharge 
section  until  a  position  is  reached  at  which  the  static  pressure  is  a  maximum. 
From  that  point  on,  the  pressure  decreases  due  to  wall  shear.  Peters 
(Ref.  16)  has  found  that  the  velocity  profile  at  the  point  of  maximum  static 
pressure  is  very  nearly  the  same  as  the  fully-developed  profile.  Results 
obtained  by  Peters  showing  the  point  of  maximum  static  pressure  recovery 
as  a  function  of  diffuser  included  angle  and  entrance  section  length  are  given 
in  Fig.  7b.  The  importance  of  the  discharge  length  on  diffuser  efficiency 
is  obvious.  In  fact,  modified  definitions  of  efficiency  are  suggested  which 
utilize  the  static  presstire  at  the  point  of  maximum  recovery  rather  than 
at  tne  exit  plane  of  the  area  change  section.  A  one -dimensional  relation 
which  corresponds  to  Eq.  (20)  is 


(22) 


Similarly  can  06  substituted  into  Eq.  (19)  for  the  true  efficiency  in 

place  of  P~.  Therefore,  the  difference  between  r]  and  rj  *  is  a  measure 

P  P 

of  the  effect  of  the  addition  of  adequate  length  of  discharge  section. 

The  difference  between  the  true  or  actual  efficiency  between 
two  sections  and  the  corresponding  one -dimensional  approximation  is  a 
measure  o:  the  nearness  to  unity  of  the  values  of  (3  at  the  two  sections. 
Figure  8  shows  an  interesting  comparison  between  the  actual  efficiency 
between  sections(l)  and  (2),  the  approximate  efficiency  between  (1)  and  (2), 
and  the  approximate  efficiency  between  (1)  and  (2').  The  actual  efficiency 
between  sections  (1)  and  (2')  is  not  shown  but  could  be  inferred  quite  easily. 
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in  the  inlet  and  exit  velocity  profiles.  What  is  most  important  about  these 
results  is  the  significance  of  even  a  seemingly  uniform  profile  (i.e.,  fully 
developed)  on  the  value  of  the  efficiency.  If  the  submerged  jet  of  a  jet 
modulator  were  to  be  diffused  directly  by  an  area  change  diffuser,  the 
corresponding  efficiency  would  be  quite  low  (even  for  so-called  optimum 
geometries)  because  of  the  high  nonuniformity  of  the  inlet  velocity  profile. 

It  was  this  belief  that  led  to  the  hypotheses  of  Sec.  2.2.2.  A  careful  study 
of  the  available  literature  and  past  experience  with  jet-pipe  valves  leads 
the  writer  to  believe  that: 


1.  the  most  efficient  receiver-diffuser  for  a  specified 
''degree  of  diffusion"  (i.e.,  total  area  change)  is  one 
having  in  succession  a  constant -area  mixing  sectio.  , 
a  "short"  .  -  8  degree  diverging  section,  and  finally 
a  sudden  expansion  section. 

2.  the  most  linear  pressure -flow  characteristics  of  a 
three  way,  axially  symmetric,  axially  aligned  jet- 
receiver-diffuser  will  result  with  a  short,  constant  - 
diameter  mixing  section  followed  by  a  sudden  expansion 
section. 

The  discussion  above  gives  only  a  sample  of  the  factors  which 
must  be  considered  in  the  efficient  design  of  area-change  diffusers.  Much 
remains  to  be  done  in  properly  evaluating  and  extending  the  existing  diffuser 
performance  data.  A  more  detailed  discussion  of  diffuser  performance 
will  be  presented  in  the  writer's  ScD.  thesis  (summer  1964). 

2.  2.  5.  Orifice  Characteristics 

In  many  applications  of  jet  modulators  (i.e.,  a  "four-way"  jet -pipe 
valve  or  jet  amplifier),  fluid  may  flow  out  of  the  receiver -dif^jser  into  a 
reservoir  or  plenum  under  certain  conditions  of  downstream  loading.  In 
such  situations,  the  receiver -diffuser  flow  is  "reversed",  so  to  speak,  and 
the  fiow  characteristics  would  be  expected  to  be  radically  different.  Indeed 
they  are.  If  the  receiver -diffuser  consists  of  a  constar.t-diameter  mixing 
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section  followed  by  an  area-change  section,  then  the  reversed  flow 
situation  is  that  of  flow  through  a  "short -tube"  orifice  with  a  converging 
or  "sudden  contraction"  entry  section.  Curiously  enough,  the  flow 
characteristics  of  orifices  of  this  type  are  not  very  well  established, 
especially  in  the  low  Reynolds  number  range.  A  cursory  discussion  of 
the  available  literature  on  short-tube  orifices  is  gi/en  by  the  writer  in 
Ref.  3. 

The  writer  reported  the  results  of  a  recent  M.I.T.  study  or 
sharp-edged  short -tube  hydraulic  orifices  in  Ref.  17.  This  study  was 
carried  out  by  D.  J.  Tapparo,  an  S.M.  thesis  student,  under  the  super¬ 
vision  of  the  writer.  The  results  of  the  study  have  ledto  an  empirical 
means  of  predicting  the  discharge  coefficient  for  sb  rp-edged  short -tube 
orifices  within  the  following  limitations: 

2  <  —  <  6,  ICO  <  Red  <  1000 

and  4000  -  Re^  -  25,  000  (where  L  =  orifice  length,  D  =  orifice  diameter, 
and  Re^  •=  diameter  Reynolds  number).  It  is  expected  that  these  results, 
even  though  somewhat  limited,  will  be  extremely  useful  to  the  writer  in 
e  tablishing  a  means  for  predicting  over -all  static  performance  of 
receiver -diffusers. 

Also  expected  to  be  c.  considerable  value  a:e  the  results  of  a 
study  of  the  friction  factor  in  the  laminar  entry  region  of  a  smooth  tube 
by  Shapiro  et  al  several  years  ago.  These  results  are  valid  for  a  short  - 
tube  with  a  bell-mouth  and  for  Re  <5  x  10^  (where  Re  is  the  Reynolds 
number  basea  on  length  along  the  tube).  With  rather  easily  defended 
assumptions,  the  orifice  discharge  coefficient  car.  be  shown  to  be  a  function 
of  the  integrated  apparent  friction  factor  multiplied  by  tne  L/D  ratio. 
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2.  2.  6.  Experimental  Study  of  an  Axially -Aligned,  Axi -Symmetric , 
Three -Way  let  Ivfoavlator  ’  ' 

It  is  anticipated  that  results  of  studies  discussed  in  previous 
sections  will  eventually  yield  techniques  useful  in  predicting  over -a  11 
static  performance  of  jet  modulators  based  on  empirically  derived 
characteristics  of  sub-parts.  Consequently,  an  experimental  study 
has  been  initiated  to  provide  a  means  for  evaluating  these  techniques. 

A  large-scale  model  of  an  axi-symmetric,  axially-aligned  th~ee-way 
jet  modulator  has  been  designed  and  fabricated.  Freliminat-y  .ests 
have  been  conducted  at  relatively  high  Reynolds  numbers  and  with  a 
few  geometrical  variations.  These  results  are  presented  in  an  S.B. 
thesis  by  H.  Leibowitz  (Ref.  18).  The  apparatus  was  designed  so  that 
a  systematic  study  can  be  made  of  the  over -all  diffusion  process  for  a 
wide  number  of  internal  and  external  receiver -diffuser  geometries. 
Provisions  have  been  made  in  the  apparatus  for  varying  the  nozzle -to - 
receiver  spacing,  the  length  of  constant -area  mixing  section,  and  the 
incluued  angle  and  area  ratio  of  the  area  change  diffuser  section  (by 
using  a  series  of  molded  plastic  diffusers).  A  schematic  of  the  ex¬ 
perimental  apparatus  is  shown  in  Fig.  9.  The  apparatus  is  sufficiently 
flexible  so  that  radically  different  geometries  could  be  incorporated 
if  desired. 

Illustrative  of  the  preliminary  experimental  results  are  the 
pressure-flow  data  shown  in  Fig.  10.  These  date  low  the  effect  of 
variation  of  included  angle  of  the  area  change  diffuser.  The  significant 
increase  of  maximum  power  recovery  (Q  x  P  )  by  use  of  a  low  angle 
diffuser  is  f  *ident. 

A  comprehensive  experimental  study  is  planned  which  will 
cover  a  wide  Reynolds  number  range  and  incorporate  a  large  number  of 
geometrical  variations.  Air  (or  water)  and  hydraulic  oil  will  be  used  as 
working  fluids.  Results  will  be  presented  when  they  become  available. 

2.3.  Jet-Receiver  Dynamic  Interaction  Studies 

2.3.1.  Preliminary  Measurements  of  Source  Dynamic 
Characte  rlstlc  sT 

In  the  previous  annual  report  on  this  contract  (Reff  1),  an 
experiment  involving  a  simple  axi -symmetric  nozzle -receiver  was  outlined. 

-  33  - 


mrartyn. 


Cutaways  in  nozzle -guide  to 


Fig.  9.  Schematic  of  Apparatus  for  Studying  Steady-State  Characteristics  of  a 
Three-Way  Jet  Modulator. 
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Pressure-Flow  Characteristics  of  a  Three-Way  Jet 
Modulator 


The  experiment  was  designed  to  measure  the  dynamic  characteristics  of 
the  jet-receiver  interaction  region  (termed  the  "source”)  as  viewed  from 
a  point  downstream  in  the  receiver.  Provisions  were  made  for  introducing 
a  pressure  pulse  at  a  point  downstream  in  the  receiver,  measuring  the 
pressure  associated  with  the  incident  and  reflected  waves,  and  damping 
out  all  extraneous  waves.  The  experiment  thus  constituted  a  pulse  or 
step  response  of  the  interaction  region.  Having  obtained  a  dynamic 
characterization  of  the  source,  stability  of  systems  incorporating  this 
source  could  be  investigated. 

Since  the  last  writing,  measurements  have  been  made  which 
indicate  that  the  reflection  of  a  pulse  or  step  from  the  source  terminal 
can  be  characterized  simply  by  a  "reflection  coefficient"  for  frequencies 
up  to  at  least  several  hundred  cycles  per  second.  This  implies  that  for 
low  frequencies,  the  steady-state  pressure-flow  characteristics  are 
adequate  to  predict  the  interaction  region  phenomena.  Measurement 
techniques  and  typical  steady-state  characteristics  are  discussed  in 
Sec.  2.3.3. 

2.3.2.  "Creation"  of  an  Active  Source 

Before  a  study  of  the  dynamic  or  static  characteristics  of  the 
interaction  region  could  be  made,  it  was  necessary  to  establish  means 
for  "creating"  an  active  source,  i.e.,  a  source  which  "feeds"  an  im¬ 
pending  oscillation.  Certainly,  past  experience  with  jet  modulators  of 
various  types  is  sufficient  to  give  one  confidence  that  a  sustained  pressure 
oscillation  in  the  receiver  can  be  produced!  It  will  be  recalled  that  the 
axi~8ymmetric  nozzle -receiver -transmission  line -load  apparatus  was 
constructed  to  eliminate  all  unnecessary  solid  boundaries  in  the  inter¬ 
action  region  and  all  potential  receiver  cross -coupling  effects.  With 
the  apparatus  so  constructed,  the  influence  of  various  types  of  nearby 
bounding  surfaces  can  be'  studied  in  a  systematic  manner. 

Consequently,  the  apparatus  was  arranged  to  have  a  source 
driving  a  transmission  line  blocked  at  its  downstream  end,  since  a  sus¬ 
tained  pressure  oscillation  is  most  likely  to  occur  in  such  a  system.  Then, 
•various  solid  boundaries  were  placed  in  the  interaction  region  in  an  effort 
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to  create  a  "large"  amplitude  pressure  oscillation  in  the  line .  r 
configurations,  shown  in  Fig.  11,  were  found  to  produce  the  desired 
oscillation.  Case  (a)  typifies  the  phenomena  observed. 


By  proper  placement  of  the  wall  in  the  jet,  the  jet  is  caused 
to  attach  and  follow  the  wall  until  the  point  is  reached  where  an  adverse 
pressure  gradient  causes  it  to  separate.  Placement  of  the  receiver 
mouth  near  the  separation  point  results  in  a  sustained  pressure  oscilla¬ 
tion  in  the  transmission  line  and  a  corresponding  oscillation  of  the  jet. 

The  nature  and  cause  of  the  oscil’atior:  can  be  explained  qualitatively  as 
follows:  As  the  jet  is  "turned  on",  the  pressure  begins  to  rise  in  the 
receiver  as  a  result  of  simple  charging  of  the  associated  capacitance; 
the  pressure  in  the  separation  region  also  rises  above  its  "equilibrium" 
value  because  of  the  Id  gh  resistance  to  flow  introduced  by  the  wall-to- 
receiver  mouth  restriction;  the  pressure  in  the  separation  region  finally 
becomes  large  enough  (but  not  necessarily  above  atmospheric)  that  the 
jet  is  caused  to  separate  further  up  the  wall  and  thus  partially  miss  the 
receiver;  the  receiver  and  separation  region  pressures  then  tend  to  de¬ 
crease  until  the  original  separation  point  is  restored,  at  which  point  the 
process  is  repeated  again  and  again.  Of  course,  wave  motion  in  the 
receiver -line  plays  an  important  role  in  this  phenomena.  As  the  jet  is 
turned  on,  a  compression  wave  travels  down  the  line  and  is  reflected 
from  the  blocked  end  as  a  compression  wave.  Vtrhen  this  wave  returns 
to  the  receiver  mouth,  it  causes  the  separation  region  pressure  to  in¬ 
crease  rapidly  and  the  jet  to  suddenly  flip  away  from  the  wall.  The 
motion  of  the  jet  causes  a  sudden  decrease  in  the  pressure  at  the  receiver 
entry  and,  as  a  consequence,  causes  a  rarefaction  wave  to  travel  down 
the  line.  This  wave  is  reflected  from  the  blocked  end  in  like  sense,  and 
on  its  return  to  the  receiver  mouth,  reduces  the  separation  region  pres¬ 
sure  enough  so  that  the  jet  can  re -attach  at  the  original  point  of  attachment 
The  receiver  pressure  oscillation  is  very  nearly  triangular  in  shape. 

It  is  important  to  note  that  is  is  not  necess?.ry  to  have  a  wave- 
like  element  (i.e.,  a  transmission  line)  connected  to  the  source  in  order 
to  observe  a  sustained  oscillation.  If  a  tank  (simple  capacitance)  is  sub¬ 
stituted  for  the  transmission  line,  the  tank  pressure  will  oscillate  as  well. 
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In  this  case,  the  oscillation  is  characterized  by  an  exponential  rise 
followed  by  an  exponential  decay,  1  .e.,  a  direct  result  of  the  inherent 
RC  characteristics  of  the  receiver -tank  combination.  The  charging  and 
discharging  time  constants  usually  are  different,  the  discharging  time 
constant  genera1  !y  being  the  smaller  of  the  two. 

The  explanation  of  the  oscillation  phenomena  given  above 
applies  qualitatively  to  Cases  (b),  (c),  and  (d)  ot  Fig.  11  as  well.  The 
most  important  features  of  the  potentially  unstable  configuration  seem 
to  be:  (11  a  boundary  be  present  onto  which  the  jet  can  attach,  and,  (2) 
the  boundary  be  placed  slightly  into  the  jet  and  near  the  receiver  mouth 
so  as  to  create  a  semi-enclosed  separation  region. 

2.3.3.  Conceptual  Representation  of  a  Single -Receiver  Jet- 

Modulator  System  ' 

It  is  convenient  to  represent  the  system  being  studied  by  a  chain 
of  four  terminal  elements  as  she  m  in  the  block  diagram  of  Fig.  12,  where 
the  W's  are  mass  flow  rates  and  the  P's  are  static  pressures.  The  "source" 
includes  the  nozzle,  jet,  and  receiver  -diffuser .  The  "line",  if  one  exists, 
includes  any  reasonably  long  constant  diameter  passage*  which  connects  the 
source  to  the  load,  and  along  which  wave  motion  effects  are  important.  The 
"load"  includes  all  active  and  passive  elements  attached  to  the  line  or  directly 
to  the  source,  i.e.,  tanks,  orifices  (laminar  or  square  law),  jet-modulator 
control  ports,  etc. 

A  characterization  of  the  source  involves  establishing  the  functional 
re’ationship  W  =  f  (P  ,  geometry). 

E  3. 

2.3.4.  Measurement  of  Steady-State  Characteristics  of  an  Active 

Source  . ~ .  ~~  . .  . 

Reasonably  early  in  the  study,  physical  intuition  and  reasoning 
led  to  the  conclusion  that  the  steady-state  source  characteristic  for  a 
potentially  unstaole  geometrical  configuration  was  probably  S -shaped. 

Such  is  true  for  most  dynamic  systems  which  exhibit  a  tendency  toward 
self -excited  oscillations.  A  sketch  of  the  "anticipated"  active  source 
characteristic  is  shown  in  Fig.  13  along  with  a  typical  passive  characteristic. 


The  line  may  also  be  made  up  of  several  lengths  of  piping,  each  having 
a  different  diameter. 
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Fig.  12.  Block  Diagram  of  Single-Receiver  Jet-Modulator  System. 


Fig.  13.  Typical  Passive  and  Active  Source  Characteristics 


It  is  well  known  that  instability  is  most  often  associated  with  a 
blocked  load  condition.  Furthermore,  if  sufficient  fluid  is  allowed  to 
flow  from  the  load,  the  system  may  be  stabilized.  Considering  a  system 
consisting  of  a  source  plus  a  resistive  load,  a  blocked  load  (i.e.  ,  an 
infinite  load)  characteristic  superimposed  on  the  source  characteristic 
consists  of  the  zero-flow  axis,  as  shown  in  Fig.  13.  The  point  of  inter¬ 
section  is,  of  course,  the  simultaneous  "solution"  or  the  equilibrium 
condition  (if  one  exists).  To  better  underetand  the  relationship  which 
must  exist  between  load  and  source  resistances  in  order  to  produce  an 
instability,  it  is  helpful  lo  refer  to  Fig,  14.  Three  typical  cases  are  shown 
for  example  only.  This  in  no  way  presents  a  complete  picture  of  system 
stability,  as  will  become  evident  in  the  following  section.  In  each  of  the 
cases  shown  in  Fig.  14,  only  the  region  of  intersection  of  the  two  charac¬ 
teristics  is  shown.  Stability  "in -the -small"  may  be  determined  by  con¬ 
sidering  the  consequences  of  a  small  positive  pressure  disturbance  relative 
to  the  apparent  equilibrium  condition.  Below  each  W  versus  P  plot  is  a 
schematic  of  a  small  control  volume  supposed  to  be  situated  between  the 
source  and  load.  The  source  flow  enters  the  control  volume  and  the  load 
flow  leaves  the  control  volume;  equilibrium  requires  that  these  two  flows 
be  equal.  The  length  of  arrow  is  indicative  of  the  relative  magnitude  of 
the  corresponding  flows  entering  and  leaving  as  a  result  of  the  pressure 
disturbance. 

Both  cases  (a)  and  (b)  represent  stable  systems  since  a  positive 
pressure  disturbance  leads  to  a  load  flow  demand  which  is  greater  than  the 
source  can  provide,  thereby  causing  the  pressure  to  decrease  once  again. 

In  contrasi.  Case  (c)  represents  a  situation  where  the  positive  pressure 
disturbance  leads  to  a  load  flow  demand  which  is  smaller  than  t^e  source 
supplies.  As  a  result,  the  pressure  increases  even  further  until  eventually 
some  mechanism  places  a  limit  on  its  magnitude-  This  is  a  static  type  of 
instability.  It  will  be  seen  in  Sec.  2.  3.  5  that  three  types  of  instabilities 
can  be  identified,  one  of  which  is  a  static  instability. 

The  conventional  technique  for  measuring  the  portion  of  the 
characteristic  curve  in  the  first  quadrant  for  a  passive  source  involves 
varying  the  load  resistance  from  zero  (vertical  line)  to  infinity  (horizontal 
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line),  in  such  a  way  that  the  resistance  is  always  positive.  To  measure 
the  portion  of  the  curve  in  the  fourth  quadrant,  it  is  necessary  to  vary  the 
load  curve  slope  from  a  large  negative  value  to  a  small  negative  value. 

The  negative  slope  load  is  produced  by  supplying  flow  from  an  external 
source.  In  all  cases,  the  load  curve  passes  through  the  origin  of  the 
coordinates.  Figure  1  5  illustrates  schematically  the  experimental 
technique.  Notice  that  there  are  two  throttle  valves,  one  of  which  is 
generally  always  closed  (although  this  is  not  necessary). 

Obviously  the  technique  employed  to  measure  the  characteristic 
of  a  passive  source  would  lead  to  a  static  instability  in  the  cate  of  an 
active  source  (i.  e.  ,  see  Fig.  14c).  If,  however,  a  family  of  load  charac¬ 
teristic  curves  can  be  generated  which  ail  have  a  relatively  low  resistance 
(i.  e.  ,  large  change  in  flow  for  a  small  change  in  pressure)  in  the  vicinity 
of  their  point  of  intersection  with  the  source  curve,  stable  operation  can 
be  achieved  (i.e.,  see  Fig.  14  b).  A  possible  set  of  load  curves  is  shown 
superimposed  on  a  source  curve  in  Fig.  16.  It  is  not  necessary  that  all 
the  curves  have  the  same  slope,  only  that  the  slopes  be  very  near  zero  in 
the  region  of  potential  static  instability.  In  contrast  to  the  passive  source 
measurement  technique,  it  is  not  possible  for  all  the  lord  curves  to  pas3 
through  the  origin. 

An  experimental  technique  which  can  be  used  to  obtain  the  active 
source  characteristic  curve  is  shown  schematically  in  Fig.  17. 

The  critical  flowmeter  is  a  simple  orifice  in  which  the  flow  is 

always  choked.  Choked  flow  is  insured  if  the  supply  pressure  to  the  orifice 

is  maintained  always  at  least  twice  the  downstream  pressure.  Constant  flow 

is  achieved  by  maintaining  the  -  pressure  constant  at  a  reasonably 

large  value  compared  to  the  maximum  expecte-.  value  of  the  downstream 

pressure.  The  flow  rate  to  be  determined,  W  ,  is  equal  to  the  total  flew 

a 

as  indicated  by  the  downstream  rotameter  minus  the  constant  flow  from  the 
external  source.  In  the  first  quadrant,  W  is  positive  and  the  rotameter 
flow  is  greater  than  the  flow  from  the  constant  flow  source.  In  the  fourth 
quadrant  Wa  is  negative  and  the  rotameter  flow  is  less  than  the  flow  from 
the  constant  flow  source.  It  is  evident  that  if  the  flow  from  the  constant 
flow  source  is  adjusted  properly,  the  flow  rate  W  can  be  made  to  undergo 

cL 
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Experimental  Technique  for  Determining  the  Characteristic 
jf  a  Passive  Source. 


w 


w 


Fig.  16.  Concept  of  Point -By -Point  Measurement  of  an  Active  Source 
Characteristic . 
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large  changes  while  the  metering  pressure  Pm  undergoes  only  small 
changes.  The  progressive  variation  from  one  load  characteristic  curve 
to  another  is  accomplished  simply  by  varying  the  amount  of  downstream 
throttling.  Notice,  for  example,  that  the  load  curve  which  intersects  the 
source  curve  at  W  =  0  is  obtained  by  setting  the  downstream  throttle 

a 

such  that  P  =  P  =  P  ,  ,  .  In  this  case,  the  rotameter  flow  exactly 

a  m  ambient  *  7 

equals  the  flow  from  the  constant  flow  source. 

The  measurement  technique  outlined  above  has  been  used 
successfully  to  measure  steady-state  source  characteristics  for  one 
particular  active  source.  There  are,  of  course,  other  techniques  which 
could  be  used.  For  example,  one  promising  technique  involves  the  use 
of  a  very  long  line  connected  to  the  source  and/or  a  terminal  load 
attached  to  the  end  of  the  line  having  an  impedance  equal  to  the  line  surge 
impedance.  Any  pressure  waves  would  then  be  damped  due  to  friction 
and/or  absorbed  by  the  "perfectly  matched"  terminal.  Although  a  reason¬ 
ably  successful  attempt  was  made  to  construct  a  "perfectly  matched"  load 
terminal,  this  technique  was  abandoned  in  favor  of  the  other  technique 
because  of  superior  performance  and  simplicity  of  the  latter.  A  typical 
measured  characteristic  is  plotted  in  dimensionless  form  in  Fig.  18. 

The  source  geometry  was  chosen  such  that,  when  coupled  with  a  load, 
a  high  degree  of  instability  existed.  It  was  the  purpose  of  this  preliminary 
experiment  to  demonstrate  technique  only.  A  systematic  study  of  the  effect 
of  geometrical  variations  on  the  character  of  the  source  is  expected  to  be 
carried  out  in  the  future  by  another  investigator. 

2.3.5.  Analytical  Prediction  of  "Stability -in -the  Small" 

Consider,  as  before,  the  single -receiver  jet-modulator  system 
consisting  of  an  active  source  coupled  to  an  active  or  passive  load  terminal 
by  means  of  a  simple  transmission  line.  In  the  discussion  to  follow,  the 
concept  of  "load"  will  be  broadened  to  include  all  elements  coupled  to  the 
source.  For  the  particular  case  being  considered  here,  the  "load"  in¬ 
cludes  the  transmission  line  and  the  load  terminal.  Such  a  system  is 
interesting  as  well  as  practical  in  that  it  serves  as  a  potential  building  block 
for  more  complicated  systems. 
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Fig.  18.  Measured  Static  Characteristic  of  an  Active  Source. 


The  presence  of  the  "wave-like"  element  (transmission  line)  in 
the  system  suggests  that  consideration  be  given  to  system  dynamic  stability 
as  well  as  the  static  stability  discussed  in  the  previous  section.  Professor 
F.  T.  Brown  has  suggested  that  the  absolute  stability  "in-the-amall"  of  such 
a  system  be  classified  according  to  the  character  of  the  load.  The  following 
types  of  stability  (or  instability)  are  then  identifiable:  Static,  wave,  and 
surge.  These  three  types  of  stability  arise  on  considering  the  effect  of  the 
length  of  the  transmission  line  coupling  the  source  and  load.  Static  stability 
is  obviously  associated  with  a  sys'em  containing  an  extremely  short  line. 
Referring  again  to  the  discussion  of  static  stability  in  the  previous  section 
and,  in  particular,  to  the  sketches  in  Fig.  14,  it  is  possible  to  establish  a 
simple  criterion  for  static  stability  "in -the -small"  (i.  e,,  for  small  dis¬ 
turbances).  The  system  is  stable  in  the  vicinity  of  the  equilibrium  point 
(i.  e. ,  point  of  intersection  of  the  source  and  load  curves)  if  the  source 
impedance,  Z  ,  and  load  impedance,  Z#,  are  related  by  the  inequality 

—  +  —  >  o  (23) 

ZS  Zi 

Notice  that  a  positive  value  of  source  impedance  implies  a  negative  slope 
on  the  W  versus  P  plot,  whereas  a  positive  value  of  the  load  impedance 
implies  a  positive  slope  on  the  W  versus  P  plot.  This  is  merely  a  result 
of  sign  convention  for  the  flow,  i.  e. ,  positive  flow  is  defined  as  flow  from 
the  source  to  the  load. 

If  the  source -to -load  capacitance  is  also  very  small,  the  static 
instability  will  result  in  a  very  high  frequency  oscillation.  In  such  a  case, 
the  concept  of  static  stability  based  on  the  measured  steady-state  charac¬ 
teristic  curve  for  the  source  may  lose  its  meaning  because  of  the  inadequacy 
of  this  characteristic  at  high  frequencies.  For  this  case,  the  source  dynamic 
characteristics  must  be  considered  and  it  is  appropriate  to  speak  of  the 
corresponding  dynamic  stability  (as  distinct  from  wave  and  surge  types  of 
dynamic  stabilities  associated  with  transmission  lines). 

If  the  line  attached  to  the  source  is  very  long,  pressure  waves 
take  a  considerable  time  to  travel  down  the  line  and  then  return  again. 
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Over  the  period  of  the  round-trip  of  a  wave,  the  effective  load  impedance 

seen  by  the  source  is  the  "surge11  or  "characteristic"  impedance  of  the 

+ 

line  rather  than  the  impedance  of  the  load  terminal.  The  idealized  charac¬ 
teristic  impedance  of  a  line  of  cross-sectional  area  A,  in  terms  of  a  ratio 
of  pressure  to  mass  flow  rate  at  a  point  in  the  line,  is  given  simply  by 
Zc  =  o/A,  where  c  is  the  speed  of  sound  in  the  medium.  The  stability 
criterion  given  by  Eq.  (23)  applies  to  this  surge  situation  if  is  replaced 

simply  by  Z  .  Thus  the  system  is  surge -wise  stable  if  the  following  in- 

c 

equality  is  satisfied 


>  0 


(24) 


The  oscillation  frequency  associated  with  a  static  or  surge-type  instability 
is  characteristically  quit  high.  The  very  high  frequency  noise  associated 
with  most  jet -modulator  systems  is  undoubtedly  a  manifestation  of  the 
surge -type  instability. 

Since  the  surge  impedance  only  represents  the  true  load  im¬ 
pedance  seen  by  the  source  for  relatively  short  periods  of  time,  it  is 
necessary  to  consider  wave  motion  effects  if  stability  over  longer  periods 
of  time  <  of  interest  (and  it  most  often  is  of  prime  interest).  Wave  motion 
effects  obviously  have  their  greatest  importance  for  intermediate  length 
lines.  Here,  consideration  of  the  reflection  of  waves  from  both  ends  of  the 
transmission  line  leads  to  the  concept  of  "wave  stability." 


The  ratio  of  the  pressure  amplitude  of  a  wave  reflected  from  the 
termination  of  a  line  to  the  pressure  amplitude  oi  the  incident  wave  pro¬ 
ducing  the  reflection  is  known  as  the  "reflection  coefficient".  An  expression 
was  derived  from  first  principles  in  the  previous  report  (Ref.  1)  which  relates 
the  reflection  coefficient  to  the  terminal  and  line  characteristic  impedance  as 
follows : 


P  - 

■  =  r  = 

P  + 


Z  -  Z 
t  c 

z .  +  z 

t  c 


(25) 


neglecting  wall  shear 
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where  is  the  terminal  impedance.  A  close  examination  of  Eq.  (25) 
reveals  that  the  reflection  coefficient,  r,  may  be  positive  or  negative, 
must  be  equal  to  or  less  than  unity  if  the  terminal  is  passive,  and  is 
unbounded  if  the  termination  is  active. 

Consider  now  a  "lossless''  transmission  line  having  two  terminals, 
a  source  and  a  load.  For  the  sake  of  generality,  let  the  associated  ter¬ 
minal  impedances  be  denoted  as  and  Z^,  without  specification  as  to 
which  refers  to  the  source  and  which  to  the  load.  If  a  wave  reflects  first 
off  of  one  terminal  and  then  the  other,  that  is,  undergoes  one  complete 
cycle  returning  to  its  initial  position  and  direction  of  motion,  its  n*w 
amplitude  is  r^ r^  time-,  its  irdtial  amplitude.  The  product  r^r^  maV  he 
positive  or  negative  and  equal  to,  less  than,  or  greater  than  unity.  If 
I  rlrz\  >  the  wave  has  picked  up  energy  (from  an  active  source 

and/or  load)  and  the  system  is  unstable.  Such  a  result  was  clearly  evident 
in  the  limited  data  (oscilloscope  pictures)  presented  in  the  previous  report. 
One  of  two  types  of  instabilities  may  exist.  If  r  >  +  1,  the  instability 
is  "exponential"  in-the-small  (Type  I)  whereas  if  Tjr^<  -  1,  the  instability 
is  oscillatory  (Type  II).  This  follows  directly  from  consideration  of  more 
than  one  complete  cycle  of  travel  of  a  given  wave. 


The  criterion  for  wave -type  absolute  stability  is  therefore 


rlr2 


(26) 


A  condition  of  marginal  stability  exists  when 


Combining  Eqs.  (25)  and  (26)  gives  the  criterion  for  absolute  stability 
in  terms  of  the  surge  and  two  terminal  impedances.  Thus 


T*  T* 

/Zj  -Zc\  /z2-zc\ 

rl  2 

VZ>+  V  \Z2+Z‘/ 

1. 


(27) 


-  51 


Consider  again  the  two  types  of  instability.  For  Type  I,  Eq.  (27)  shows 
that  for  absolute  stability 

Zj  +  Z2  <  0,  if  Z2(or  Zj)  <  Zc  <  Zj  (or  Z£) 

Z,  +  Z0  •>  0,  if  Z  >  Z,  and  Z^, 

1  2  c  1  -  2 

and  for  marginal  stability,  Z^  +  Z2  =  0. 

For  Type  II,  Eq,  (27)  shows  that  for  absolute  stability 

Z,  Z0  +  Z  ^  <  0,  if  Z  >  Z,  and  Z-> 

12  c  c  l  -  2 

2 

Z.Z-,  +  Z  >  0,  if  Z  <  Z,  or  Z,  but  not  both. 

12c  cl  —  2 

and  for  marginal  stability, 

Z ,  Z.,  Z  =  0. 
i  2  c 

A  summary  of  the  stability  types  and  criteria  for  all  possible 
permutations  of  real  terminal  impedances  (resistances)  is  given  in 
Table  1.  This  summary  was  originally  prepared  by  F.  T.  Brown  and 
is  reproduced  from  Ref.  (19).  The  reader  interested  in  a  detailed 
treatment  of  stability  considerations  in  "two-line"  systems  should  refer 
tc  Ref.  19.  In  that  reference,  coupling  between  receivers  and/or  loads 
is  also  considered. 

2.3.6.  Graphical  Techniques  for  Prediction  of  System  Relative 
Stability 

The  stability  criteria  developed  in  the  previous  section  are 
particularly  useful  when  information  concerning  the  absolute  stability 
of  a  system  is  desired;  that  is.  is  the  system  stable  or  unstable  "in-the- 
small?"  Nothing  can  be  inferred  from  these  criteria  concerning  the 
relative  stability  of  the  system.  In  addition,  these  or  similar  criteria 
are  difficult  to  apply  in  cases  involving  non-real  terminal  impedances. 
Very  often,  a  knowledge  of  the  frequency  and  amplitude  of  a  pressure 
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Surge  Stability  Static  Stability 


oscillation  is  required.  For  example,  such  information  is  vital  in  the 
design  of  oscillators.  As  another  example,  consider  a  system  which  ex¬ 
hibits  an  instability.  The  over-ail  performance  of  the  system  may  not  be 
impaired  significantly  providing  the  oscillation  frequency  is  high  enough 
and  the  amplitude  is  low  enough;  again  a  knowledge  of  the  detailed  character 
of  the  oscillation  is  required  before  its  effect  can  be  fully  evaluated. 

The  suitability  of  graphical  techniques  such  as  the  method  of 
characteristics  (Ref.  20)  and  slope -line  integration  (Ref3.  21  and  22)  for 
dealing  with  transient  phenomena  in  systems  incorporating  simple  non- 
linearities,  wave -like  elements,  and  other  simple  active  and  passive 
elements,  is  well  ’.crown.  It  follows,  therefore,  that  these  graphical 
techniques  should  be  extremely  useful  in  studying  the  relative  stability 
of  simple  jet-modulator  systems  comprising  a  single  source,  transmission 
line,  and  load.  As  will  become  evident  in  the  discussion  which  follows, 
these  graphical  techniques  afford  an  amazingly  accurate  and  rapid  approach 
to  the  study  of  jet-modulator  system  stability.  Asa  demonstration  of  the 
utility  of  these  techniques,  three  interesting  examples  are  discussed  below. 

In  each  case,  the  analysis  is  semi -empirical  in  that  experimentally  derived 
(or  otherwise  known)  source  and  load  pressure-flow  characteristics  are 
utilized.  The  solutions  presented  are,  therefore,  limited  to  the  range  of 
frequencies  for  which  the  steady -state  source  and  load  characteristics 
are  valid.  In  all  cases,  the  basic  equations  describing  system  behavior 
are  presented  in  finite  difference  form,  since  the  graphical  techniques 
involve  a  type  of  numerical  step-by-step  solution. 

Example  1  -  Active  source  coupled  to  a  uniform  rigid  transmission 
line  blocked  at  its  downstream  end.  A  schematic  of  the  system  under  con¬ 
sideration  is  shown  in  Fig.  19.  The  solution  sought  is  as  follows:  given 
a  small  disturbance  in  the  load  end  pressure,  P^,  what  is  the  resulting 
variation  of  with  time,  i.  e.  ,  does  the  disturbance  die  out  or  does  the 
presence  of  the  active  source  lead  to  a  wave  instability.  Furthermore,  the 
detailed  nature  of  any  limit  cycle  is  to  be  determined. 

Before  proceeding  with  the  example  a  slight  digression  to  consider 
pertinent  details  of  the  transmission  line  portion  of  the  system  only  is  deemed 


-  54  - 


Active  source 


Active  Source  Coupled  to  a  Transmission  Line  Blocked 
at  its  Downstream  End. 


necessary.  The  line  may  be  considered  as  a  distributed  parameter 
system  but  only  small  amplitude  waves  are  admitted.  It  is  further 
assumed  that  the  gross  fluid  velocity  is  very  small  compared  to  the 
speed  of  sound  in  the  medium  and  that  wave  dispersion  effects  are 
unimportant.  The  assumptions  permit  use  of  the  well-known  linearized 
form  of  the  wave  equations.  Of  the  various  forms  of  the  solution,  the 
most  convenient  is  that  given  in  terms  of  the  conditions  at  the  ends  of 
the  line  as  follows: 


where 


P, 

cosh  P  Z  sinh  P 

b 

_ 

»  c 

W. 

1  ,  cosh  r 

b 

-  - sin  h  r i 

L-  J 

Z 

c  ' 

F=  propagation  operator 

^  characteristic  impedance  of  the  line;  defined  as  the 
complex  ratio  of  the  instantaneous  pressure  to  the 
instantaneous  flow  rate  at  a  given  point  in  the  line. 


(28) 


Both  r  and  Zc  are,  in  general,  frequency  dependent  properties.  Lu 
terms  of  the  physical  characteristics  of  the  line. 


and 


r(s)  =  L  ]Y(s)  .  Z(s) 


Zc(s) 


Z(s) 

Y(s) 


where 

L.  =  distance  separating  two  stations, 

Y  =  line  shunt  admittance  per  unit  length, 

Z  =  line  series  impedance  per  unit  length,  and 
s  =  LaPlace  operator. 

Brown  (Ref.  ?.3)  has  she  a  that,  for  the  case  of  uniform  attenuation  (i.e. , 
no  dispersion),  the  propagation  operator  can  be  approximated  as  follows: 

r  =  -  a  +  Ts  (29) 


-  56  - 


where 


T  =  characteristic  time  of  the  line, 

eQ  =  the  attenuation  factor;  a  is  termed  the  attenuation 
constant. 

For  the  special  case  of  an  ideal,  lossless  pneumatic  line  it  follows  that 


r  =  Ts 


where 


Z  = 
c 


=  Z 

CO 


cq  =  acoustic  speed  of  sound;  and 
A  =  cross-sectional  area  of  the  line. 


Experiment  and  analysis  shows,  however,  that  line  friction  is  important 
even  in  seemingly  short  lines.  As  will  become  evident  in  Sec.  2.  3.  9, 
reasonable  correlation  between  theory  and  experiment  requires  con¬ 
sideration  of  the  effects  of  line  friction  (and  corresponding  attenuation). 

Two  convenient  teciiniques  are  available  for  determination  of 
the  characteristic  parameters  T  and  Zc;  the  constant  R-L-C  model 
and  a  more  rigorous  and  exact  model  based  on  the  direct  solution  of 
the  continuity,  momentum,  and  energy  equations  which  describe  the 
process.  Although  considerably  more  complex,  the  "exact"  model  is 
generally  superior  since  it  includes  varying  velocity  profile  and  heat 
transfer  effects.  Nevertheless,  in  certain  cases  the  two  models  yield 
very  nearly  the  same  results  and  the  simpler  R-L-C  model  is  adequate. 
Brief  summaries  of  the  techniques  available'  for  computation  of  the 
characteristic  parameters  are  given  below. 

Constant  R-L-C  model  (Refs.  23,  and  26) 

The  series  impedance  and  shunt  admittance  are  given  by 


and 


Z  -  ID  +  R 

Y  =  CD 
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where 


I  =  inertance  per  unit  length  of  line, 

C  =  capacitance  per  unit  length  of  line, 

R  =  resistance  to  flow  in  line  per  unit  length 
D  =  d/dt 


Thus 


■  lV 


IC  D  +  RCD 


But  if  R  is  sufficiently  small  (as  it  usually  is)  or  if  R/lco  «  1 ,  where 
co  is  the  lowest  frequency  component  of  interest,  then 


T  = 


-  a  +  TD  « 


D 


Z  5I  (1  +  -5-  ) 

C  C  2  ID 


For  sinusoidal  excitation  the  following  limiting  forms  are  useful: 


1 


CO  I 

Z  a] 

"  A 

1  R 

jco  C 

r  -  LY 

RC 


CO 


Thus,  for  uniform  attenuation,  the  attenuation  constant  is 
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Q  = 


and 


for 


R 
to  I 


«  1 


_R_ 

u>I 


»  1 


If  it  is  assumed  that  the  flow  is  fully  developed  and  laminar  and 
that  the  line  friction  is  the  same  as  for  steady  flow,  the  attenuation 

-f 

constant  can  be  expressed  in  terms  of  line  and  fluid  properties  as 
follows: 
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The  assumption  that  the  resistance  is  the  same  as  for  steady,  fully 
developed  laminar  flow,  can  only  be  defended  for  cases  where  the  fre¬ 
quency  is  low  and  the  line  length-to-diameter  ratio  is  large.  At  the 
higher  frequencies  the  unsteadiness  of  the  flow  results  in  an  effectively 
higher  value  of  the  resistance.  Rohmann  and  Grogran  (Ref.  26)  suggest 
that  R  be  increased  by  20  per  cent  over  the  calculated  value.  Equally 
important  is  the  fact  that  in  many  practical  situations  (and  perhaps  the 
majority)  the  line  length-to-diameter  ratio  is  not  large  enough  to  discount 
the  effect  of  increased  resistance  to  developing  flow,  i.  e. ,  the  apparent 
friction  factor  (and  therefore  the  resistance)  is  significantly  greater  in  the 
’’entry  length"  of  a  tube  than  in  the  fully  developed  region  many  diameters 


y  =  kinematic  viscosity 
a  =  line  radius 
co  =  frequency 
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downstream.  Resistance  estimates  based  on  the  friction  factor  correla¬ 
tions  summarized  by  Shapiro,  et  al  (Ref.  27)  should  provide  corrections 
which  are  adequate  for  most  engineering  purposes.  These  correlations 
allow  computation  of  the  integrated  apparent  friction  factor  in  the  laminar 
entry  of  a  tube  as  a  function  of  diameter  Reynolds  number  and  the  tube 
L/D  ratio.  The  resistance  R  is  computed  easily  from  friction  factor 
information. 

"Exact'Model.  -  Derivations  for  the  propagation  operator  and 
characteristic  impedance  are  given  by  Brown  (Ref.  24).  The  model 
assumes:  cylindrical  tube,  one -dimensional  flow,  isothermal  walls, 
small  amplitude  laminar  disturbances,  and  continuous  medium.  Results 
pertinent  to  the  present  study  have  been  extended  in  range  and  replotted 
in  Figs.  20  and  21.  Attenuation  faxtors  as  a  function  of  line  and  fluid 
properties  may  be  obtained  from  Fig.  20.  Since  dispersive  effects  are 
neglected  in  the  present  example,  only  the  magnitude  of  the  characteristic 
impedance  need  be  considered.  Figure  21  may  be  used  to  obtain  a  correc¬ 
tion  to  the  nominal  value  of  the  characteristic  impedance  (i.e. ,  Zcq  =  cq/A). 

Having  established  means  of  computing  r  and  Z^  it  is  necessary 
now  to  rewrite  the  line  equations  in  a  form  suitable  for  step-by-step 
gtap’ deal  or  numerical  solution.  The  method  of  solution  requires  that 
the  equations  be  written  in  time  difference  form.  By  noting  that 


cosh  r  = 


r  -r 

sinh  r  =  — - — — 

2 


and 

TD 

e  F(t)  =  F(t  +  T) 

a  combination  of  Eqs.  (28)  and  (29)  yields  the  following  two  equations: 


Pb«  +  2T>  2a  .  Pb(t! 
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Attenuation  per  Cycle  of  Small  Amplitude  Sine  Waves  in  Rigid,  Isothermal  Wall, 
Laminar-Flow,  Circular  Air  Line. 


o 
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.Fig.  21.  Magnitude  of  Line  Characteristic  Impedance. 
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where 
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time  variable 

L/c  =  characteristic  time  of  the  line, 
acoustic  speed  of  sound  in  the  fluid  medium, 
characteristic  impedance  of  the  line, 
nozzle  mass  flow  rate,  and 
nozzle  supply  pressure. 


It  is  seen  that  events  occur  in  increments  of  the  characteristic  time  of 
the  line  T,  that  is,  the  time  required  for  a  wave  to  travel  from  one  end 
of  the  1  \e  to  the  other.  One  period  is  equal  to  2T. 

Returning  now  to  consideration  of  the  entire  system  of  Fig.  19, 
the  describing  equa^ons  in  addition  to  the  line  equations  are 
Source: 

W  /P 

-i-  =  f(~^ 

W  P 

s  Vs 

Load: 

W 

—  =  0;  Z  =  oo  (33) 

W 

s 

The  required  solution  i3  then  a  simultaneous  solution  of  Eqs.  (30)  through 
(33). 


,  geometry! 


(32) 
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Consider  first  the  special  case  of  a  lossless  line  for  which 
ea  =  1.  If  the  solution  is  assumed  to  "start"  at  an  arbitrary  non¬ 
equilibrium  condition  at  t  =  0,  a  combination  of  Eqs.  (30),  (31)  and 
(33)  yields : 


'6  M  =  Pa(T>  _  = 

>  P  /  P  P  P  w 

\  s  1 1  8  8  s  s 


for  the  period  t  = 


0  to  t  =  T,  whereas 


(34) 


P  (T)  Z  W  W  (T) 

3  ^  C  S  3 

p  p  w 

s  s  s 


for  the  period  t  =  T  to  t  =  2T. 

Analogous  results  for  the  more  general  case  (i.  e.  ,  ea  ^  1) 

are: 


(35) 


Pa(T> 


Pb(°) 


Z  W 

c  s  -a 
-  e 


W  (T) 
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n  Pa(T) 

+  (1  -  e  ~a)  ~ - 


(36) 
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Pb(2T) 


P  (T)  Z  W  W  (T)  P  (T) 

av  c  s  a  av  '  .  a  ,,  av 

-  =  -  e  - f-  (e  -  1)  — - 
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(37) 


Although  Eqs.  (34)  and  (35)  for  the  lossless  case  admit  simple 
graphical  interpretation,  Eqs,  (36)  ana  (37)  do  not.  Since  the  usual  case 
is  such  that  -  a  «  1,  Eqs.  (36)  and  (37)  may  be  approximated  by  the 
following  more  tractable  equations: 
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(39) 


A  simple  graphical  interpretation  of  the  various  system  equations 
is  presented  in  Fig.  22.  The  limiting  case  tor  which  ea  =  1  is  shown  by  the 
dotted  line  solution.  The  line  equations  may  be  rewritten  to  cover  the  time 
period  t  =  2T  to  t  =  4T.  The  graphical  solution  is  contJ.-u~d  in  the  same 
fashion  as  for  the  previous  period  of  time.  A  complete  s<  I  ation  for  the 
case  of  a  "lossless"  line  is  carried  out  using  an  actual  measured  source 
characteristic  in  Fig.  23a.  These  results  are  presented  in  terms  of  a 
phase -plane  plot  of  P,/P  versus  W  /W  .  The  corresponding  plot  of 
Pb/Pg  versus  t/T  is  shown  in  Fig.  23b.  It  is  evident  that  an  exponential 
wave  instability  exists  but  that  the  nonlinearity  associated  with  the  source 
characteristic  places  an  upper  bound  on  the  oscillation  amplitude.  Thus  a 
"limit  cycle"  or  "standing  wave"  exists  whose  frequency  is  the  quarterwave 
frequency  of  the  line  and  whose  amplitude  is  equal  to  3.  84  times  the  nozzle 
supply  pressure  for  the  case  of  a  supply  pressure  equal  to  10  inches  of  H^O. 
A  complete  solution  showing  the  effect  of  line  friction  in  a  3  foot  line  (0.  518 
inch  in  diameter)  is  shown  in  Fig.  24  for  comparison.  It  is  clearly  evident 
that  line  friction  is  significant  even  in  a  comparatively  short  line.  The 
resulting  limit  cycle  amplitude  is  reduced  to  Pb/Pg  =  1*94,  whereas  the 
oscillation  frequency  remains  the  quarter  wave  frequency  as  expected.  The 
accuracy^  of  the  graphical  method  of  solution  as  presented  here,  depends 
mainly  on  the  validity  of  the  assumption  that  -  a  «  1. 


The  term  accuracy  here  referr  to  the  method  of  solution  of  the  describing 
equations  and  not  to  the  quanti1  itive  comparison  of  the  results  with  actual 
experiment. 
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P'ig.  22.  Graphical  Interpretation  of  Eq.  (38). 


Fig.  23.  Graphical  Stability  Analysis  of  an  Active  Source  Coupled  to  a 
Lossless  Transmission  Lice. 


Graphical  Stability  Analysis  of  an  Active  Source  Transmission  Line 
System  Including  Effects  of  Line  Friction. 


For  the  sake  of  generality,  the  solutions  of  Figs.  23  and  24  were 
carried  out  on  a  nondimens ional  basis.  However,  although  the  dlmens ionle s s 
source  curve  (for  the  particular  geometry  tested)  is  relatively  independent 
of  nozzle  supply  pressure  over  the  range  investigated  (5  <  P  -  20  inch  K?0), 
the  corresponding  dimensionless  characteristic  impedance  of  the  line  is  not. 
The  dimensionless  characteristic  impedance  is  defined  as 


jL  - 
c 


Z  W 
c  s 


Thus  the  only  way  in  which~Z"  can  be  independent  of  P  is  if  W  /P  =  constant, 
that  is,  the  nozzle  is  a  laminar  flow  resistance.  In  most  practical  cases, 
the  nozzle  is  more  nearly  a  square-law  resistance.  For  the  nozzle  used  in 
the  experiments,  the  following  relation  can  be  established  from  the  nozzle 
calibration  data: 


14.6 

1  0.  515 


where  Pg  is  given  in  inches  of  H^O. 

Solutions  showing  the  dependence  of  the  predicted  limit  cycle 
amplitude  on  the  value  of  the  dimensionless  characteristic  impedance 
are  presented  in  Sec.  2.3.9. 

The  graphical  technique  demonstrated  here  is,  of  course,  not 
limited  to  an  infinite  load  terminal  impedance.  Solutions  may  be  carried 
out  for  any  real  finite  load  terminal  impedance,  whether  it  be  passive  or 
active.  It  is  interesting  that  a  relatively  small  bleed  orifice  placed  at  the 
end  of  the  line  can  completely  eliminate  the  limit  cycle  instability. 

Example  2  -  Active  source  coupled  to  a  load  chamber  (lumped 
compliance)  with  and  without  flow  out  of  the  chamber;  no  transmission  line. 
A  schematic  drawing  of  the  system  considered  is  shown  in  Fig.  25. 
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The  basic  nondimensional  equations  which  describe  the  system  behavior 
include 


Source: 


geometryj 


Chamber: 


W  W. 
a  b 


W  W 
s  s 


or  in  finite  difference  form 


6 


6t 


ave 


(40) 


(41) 


(42) 


Load  resistance  (assumed  to  be  linear): 


W  R  W  P 

S  (8  8 


where  r  is  the  dimensionless  time  variable.  These  three  equations 
contain  four  unknowns,  so  one  additional  equation  is  required.  Con¬ 
sideration  of  energy,  momentum,  and  continuity  for  the  sudden  expansion 
from  section  a  to  b  leads  to  the  conclusion  that  Pa/P3  =  f(P^/Ps,  Wa/W8, 

A  /A.)  where  A  =  cross-sectional  area  of  inlet  and  A,  =  cross-sectional 
a'  b  a  b 

area  of  chamber.  This  functional  relationship  follows  directly  from  the 
earlier  discussion  of  "sudden  expansion  diffusers."  If,  for  the  purposes 
of  this  example  only,  it  is  assumed  that  A^»  Aa,  then  it  follows  that 
p  /p  ~  P  /p  .  Therefore,  Eq.  (40)  for  the  source  characteristic  can 
be  written 
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W  /p, 

— =  { L± 

W  \  P 

s  \  s 

and  the  original  measured  characteristic  (W  versus  P  )  would  still  apply. 

EL  at 

Equations  (42),  (43),  and  (44)  may  be  interpreted  graphically  (for  the  first 
increment  in  dimensionless  time),  as  shown  in  Fig.  26. 

Figure  27  shows  complete  graphical  solutions  for  three  selected 
values  of  the  lead  resistance.  A  limit  cycle  is  seen  to  exist  for  the  cases 
where  the  load  resistance  is  high;  for  low  enough  values  of  load  resistance 
the  system  is  totally  stable.  The  limit  cycle  frequency  is  easily  determined 
from  the  plots. 

Since  the  method  of  solution  requires  the  use  of  an  approximate 
finite  difference  equation  to  represent  the  original  first  order  differential 
equation,  the  solution  is  expected  to  be  approximate.  Accuracy  is  generally 
assured  if  the  slope  of  the  "slope  line"  (1:6t/2)  is  maintained  sufficiently 
large.  In  other  words,  since  the  parameter  6t  contains  both  t  (the  time 
variable)  and  (the  chamber  volume)  it  suffices  to  always  maintain  At 
small  compared  to  the  integrating  time  constant  (i.  e. ,  pVPg/(3Wg).  If 
the  solution  is  carried  out  using  dimensionless  parameters  the  actual 
oscillation  frequency  or  period  for  a  given  volume  Vc  is  simple  determined 
by  appropriately  interpreting  the  oscillation  period  obtained  in  terms  of  r  . 

Consideration  of  the  limiting  case  for  which  Rj,  =  oo  and  V — 0 
yields  interesting  results.  As  the  volume  approaches  zero,  6r/2 — »0 
(i.  e. ,  the  "slope  line"  approaches  a  vertical  line).  It  is  then  obvious 
that  the  limiting  condition  results  in  an  infinite  frequency  limit  cycle 
bounded  by  the  "knees"  of  the  source  characteristic.  (See  the  dotted  lint 
in  Fig.  26).  This  result  is,  of  course,  impossible  and  therefore  of  academic 
interest.  Certainly  the  source  dynamics  must  become  significant  long  before 
the  theoretical  limiting  case  is  approached. 

Actually  the  mathematical  formulation  presented  above  is  over¬ 
simplified  for  some  practical  cases.  It  is  not  necessarily  valid  to  assume 
that  the  dynamic  behavior  of  the  system  is  governed  entirely  by  the  charging 
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Fig.  26.  Graphical  Interpretation  of  Eqs,  (40),  (42),  and  (43). 


Plus  Resistive  Load. 


and  discharging  of  the  lumped  capacitance.  Even  though  no  line  precedes 
the  pure  capacitance,  there  may  be  a  significant  inertial  effect  at  the 
entrance  due  to  the  acoustic  impedance  of  the  entrance  aperture.  In  fact, 
it  can  be  shown(Ref.  26)  that  the  "effective"  length  of  a  sharp-edged 
orifice  is  1.7a  where  -a  is  the  radius  of  the  orifice.  The  effect  of  the 
presence  of  a  line  preceding  the  lumped  capacitance  is  treated  in 
Example  3. 

Example  3  -  Active  source  coupled  to  a  transmission  line 
terminated  by  a  load  chamber  and  load  resistance.  Figure  28  shows 
a  schematic  of  the  system  considered. 


The  describing  time -difference  equations  are  given  below: 

Source: 


W  /  P 

a  _  f  I  a 

W  l  P 
s  y  $ 


geometry 


J 


Line:  (A  lumped  parameter  model  could  be  used  to  accurately  model  the 
line  for  frequencies  which  are  small  compared  to  the  fundamental 
natural  frequency  of  the  line  (See  Sec.  2.  3.  8) 
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W 

s 


2a  Wb<'> 
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(46) 


Wb(t  +  2T)  2a  Wb(t)  q  W  <t  +  T) 

- +  e  -  =  2e  - 

WWW 
s  s  s 

Pb(t+  2T)  2q  Pb(t) 
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(47) 


S5 
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Chamber  (See  Example  2): 


Load  Resistance: 


W  R  W  P 

S  [S3 


The  graphical  technique  involves  a  simultaneous  solution  of  the  describing 
eqiiations  as  in  the  previous  examples.  There  is,  however,  a  fundamental 
difference  between  this  example  and  the  previous  ones.  Only  one  "dynamic'' 
element  was  present  in  the  systems  of  Examples  1  and  2.  In  the  present 
example,  there  are  two  coupled  "dynamic"  elements  -  the  line  and  the 
chamber.  Consequently,  two  graphical  solutions  must  be  carried  on 
simultaneously,  with  the  resulting  "common  point  of  intersection"  being 
the  "answer".  A  graphical  interpretation  of  the  technique  is  given  in 
Fig.  29.  It  is  evident  that  the  minimum  value  allowed  for  6t/2  is  TW^/K,Po. 

A  typical  complete  graphical  solution  for  the  lossless  line  case 
is  shown  in  Fig.  30.  For  the  particular  conditions  chosen,  it  is  seen  that 
a  limit  cycle  instability  exists. 

The  limit  cycle  or  "surge"  frequency  can  be  determined  directly 
from  the  plot  by  counting  the  number  of  time  increments  required  to  complete 
one  encirclement  on  the  phase  plane  plot,  or  more  easily,  by  observing  the 
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Fig.  29.  Graphical  Interpretation  of  Eqs.  (45),  (46),  (47),  (48)  and  (49). 


Giaphical  Stability  Analysis  of  an  Active  Source  Coupled  to  a 
Lossless  Transmission  Line  Terminated  by  Capacitive  Plus 
Resistive  Load. 


period  on  the  time  plot.  The  rameters  which  determine  the  limit  cycle 
frequency  and  amplitude  for  .  given  source  characteiistic  are  obviously 
the  line  length  and  characteristic  impedance,  the  chamber  volume,  and 
the  terminal  impedance.  The  nature  of  the  solution  in  normalized  form 
reduces  the  '’free"  parameters  to  the  following  pair  of  ratios 


R, 

— — and 
Z 

c 


V 


V 

c 


Simple  analytical  means  for  computing  tne  "surge"  or  limit  cycle 
frequency  for  the  special  case  of  R,  =  oo  are  presented  in  Sec.  2.3.8. 
Unfortunately,  accuracy  of  solution  using  this  particular  techniovue  is 
not  assured  in  all  cases.  Since  the  system  contains  a  single  delay 
element,  (i.e.,  the  transmission  Une)  computations  yield  results  only 
at  the  points  in  actual  time  represented  by  t  =  2T,  4T,  6T,  .  .  .etc. 

As  in  Example  2,  good  accuracy  requires  that  the  slope  1:6t/2  be  large. 
However,  in  contrast  to  the  previous  exs  nple,  there  is  now  an  important 
restraint  placed  on  the  minimum  sise  of  6r/2  due  to  the  presence  of  the 
line.  It  easily  is  shown  that  unless  V-  /V  is  very  small  compared  to 
unity,  that  the  slope  1:  6r/2  may  be  tco  small  (or  1:Z  may  be  too  large) 
to  insure  adequate  accuracy  of  solution.  Obviously  simple  scaling  mani¬ 
pulations  on  the  abscissa  and  ordinate  do  not  help.  What  is  actually  re¬ 
quired  is  the  introduccion  of  an  artifice  that  will  effectively  reduce  the 
minimum  relative  value  of  6r/2.  Such  an  artifice  is  achieved  simply 
by  dividing  the  uniform  transmission  line  into  n-sections,  each. having 
a  length  L!  =  h/n  and  a  characteristic  time  T*  =  T/n.  The  chamber 
"slope  lines"  are  increased  in  number  and  slope  by  the  factor  n.  For 
example,  suppose  the  line  is  divided  into  4  equal  sections  as  shown  in 
Fig.  31.  Equations  identical  to  Eqs,  (46)  and  (47)  except  for  obvious 
changes  in  subscripts,  can  be  written  for  each  section.  The  chamber  slope 
lines  have  a  slope  1:5t/8.  It  is  possible  now  to  compute  the  values  of  all 
the  state  variables  of  interest  at  points  in  actual  time  represented  by 
t  =  2T',  4T1,  6T’...etc.  This,  of  course,  requires  that  the  state 

variables  at  the  line  points  1,  2,  and  3  also  be  computed  every  2T‘.  A 
graphical  interpretation  of  the  solution  teclmique  is  shown  in  Fig.  32 
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Fig.  32.  Graphical  Interpolation  of  Example  3  for  a  Line  Divided  into  Four  Sections 


Graphical  Interpr  .ation  of  Example  3  for  a  Line  Divided  into  Four  Sections 


for  the  case  of  n  =  4.  It  is  interesting  that  the  technique  yields  the  values 
of  the  state  variables  at  points  1,  2,  and  3  rather  "Automatically";  in  the 
numerical  solution  to  be  discussed  in  Sec.  2.3.7  these  values  must  each 
be  computed  individually  as  time  progresses. 

2.3.7.  Computer-Aided  Numerical  Prediction  of  System  Relative 
Stability -  - - - 

The  elegance  and  general  utility  of  the  graphical  techniques  pre¬ 
sented  in  the  previous  section  are  obvious.  These  techniques  provide  in¬ 
sight  into  problem  formulation  and  solution  which  is  exceedingly  difficult 
to  realize  with  other  numerical  methods.  Nevertheless,  for  detailed  design 
studies  involving  extensive  computations,  only  machine  computation 
affords  the  necessary  speed  and  over-all  flexibility.  Techniques  for  solving 
the  previous  three  example  problems  numerically  with  the  aid  of  a  digital 
computer,  are  demonstrated  below.  Results  of  actual  computations  are  given 
in  Sec.  2.  3.9.  Fortran  program's  for  each  cf  the  Examples  are  available 
from  the  writer  and  are  to  be  presented  in  the  writer's  ScD.  thesis 
(Summer  1964). 

Example  1  -  Active  source  coupled  to  a  transmission  line  blocked 
at  its  downstream  end.  Equations  (30)  through  (33)  fully  describe  the 
behavior  of  this  simplified  system.  The  source  characteristic  written 
in  functional  form  in  Eq.  (32),  is  experimental  in  nature  and  therefore 
may  be  presented  in  graphical  or  tabular  form.  Tabular  representation 
is  preferred  (and  assumed)  here. 

The  results  of  the  previous  graphical  solutions  indicate  that  events 
occur  in  increments  of  the  characteristic  time  of  the  line,  T.  Accordingly, 
the  simultaneous  solution  of  the  governing  equations*  is  most  easily  realized 
by  "computing"  changes  every  2T  units  of  time.  If  it  is  assumed  that  the 
state  of  the  system  is  known  at  time  t  =  0,  then  subtraction  and  addition 
of  Eqs.  (30)  and  (31)  yield  the  following  two  relations  (W^/Wg  =  0). 

(M 

CO 

oo  *  The  "exact"  form  of  the  line  equations  may  be  used  here  with  no  in- 

£  crease  in  difficulty. 
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It  is  seen  that  the  pressure  at  the  downstream  end  of  the  line  at  time  t 
and  t  +  2T  is  related  to  the  pressure  and  flow  at  the  upstream  (source) 
end  at  time  t  +  T.  The  source  characteristic  is,  of  course,  good  for 
all  time.  Thus  from  Eq.  (32)  we  have 

(32a) 


Wa  (t  +  T) 
W 


=  f 


P  (t  +  T) 
a' 


Starting  with  a  given  set  of  initial  conditions  (i.e.  ,  at  t  =  0),  the  solution 
involves:  1)  simultaneous  solution  of  Eqs.  (32a)  and  (50)  to  obtain  P  (T)/P 

3.  S 

and  W  (T)/W  ,  and  then  Z)  solution  of  Eq.  (51)  to  obtain  P,  (2T)/P  .  The 
solution  for  each  succeeding  increment  in  time  2T  follows  in  exactly  the  same 
way.  Due  to  the  nature  of  the  source  characteristic,  it  is  obvious  that  the 
simultaneous  solution  of  Eqs.  (32a)  and  (50)  involves  iteration  and  inter¬ 
polation  (unless,  of  course,  an  algebraic  relation  is  established  for  the 
source  characteristic).  A  useful  and  completely  general  approach  is  as 
follows : 


1.  From  Eq.  (50)  define  a  new  function,  WP,  such  that 


WP  =  Z 
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Thus,  an  exact  simultaneous  solution  of  F.qs.  (32a)  and 
(50)  would  give  WP  =  0. 

2.  Starting  with  the  first  entry  in  a  table  of  values 
representing  the  source  characteristic,  compute 
WP  for  successive  points,  noting  that  P.  (t)/Pg  is 
known  and  W  (t  +  T)/W  and  P  (t  +  T)/P  are  the 
values  from  the  table.  The  two  points  in  the  table 
which  bracket  the  exact  solution  is  found  by  noting 
the  change  in  sign  of  WP. 

3.  Employ  any  standard  random  or  systematic  iteration 
technique  to  find  the  unknown  point  within  an  acceptable 
tolerance.  One  useful  and  efficient  iteration  technique7^ 

nvolves  successively  halving  the  interval,  computing 
W  /W  ,  interpolating  to  find  P  /P  computing  WP,  and 
testing  WP  against  an  acceptable  tolerance  B.  The 
process  is  repeated  (usually  6  or  7  times  for  a  reasonable 
tolerance)  until  jwp{  -  B.  For  the  step  involving  inter¬ 
polation,  the  well-known  parabolic  interpolation  technique 
(Ref.  29)  should  be  adequate  for  most  purposes.  The 
final  results  of  the  above  iteration  can  then  be  utilized 
in  Eq.  (51)  to  obtain  P^(2T).  Then  the  encire  process 
can  be  repeated  in  steps  of  2T  until  the  complete  solution 
is  obtained. 

Example  2  -  Active  source  coupled  to  a  load  chamber.  The 

approach  is  fundamentally  the  same  as  in  Example  1  ;  however,  the 
governing  equations  are  different  and  a  slight  modification  is  necessary 

in  the  source  characteristic  searching  procedure.  Equation  (42)  des¬ 
cribing  the  charging  process  for  the  chamber  can  be  rewritten,  and  sub¬ 
divided  into  two  parts  as  follows: 


The  table  is  conveniently  arranged  so  that  the  column  of  W  /W 
in  equal  increments.  '  3 


changes 
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The  source  and  load  characteristics  become 
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— -  a  f  _ 


n  [~  P.(k+1) 

-  =  f  — - - 

P 

j  L  3  j 


Wb(k) 


R  W 
(  s 


Ph(K) 
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Wfe(k+  2) 


R  W 
l  s 


Pb(k+2) 


(55b) 


Starting  with  a  set  of  initial  conditions  (k  =  0),  the  solution  involves: 

1)  simultaneous  solution  of  Eqs.  (52),  (54),  and  (55a)  to  obtain  P,  (1)/P 
and  W  (1)/W  ,  and  2)  solution  of  Eqs.  (53),  and  (55b)  to  obtain  P,  (2)/P  . 

The  solution  for  each  succeeding  value  of  k  follows  in  exactly  the  same 
way.  An  examination  of  the  graphical  solution  of  this  particular  problem 
indicates  that  a  range  exists  for  which  the  "slope  line"  crosses  the  source 
curve  in  three  places;  thus  three  solutions  are  possible  but  only  one  is  the 
desired  solution.  The  same  logic  which  leads  to  a  selection  of  the  correct 
value  when  employing  the  graphical  technique,  must  obviously  be  incorporated 
into  the  numerical  solution.  The  general  approach  is  as  follows: 
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1.  From  Eq.  (52)  define  a  function,  WP,  such  that 


6  T 

wa(k+ 1) 

+ 

1  J  £i )  P8 

pb(k) 

Pb(k+1) 
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P 
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2.  Starting  with  the  first  entry  in  the  t  ible  of  values 
representing  the  source  characteristic,  compute  WP 
for  successive  points  until  the  first  "knee''  of  the  curve 
is  exceeded.  Care  should  be  taken  to  stop  before  the 
second  "knee"  is  reached.  A  change  in  sign  of  WP 
indicates  the  bracket  about  the  desired  solution..  If 

a  change  in  sign  is  not  found,  start  with  the  last  entry 
in  the  table  and  proceed  backwards  through  the  table 
computing  WP  for  each  point. 

3.  Employ  the  earlier  mentioned  iteration  and  inter¬ 
polation  routine  to  find  the  desired  point  within 
acceptable  tolerance. 

4.  For  the  next  complete  computation  reverse  the 
procedure  in  step  (2)  above. 

In  Step  (3)  of  the  above  procedure,  logic  decisions  must  be  employed 
which  indicate  at  which  end  of  the  table  the  computations  are  being  made. 
Furthermore,  certain  of  the  basic  formula  must  be  appropriately  changed 
to  account  for  the  direction  of  travel  and  placement  in  the  table. 

Example  3  -  Active  source,  coupled  to  a  transmission  line  ter¬ 
minated  by  a  load  chamber  and  load  resistance.  Here  the  most  interesting 
and  useful  solution  is  the  one  which  allows  the  division  of  the  line  into 
n-sections.  The  solution  proceeds  in  much  the  same  manner  as  for  the 
previous  Examples  except  that  it  is  more  lengthy  and  does  not  involve 
any  possible  triple  value  results  as  in  Example  2.  Again  it  is  necessary 
to  start  with  a  set  of  initial  conditions  for  all  state  variables  and  then 
proceed  with  a  simultaneous  solution  of  the  governing  equations.  The  pro¬ 
cedure  involves  the  following  successive  steps: 


1.  Define  a  function,  Wt-',  such  that 


WP  = 


Z  W 
c  s 


W  (t  +  T) 
a  ' 


wb(t)  Pb(t) 

W  P 

3  -J  S 


P  (t  +  T) 
a 

P 

s 


2.  Starting  with  the  first  entry  in  the  table,  compute 
WP  for  each  successive  entry  and  find  the  points 
in  the  table  which  bracket  the  desired  point,  i.e,, 
W  (nT'}  and  P  {nT”).  Note  that  T'  =  T/n. 

cl  cl 


3. 

4. 

5. 


Knowing  the  initial  values  of  the  n-1  intermediate 
state  variables,  compute  in  succession  P^j^fn  +  lTj T 1  , 

P2f (»T')J  ,  P3|jn-I)T«]  ....  _ j[J( 3T ’ )H . 

Compute  Pb(2T')  and  Wb(2T'). 

Return  to  step  (2)  to  find  W^j^jn  +  2)T^j  and 

P  r~(n  r  2)  T '  I  and  then  continue  on  in  the  same  manner 
a-L  _j 

to  find  successive  values  of  P,  and  W,  . 

b  b 


2 .  3.8.  Analytical  Prediction  of  Surge  Frequency  for  a  Transmission 
Line  I’erminated  by  a.  Chamber. 

Of  the  three  examples  presented  in  the  previous  sections,  Example  3 
represents  the  general  case.  Examples  1  and  2  are  merely  limiting  cases 
of  Example  3.  It  is  possible,  without  carrying  out  any  graphical  or  numerical 
solutions,  to  predict  the  surge  frequency  of  the  system  in  Example  j  for  two 
important  special  cases:  1)  =  oo  emd  2}  R^  =  0.  The  surge  frequency  is 

the  lowest  natural  frequency  of  the  system  and  is  the  frequency  at  which  the 
limit  cycle  is  most  likely  to  exist.  Consider  the  first  case  of  infinite  load 
terminal  impedance.  Two  approaches.  are  available,  depending  on  the  nature 
of  the  model  for  the  line.  Either  the  lumped  or  distributed  parameter  model 
may  be  used. 


Distributed  Parameter  Solution.  Assume  that  the  line  is  lossless 
and  that  the  effective  bulk  modulus  of  the  fluid  in  the  line  is  equal  to  that 
in  the  chamber.  The  governing  equations  for  the  lossless  line  are: 
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P,  =  cosh  TD  P  -Z  sinh  TD 

b  L  J  a  c  L  J 


1  I  f -  - 

W,  a  -  sinh  TD  P  +  cosh  TD  W 

b  Z  I  a  L 

cl— 


Considering  the  terminating  chamber  as  a  pure  capacitance  gives 


wb  = 


Combining  Eqs.  (56),  (57),  and  (58)  and  eliminating  W  and  W,  gives 


Z  sinh  TD 
c 


(3  j  Zc  tan  u)  t 


where 


<jj  =  frequency 


The  surge  frequency  can  be  determined  by  setting  the  characteristic 
equation  of  the  above  transfer  function  equal  to  zero.  Thus 


Since 


(3T  _ 

V  Z 
c  c 


uT  tan  'j  T 


T  =  JZ  and  Z 

V  fl  c  V  A 


we  find  that 


-  r|  tan  r) 
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where 


Vf  =  volume  of  line  =  Ai 
rj  =  dimensionless  surge  frequency  -  uT 

By  noting  that  the  fundamental  natural  frequency  of  the  line  with  zero 
terminating  volume  is 


n 


l 


speed  of  sound 
4< 


1 

4T 


and  setting  u>  =  2ttQ  „  ,  we  obtain 
3  cf 


Q 

c( 


Tf 


(62) 


Equations  (61)  and  (62)  relate  the  dimensionless  surge  frequency  to  the 
line-to-chamber  volume  ratio.  The  resuPinr  relation  for  the  fundamental 
is  plotted  in  Fig.  33.  For  the  ca3e  when  =  0,  it  is  necessary  to  consider 
the  multi-valued  character  of  Eq.  (61) 

Lumped  Parameter  Solution.  When  the  surge  frequency  is  small 
compared  to  the  fundamental  natural  frequency  of  the  line  alone,  a  simple 
lumped  parameter  model  of  the  line-tank  system  is  adequate.  For  example, 
consider  the  line  as  a  ’’pure"  inertance  (no  capacitance). 


P 

a 


A 


d  W, 
b 

dt 


(63) 


'ombining  Eqs.  (58)  and  (63),  and  solving  the  resulting  characteristic 
equation  for  the  surge  frequency  gives 


from  which  follows 
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(64) 


tw.1 


Lumpe  d 


quency  of  a  Lossless  Transmission  Line  Terminate  .  by  a  Simple  Tank 


This  expression  is  compared  to  the  distributed  parameter  solution 
in  Fig.  33. 

2.  3.  9.  Comparison  of  Analytical  and  Empitica.!  Predictions 
With  Fxperimental^ieasurements 

In  order  to  demonstrate  the  validity  of  the  graphical  and  numerical 
solution  techniques  presented  in  Secs.  2.3.6  and  2.3.7,  results  of  a  re¬ 
presentative  set  of  digital  computer  solutions  and  experimental  measurements 
are  presented  below. 

Example  1  -  Active  source  coupled  to  a  transmission  line  blocked 
at  ns  downstream  end.  Two  types  ot  experimental  measurements  were 
made:  1)  pressure  oscillation  amplitude  at  the  downstream  end  (point  of 
maximum  pressure  for  the  standing  wave),  and  2)  oscillation  frequency  and 
waveform  characteristics.  Figure  34  shows  a  set  of  oscilloscope  traces 
which  indicate  the  character  of  the  waveform.  If  there  were  no  wave  dis¬ 
persion  and  line  friction  effects,  the  waveform  would  be  a  square  wave 
since  all  odd  numbered  harmonics  must  be  present;  this  follows  directly 
from  the  multivalued  character  of  the  line  equations.  Higher  harmonics 
ai  j  seen  to  be  present  in  the  pictures  of  Fig.  34,  but  they  are  highly 
attenuated  because  of  friction  and  dispersion.  Oscillation  frequency 
estimated  from  the  apparent  fundamental  in  the  pictures  compares  very 
well  with  the  expected  q\  arter-wave  freouency. 

Oscillation,  amplitude  measurements  were  particularly  difficult 
for  this  system  because  of  the  sporatic  presence  of  higher  harmonics  and 
a  very  low  frequency  modulation.  It  is  not  immediately  clear,  in  fact,  what 
amplitude  should  be  measured.  A  single  frequency  spectrum  analysis  was 
carried  out  which  indicated  the  presence  of  3rd  and  5th  harmonics  with 
approximately  the  expected  amplitudes.  This  suggested  that  the  most 
meaningful  measurement  would  be  the  amplitude  of  the  fundamental.  A 
square  wave  having  an  amplitude  tt/4  times  this  fundamental  should  compare- 
with  the  predicted  "square -wave"  solution  within  the  validity  of  the  original 
model.  However,  instrumentation  was  not  available  which  could  be  used  for 
making  this  measurement  accurately.  If  no  dispersive  eifects  were  present, 
a  measurement  of  the  true  rms  amplitude  of  the  waveform  would  yield  the 
amplitude  of  the  desired  square  wave.  Since  a  true  rms  voltmeter  was 
available,  measurements  were  carried  out  in  this  manner.  A  few  spot 
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b)  L  *  1  FT  el  L  =  4  FT 


c'l  L  =  ! 5  n 


Legend- 

V  =  vertical  scale  in  inches  H^O/division  (approx.) 

H  =  horizontal  scale  in  milliseconds/division 
L  =  Line  length  in  feet 

P  =  nozzle  supply  pressure  =  10  inches  H-,0  for  pictures  (a) 
through  (e). 

D  -  line  diameter  =  0.  518  inch 


Fig.  34.  Pressure  Oscillations  at  the  Blocked  End  ( Downs  tr earn)  of  a 
Constant  Diameter  Transmission  Line  Coupled  to  an  Active 
Source.  (Conf.  3). 
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times  the  fundamental  amplitude  and  the  true  rms  amplitude  was  easily 
within  5  per  cent.  Figure  35  shows  the  results  of  the  rms  amplitude 
measurements.  The  results  of  digital  computer  predictions  which  in¬ 
cluded  the  fraction  correction  from  Fig.  20,  are  shown  for  comparison  in 
Fig.  35b.  Although  the  comparison  appears  generally  poor,  several 
important  results  are  apparent.  The  model  is  adequate  to  predict  the 
qualitative  behavior  for  the  longer  length  lines;  however,  the  friction 
correction  seems  too  severe  to  provide  good  quantitative  comparison. 

The  solution  ls  highly  sensitive  to  the  friction  factor;  for  example,  eQ 
for  the  one  foot,  line  is  0.  977  which  results  in  a  factor  of  two  reduction 
of  the  amplitude  as  compared  with  the  results  of  the  lossless  solution. 


Frequency  limitation  in  the  static  source  characteristic  is, 
perhaps,  a  contributing  factor  to  the  poor  correlation  at  the  higher 
frequencies  (i.e.,  shorter  lines).  Several  cf  the  basic  assumptions 
underlying  the  line  equations  are  open  to  some  question  for  the  shorter 
lines.  The  inclusion  of  dispersive  effects  might  improve  the  correlation. 
Nevertheless,  the  correlation  is  not  actually  very  bad  considering  the 
particular  difficulties  of  experimental  measurement  and  the  extreme 
sensitivity  to  the  value  of  the  attenuation  factor.  Several  improvements 
in  the  model  which  would  lead  to  better  correlation  with  experiment,  are 
presently  under  study. 

Example  2  -  Active  source  coupled  to  a  simple  chamber.  All 
measurements  wei  e  made  with  line  lengths  greater  than  6  inches  preceding 
the  chamber.  Pertinent  experimental  results  are  presented  in  connection 
with  Example  3.  The  apparent  transition  between  the  systems  of  Examples  2 
and  3  w'ill  be  explored  in  more  detail  in  the  writer's  ScD.  thesis. 

Example  3  -  Active  source  coupled  to  a  transmission  line  terminated 
by  a  chamber  and  load  resistance.  It  was  established  earlier  that  the  pertinent 
free  parameters  are  /V  and  R^/Z^.  Only  the  infinite  resistance  case 
(R^.  =  oo)  is  considered  here.  Figures  36  and  37  show  typical  waveforms 
for  various  values  of  V., /V  .  Unlike  the  results  for  Example  1  the  waveforms 
are  very  nearly  sinusoidal.  The  chamber  coupled  to  the  line  introduces  a 
type  of  ’’filtering"  effect,  such  that  the  higher  harmonics  do  not  generally 
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Pressure  Oscillation  Amplitude,  Pu/P„  Pressure  oscillation  amplitude , 


Vvc 


VVC 

0.50 


0.40 


All  Pictures: 

Vertical  scale  =  10.2  inches  H^O/division 

Horizontal  scale  =10  milliseconds /division 

Nozzle  supply  pressure  =  10  inches  H2O 

Line  length  =  12  inches 

Line  diameter  =  0.  518  inch 

V./V  =  Line -to -chamber  volume  ration 
f. 1  c 

Fig.  36.  Pressure  Oscillations  in  a  Blocked  Chamber  Coupled  to  an 
Active  Source  (Conf.  3). 
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All  Pictures : 

Vertical  scale  =  10.  Z  inches  H^C/di vision 

Horizontal  scale  =  10  milliseconds/division 

Nozzle  supply  pressure  =  10  inches 

Line  length  =  6  inches 

Line  diameter  =  0.  518  inch 

V./V  =  Line -to -chamber  volume  ratio 
l  c 

Fig.  37.  Pressure  Oscillations  in  a  Blocked  Chamber  Coupled  to  an 
Active  Source  (Conf.  3). 
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appear.  As  /V  becomes  very  srriall  the  situation  of  Example  1  is 
approached  and  higher  harmonics  begin  to  appear.  For  longer  lines  the 
waveforms  become  less  sinusoidal  and  higher  harmonics  begin  to  be  more 
apparent.  A  comparison  of  measured  and  theoretical  limit  cycle  frequencies 
is  given  in  Fig*  38.  The  comparison  is  remarkably  good  even  for  such  short 
lines,  thereby  showing  the  validity  of  the  mathematical  formulation  for  the 
range  of  conditions  considered. 

Oscillation  amplitude  measurements  were  carried  out  using  the 
true  rms  voltmeter.  In  this  case  the  rms  value  and  the  amplitude  of  the 
fundamental  are  assumed  to  be  the  same.  Typical  measured  results  are 
shown  in  Fig.  39.  As  expected,  the  amplitude  decreases  with  increases 
in  both  chamber  size  and  supply  pressure.  The  decrease  in  amplitude 
with  supply  pressure  is  a  result  of  an  increase  in  the  stiffness  of  the  jet 
with  supply  pressure  increase. 

A  typical  comparison  of  predicted  and  measured  results  is  shown 
in  Fig.  40.  The  correlation  is  remarkably  good  considering  the  difficulties 
of  measurement  and  the  sensitivity  of  the  solution  to  the  friction  correction. 
Since  the  frequencies  were  low  compared  to  the  natural  frequency  of  the  line 
alone,  the  attenuation  factor  was  computed  from  the  lumped  parameter 
approximation  (see  Sec.  2.3.6).  A  minor  improvement  (in  the  favorable 
direction)  could  be  achieved  by  using  the  more  exact  values,  but  this 
requires  lengthy  computation.  Typical  values  of  the  attenuation  factor 
are:  ea  =  0.  9988  for  L.  -  1  ft.  and  ea  =  0.  9953  for  L  =  4  ft.  Sensitivity 
to  changes  in  the  fourth  digit  are  apparent! 

Cone lusions  -  It  is  belie  »<~d  that  the  experimental  measurements 
conducted  to  date  adequately  demonstrate  the  importance,  general  usefulness, 
and  limitations  of  the  prediction  models.  Certain  refinements  are  needed  in 
the  models,  but,  for  the  majority  of  practical  cases,  the  predictions  are  in 
remarkably  good  agreement  with  the  experiments. 
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Line  Terminated  by  a  Simple  Chamber. 
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Nozzle  supply  pressure.  s 

.  Exne:  iroental  Measurements 
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3.  DYNAMIC  ANALYSIS  OF  CONTROLLED  JET  DEFLECTION 


3.1.  LaPlace  Transform  Solution  of  the  Jet  Equations 

In  an  earlier  report  ,  Ref.  1,  the  equations  describing  the  motion 
of  the  jet  were  written  directly  in  terms  of  the  values  of  variables  at  dis¬ 
crete  points.  These  equations  in  actuality  assumed  certain  types  of 
approximations  for  the  space  derivatives  of  the  transverse  jet  velocity 
and  position.  There  was  no  standard  solution  against  which  co-  1  ■'e 
compared  vhe  digitally  computed  solution  and  consequently  no  way  of 
knowing  which  derivative  approximations  w'ould  give  the  best  results. 

The  purpose  of  this  section  is  to  compare  the  solution  of  the  jet  equations 
using  various  space  derivative  approximations  against  two  exact  solutions. 

3.1.1.  Rederivation  of  Jet  Equations 

Consider  a  segment  of  the  jet  with  stagnant  fluid  around  it  as 
in  Fig.  41 . 

The  momentum  equation  gives: 

Phdf  -  (P  +  -  —  -  dn)  hd/  =  — - —  (phd  £  dnV) 

3  n  3  t 


+  (V  +  d t)  V.phdn  -  VV.phdn  ('  ">) 

3 1  J  J 


ap  3V  .  -r  OV 

-  =  P  -2—  x  P  V.  - 

3  n  3 1  ^  3 1 


The  appropriate  nondimensionalized  variables  are: 


V. 

3 


(66) 
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h  =  dVpth 
p  =  density 

Vj  =  axial  jet  velocity 
V  =  transverse  jet  velocity 
P  =  pressure 

Fig.  41.  Control  Volume  Around  a  Segment  of  Jet. 
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The  quantity  l/t  can  be  defined  as  the  Strouhal  number  for  the 
modulator.  Equation  (66)  becomes 
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The  jet  position  equation  is(Ref.  1): 


+  Y 


We  can  combine  Eqs.  (68)  and  (69)  and  obtain  (dropping  the  (  )  notation): 


i*  ♦ 

at2  at9i 


This  is  a  second-order  partial  differential  equation  of  the  parabolic  type. 
The  forcing  term  of  Ea.  (70)  can  be  evaluated  as  follows: 
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=  F 


(71) 


where 

P  -  bubble  pressure 
P&  =  atmospheric  pressure 

Thus  Eq.  (70)  becomes 


2  2  2 

O  +  2  afx  + 

at2  ata;  a;2 


(72) 


3.  1. 2.  Step  Change  in  Pressure  Across  Jet 
Consider  a  step  change  in  F  which  acts  over  a  finite  length 
of  jet  as  shown  in  Fig.  42. 


The  expression  for  F  is: 


F  =  FoUs(t)  (Vg(  f)  -  U3(f  -  lQ)  (73) 

where  Ug  is  the  unit  step  function.  The  LaPlace  Transform  with  respect 
to  (  of  Eq.  (72)  with  Eq.  (73)  inserted,  is: 


+  2q 


dY 

dt 


+  q2  Y  =  F  U  (t)  I- 
o  s  i 

\q 


~*oq  \ 

-  =  F  U  (t) 

/  os 

q  ) 


(74) 


where 


(75) 


Equation  (74)  is  now  an  ordinary  differential  equation  which  can  be  solved 
in  a  straightforward  manner  by  finding  the  homogeneous  and  particular 
solutions,  then  superposing  the  results,  and  evaluating  the  arbitrary  constants. 
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Particular  Solution 


F  U  (t)  Q 

Y  -  °  s 

P  “  3 

q 

Homogeneous  Solution 
The  roots  are  equal,  hence: 


(76) 


A^e 


-qt 


(77) 


T otal  Solution 


Y  =  A 1  e  ~qt  +  A2te"qt  + 


F  U  (t)Q 
o  s' 


3 

q 


(78) 


Initial  Conditions 


Y  =  0  a*  t  =  0 


dY 

dt 


0  at  t  =  0 


Thus 


A, 


Fous(t)Q 

"  "  * 

3 

q 


F  U  (t)  Q 
o  s 


2 

q 


us(t)  (  - 


The  Inverse  LaPlace  Transform  of  Eq.  (79)  is: 


-  I 


u3(0 
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(79) 


(80) 


"  mta.ni  lagged  - 


or 


6) 


_Y _ 

F 

o 


r~ 

U8(<) 


(t 


') 


(f  -  t  ) 

Us({  '  *o>  - —  -Us^  "'o'** 


/  2  \ 
A'  - V  t2' 


(81) 


In  addition  we  can  obtain  V  by  using  Eq.  69). 


~  -  Us<f>  H  +  Us<(  -*>[*-  <]-  Us>'  -  fo>[f  -  *ol 
o 

+  us<f  <82> 

Both  Y  and  V  are  shown  in  Fig.  43. 

3.1.3.  Exp<~  lential  Change  in  F  res  sure  Across  Jet 
It  is  apparent  from  an  earlier  report  (Ref.  1,  page  13)  that  the 
control  region  pressure  looks  like  the  step  response  of  a  simple  first-order 
system  with  a  two  unit  time  constant.  The  step  response  of  such  a  system  is 
a  simple  exponential  function  of  time  which  can  be  easily  generated  for  the 
purpose  of  separately  testing  the  jet  equations.  The  simple  first-order 
response  will  be  a  good  test  input  since  it  is  reasonable  to  expect  that  the 
bubble  pressure  will  always  be  describable  by  a  first  or  higher-order 
differential  equation.  The  jet  equations  are  solvable  with  this  type  of  input. 
The  step  response  as  derived  in  Sec.  3.  1. 2  is  the  special  case  where  the 
time  constant  of  the  forcing  function  is  zero. 
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Fig.  43.  Response  of  Jet  Equations  to  Step  and  Exponential  Inputs 


(83) 


The  forcing  function  is 


t 

F  =  F  U  (t)  (1  -  e  T  )  (U  {£)  -  U  (t  -  £)) 

os  S  3  O 


The  LaPlace  Transformed  jet  equation  is 


2q 


t 

—  +  q2Y  =  F  U  (t)  <1  -  e  T 
dt 


(84) 


Particular  Solution 


F  U  (t)  Q  F  U  (t)  Q 
os  o  s 


T 


The  homogeneous  solution  is  the  same  as  Eq.  (77) 


T otal  Solution 


(85) 


Y  = 


*  -qt 

e  "f 


A2te”qt  + 


FoUs(t)  Q 


3 

q 


Fous(t)  Q 


With  the  initial  conditions  as  before,  the  arbitrary  constants  are: 


A 


1 


,  1  ,2  j 

q(  -  -  q)  / 


A„  = 

c. 


(86) 


-  109  - 


r 

r'v 


The  solution  can  be  obtained  by  inserting  the  expressions  for  A  j  and 
in  Eq.  (86)  and  then  taking  the  inverse  LaPlace  transform.  The  resulting 
expression  for  Y/F  is  long  and  is  best  presented  in  parts  for  the  various 
possibilities  of  X  and  t.  The  velocity  V  is  also  included. 


t  ti  JS  i.\  I 


'  +  M  e 
o 


V  ( t  -  f  -  t ) 

—  =  U8(t)  (t  -  t)  -  (if  -  iQ)  +  re  _ 2 _ 

F  '  T 

o  r 


t  >  f 


us(t)  - 


Z  2 

/ 4  -  (/  - 1  r 

o 


/ 


O  \  (  ^  -t) 


-  / )  1  -  e  T  -  /  e  T 
’  o 


/  "_2_\  J<_2> 

=  Us(t)  <Q  -  T  i  1  -  a  T  e  T 


/  >  l  +  t 
o 


Both  Y  and  V  are  shown  in  Fig.  43 


3.  2.  Solution  of  the  Jet  Equations  Using  Various  Finite  Difference 
Approximations  to  tne  Space  Derivatives 

The  spa<_e  derivatives,  8V/8 /  and  BY/8/,  in  Eqs.  (68)  and  (69) 

can  be  approximated  at  a  fixed  point  space  by  a  finite  difference  expression 

ci  the  form 


-  Ill 


AX 

v  a.f.;  1  -  k  -  n 

i=  1  11 


The  coefficients  of  the  differentiation  formula,  Eq.  (88)  can  be 
obtained  from  almost  any  numerical  analysis  text,  (Ref.  30  and  31)  though 
often  not  directly.  A  digital  computer  program  has  been  written  which 
calculates  these  coefficients  for  any  n  and  all  k's. 

In  the  following  ccses,  the  space  increment,  Af,  is  unity.  The 
jet  is  divided  into  3  segments;  the  first  four  of  which  experience  the  driving 
force,  F,  and  the  last  four,  no  driving  force.  Both  V  and  Y  are  computed,  and, 
in  the  three  representative  cases,  compared  with  the  theoretical  response 
by  means  of  a  graph. 

Case  I  -  Fig.  44  -  Step  Input 


Derivative  formulae: 


=  Y  -  Y 
*k  *k-l 


■  vk  -  vk.i 


Comments : 


1.  Gives  correct  steady-state  for  V. 

2.  Does  not  give  the  correct  steady-state  for  Y. 

3.  Computed  response  leads  the  theoretical  response 
for  segments  5-8. 


Case  II  -  Step  Input 


Derivative  formulae : 


=  °- 5  Yk+ 1  -  5Yk.i 
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(92) 


dl 


0.  5  V 


k+  1 


0.  5V 


k-1 


T erminating  Conditions : 

Same  as  Eqs.  (89)  and  (90) 

Comments : 

1.  Very  little  damping.  Damping  is  introduced  only 
by  the  terminating  conditions. 

Case  III  «-  Step  Input 

Deri'  ative  formulae: 


df 


J-5Yk  -2Yk-l+  °’5Yk-2 


(93) 


d  ! 


1.5V,  -  2  V,  .  +  0.5V.  , 
k  k-1  k-2 


Starting  conditions: 


(94) 


Y_x  =  Yj  (95) 

V_j  =  Vj  (96) 


Comments : 

1.  Steady  state  is  not  correct  for  either  Y  or  V. 

2.  Computed  response  leads  the  theoretical  response 
for  segments  5-8. 

Case  IV  -  Step  Input 

Derivative  formulae: 


dY, 


d  l 


Y 

o 
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\ 


1. 5Y-  -  2Y .  +  0.  5Y 
2  1  o 


dYk  'S’ 

— —  =  2a  a,Y.  (97) 

di  i=0  * 


Same  for  V 

C  omments : 

Same  as  for  Case  EL. 

Case  V  -  Fig.  45)  -  Step  Input 
Derivative  formulae: 
for 

Yk  ,  Eq.  (93) 

for 

Vk,  Eq.  (90) 


Comments : 

1.  Gives  correct  steady  state  for  Y  and  V. 

2.  Computed  response  leads  for  segments  5-8. 

Case  VI  -  Fig.  46  -  Exponential  Input 
Derivative  formulae: 

Same  as  for  Case  V. 

Comments : 

Same  as  for  Case  V. 

Conclusions: 

Only  the  derivative  formulae  used  in  Case  V  give  the  correct 
steady  state.  The  derivative  formula  for  V  is  a  two-point  linear  type. 

It  is  reasonable  to  expect  that  the  formula  for  Y  should  be  a  three- 

point  parabolic  type.  However,  it  is  also  seen  that  the  computer  response 

is  too  slow,  especially  for  Y.  Lead  action  still  occurs  for  segments  5-8. 
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Fig.  45.  Computed  Jet  Equation  Response  -  Case  V  -  Step  Input 


Position 


Segment  No. 


Velocity 


■Vih'' 


i 


This  is  to  be  expected  because  a  pure  transport  delay  is  impossible 
to  approximate  with  a  finite  set  of  continuous  differential  equations. 

The  general  conclusion  is  that  none  of  these  approximations 
work  very  well.  A  good  approximation  should  at  least  produce  a 
correct  steady  state  and  a  transient  which  is  a  reasonable  facsimile  of 
the  theoretical  one. 

3.  Characteristic  Curves  of  Partial  Differential  Equations 

3.3.1.  Significance  of  Characteristic  Ct  .  r e ? 

Certain  types  of  partial  differential  equations  possess  what  are 
called  characteristic  curves.  (Refs.  29  and  32).  As  an  example, 
consider  the  wave  equation, 


V 


2  dzp 


(96: 


where  P  =  f(f )  and  P^  -  0  at  t  =  0.  The  :olution  is  known  to  be  of  the 
form 


P(f)t)  =  I 
2 

which  can  be  represented  as  in  Fig.  47.  The  value  of  P  at  point  C 
depends  only  on  the  values  of  P  at  A  and  D.  Furthermore,  every 
point  within  the  shaded  area  is  determined  by  a  pair  of  points  on  the 
I  axis  in  the  interval  AB.  The  two  lines,  AC  and  BC,  are  called 
"characteristic  curves".  Disturbances  at  A  and  B  prroagato  along 
these  curves  and  combine  to  form  the  value  of  P  at  point  C.  The 
shaded  area  is  called  the  "region  of  determination"  of  point  C. 

Another  region  of  determination,  { Ref.  32),  DEC,  can  be 
defined  in  r  ig.  48  when  Eq.  (98)  is  converted  into  a  difference  equation. 
The  value  of  P  at  point  C  is  now  determined  by  the  values  of  P  on  the 
last  two  horizontal  lines. 


f(f  -  Vt)  +  f(i  +  Vt) 


) 


(99) 
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If  the  solution  to  the  difference  equations  is  to  converge  to  the 
true  solution  then  the  region  of  determination  for  the  differ  ice  equations 
must  include  or  coincide  with  that  of  the  differential  equations.  Thus  DE 
must  be  greater  than  AB.  If  DE  were  smaller  than  AB  th  n  the  values  of 
P  at  points  A  and  B  would  not  influence  the  value  of  P  at  point  C;  in  which 
case  the  solution  does  not  converge  to  the  true  solution  when  the  mesh 
spacing  shrinks  to  zero. 

3.3.2.  Characteristic  Cu  ves  for  the  Jet  Equation 
The  characteristic  curves  »  ^q.  (72)  can  be  obtained  from 

the  quadratic:  (Ref.  29). 


ad!2  -  bdtdi  +  cdt2  =0  (100) 

where  a,  b,  and  c  correspond  to  the  first,  second,  and  third  coefficients 
of  Eq.  (72)  respectively.  With  a  =  c  =  1,  and  b  =  2,  Eq.  (100)  becomes: 

(d!  -  dt)2  =  0  (101) 

Hence, 

J*L=  l  (102) 

d! 

The  single  characteristic  curve  for  the  jet  equation  is  shown  in  Fig.  49. 

The  value  of  Y  at  B  is  determined  by  the  value  of  Y  an  equal  number  of 
space  and  time  units  previous.  With  F  =  0,  the  solution  to  Eq.  (72)  is: 

Y(f , t)  =  Y(f  -  A,  t  -  A)  (103) 

where  A  is  a  constant.  If  F  is  not  zero  it  would  appear  that  Eq.  (103)  is 
in  need  of  an  additional  term  on  its  right  side.  This  term  should  represent 
the  cumulative  effect  of  F  on  the  jet  over  the  intervals  (!  -  A)— *>!  and 
(t  -  a  )  *“•  t. 

All  of  the  above  indicates  that  the  jet  equation  should  be  solved 
by  a  technique  which  accumulates  the  effects  of  F  and  uses  the  values  of  Y 
only  at  previous  t’s  and  i’s. 
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Fig.  49.  Characteristic  Curve  for  Jet  Equation. 
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3.4.  Generalized  Transmission  Processes 
3.4.1.  Definition  of  a  Transmitter 
The  two  equations. 


(104) 


9q  _ 


(105) 


can  be  combined  to  obtain  the  well-known  wave  equation  in  terms  of  P  or  q: 


(106) 


2  2 
9  q  _  9  q 


(107) 


One  can  write  r.qs.  (104)  and  (105)  in  a  convenient  matrix  form  after 
defining  the  operator 


(108) 


Thus 


O  •  -D 

- H - 

-D  J  O 


(109) 


=  NS 


(110) 


The  variables,  P  and  q,  are  called  state  variables,  and  g  ,  a  state  vector. 
Processes  which  can  be  described  by  Eq.  (109)  are  w.-velike  or  diffusive  in 
nature.  Thus  Eq.  (110)  describes  the  generalized  transmission  process. 


(Ref.  33). 
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Equation  (104)  is  the  result  of  applying  the  momentum  equation. 
If  body  forces  are  present,  then  an  additional  term  is  necessary. 


dP 

dx 


9q 

at 


+  Fb  • 


(mi 


Equation  (105)  is  the  result  of  applying  the  continuity  equation.  Likewise, 
if  sources  or  leaks  are  present  then: 


9q_  _ 

a.x 


ap 

at 


+  Q  . 


Thus  Eq.  (1 10)  becomes 


^-2—  =  NS  +  F  • 

dx 


(112) 


(113) 


where  f  is  a  vector  containing  and  Q.  Thus  Eq.  (112)  describes  the 
generalized  transmiss  on  process  with  external  forcing. 

3.4.2.  Solution  of  the  Generalized  Transmitter  Differential 
Equations 

The  assumed  form  of  solution  to  Eq.  (113)  is: 


Sx  =  N  So  +  GFo  . 

where  M  and  Q  are  square  matrices.  Differentiating  Eq.  (114) 


(114) 


dS 


dM 


dG 


dx 


dx 


dx 


(115) 


and  substituting  into  Eq.  (113)  gives: 


dM 

dx 


C  + 
uo 


dG 

dx 


c  =  NMS  +  NGF  +  c  (116) 

•  o  o  o  ro 
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Collecting  terms; 


-  NM  )  SQ  +  —  -  N6  j  po  =  f  .  (H7) 

dx  J  \  dx  J 

If  p  is  assumed  constant  over  the  interval  of  space,  x,  (though  not  time) 
then 


d  M 


ax 


\ 

-NM  1  SD  + 


d_G 

dx 


-NG  -I  Ft  =  o 


(118) 


Since,  in  general,  and  FQ  are  not  equal  to  zei  then  it  must  be  that 


d  M 

dx 


-NM  =  o, 


and 


(119) 


d  G 

dx 


-NG  =  H. 


(120) 


The  solution  to  Eq.  (119)  is  the  well-known  exponential  given  by 


M  =  eN  X  M  • 


(121) 


The  solution  to  Eq.  (120)  contains  two  parts: 
Homogeneous  Solution 


a  J 


n 


G  o 


(122) 


Particular  Solution 


G  =  -  N 

p 


-l 


(123) 
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T otal  Solution 


G  =  eN  x  G  -  N  -1  (124) 

o 

The  matrices  M  and  G  can  be  determined  from  initial  conditions, 
o  o 

At  x  ■=  0,  it  is  obvious  that 

SjT  S0  •  <125> 

Therefore,  from  Eq.  (114)  it  follows  that 

S0=  M  SQ  +  GF0  ■  <IZ6> 

which  demands  that 

M  =  H  and  G  =0. 

Hence,  it  follows  that 

M  0  =  IT  and  G  0  =  N 

Equation  (114)  becomes 

Sx=  aN  *  SD+  (*N  *  -E)N  -1  F0  (129) 

which  is  a  familiar  result  (see  page  11,  Ref.  i). 

3.4.3.  Application  to  the  Jet  Equations 
Equations  (69)  and  (68)  can  be  rewritten  as 

dY 

—  =  -  DY  +  V, 
dx 

and 

_dV_  _  _  DV  +  F> 
dx 


(130’ 

(131) 


(127) 


(128) 
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The  matric  equation  is 


Thus 


f  -D 


-D  V 


(132) 


(133) 


K.  I  -xD  -xD 

N  x  !  e  i  xe 


0  * 


(134) 


N-1  =  - 

O 


-  -D. 


(135) 


(eNx-n)N-‘ 


1  . .  -xD, 

—  (1  -  e  ) 


. .  -xD.  -xD 
( 1  -  e  )  -  xe 


i  ,,  xD. 
—  (1  -  e  ) 


(136) 


The  solution  for  Y  and  V  are: 


_  ( -xD.  -xD  \ 

Y(x)  =  a-xD  (Y0  ♦  xVo)  ♦  f 1  -  *-  ->  -  ■SS—. J  F. 


(137) 


V(x)  =  e"xD  VQ  + 


I  -xD  \ 

II  ~e -  ]f. 

D  / 


(138) 
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Making  use  of  the  fact  that  e~  is  the  time  delay  operator  and  replacing 
x  by  Ax  the  above  equations  become. 


x-f  Ax,  t+  Ax 


■  (Y  f  AxV>x,t  + 


..  -AxD . 
(1  -  e  ) 


D 


Axe 


-AxD 


D 


F 


(139) 


V  =  V  + 

x+Ax,t+Ax  x,  t 


1  -  e 


-AxD 


D 


)  F’ 


(140) 


By  comparing  with  Eq.  (103)  it  is  seen  that  indeed  Eqs.  (139)  and  (140)  do 
represent  the  desired  form  of  the  solution.  The  right-most  portions  of 
Eqg.  (139)  and  (140)  represent  the  integrated  effects  of  F  over  the  space 
interval  x — •»  x  +  Ax  and  time  interval  t — •*- 1  +  Ax. 

Equations  (139)  and  (140)  are  of  a  mixed  form  involving  both 
difference  and  differential  (or  integral)  equations.  Both  types  of  equations 
are  easily  solved  by  means  of  a  digital  computer.  The  method  of  solution 
of  the  portions  of  Eqs.  (139)  and  (140)  due  to  F  is  discussed  in  the  next 
section. 

3.  5.  Monotone  Processes 
3.5.1.  Definition 

Monotone  dynamic  processes  are  best  described  by  Paynter, 

(Ref.  33,  pp.  280  -  286),  of  which  a  small  section  is  reproduced  below. 

"A  large  number  of  fluid,  thermal,  chemical  and 
other  industrial  and  organic  processes  are  characterized 
by  a  step  response  which  is  monotonic  nondecreasing  in 
time  as  indicated.  (See  Fig.  50) 

The  corresponding  frequency  response,  at  least 
for  most  continuous  processes,  would  have  a  non- 
increasing  amplitude  and  non-decreasing  phase  lag 
with  increasing  frequency  as  follows:  (See  Fig.  51) 

All  linear  systems  giving  to  such  response  can  be 
called  monotone  processes. 

Monotone  response  is  manifested  through: 

(a)  Time  delay  or  dead  time. 

(b)  Dispersion  or  rise  time. 

In  physical  processes,  time  delay  is  usually 
associated  with  propagation  or  transport  phenomena  as 
measured  by  the  ratio  of  travel  distance  to  propagation 
or  transport  velocity.  The  dispersion  in  any  process  can 
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Fig.  50.  Step  Response  of  a  Monotone 


Fig.  51.  Frequency  Response  of  a  Monotone. 


Fig.  52.  Step  Response  of  a  Simple  Lag. 


ultimately  be  attributed  to  the  lav*  of  increasing 
entropy,  whereby  the  distributed  resistances  in 
any  system  cause  an  attenuation  increasing  with 
frequency.  Such  scattering  action  is  reduced  by 
isolation  and  relaying  methods,  but  is  always 
present  to  some  degree. 

Oscillatory  processes,  characterized  by  the 
presence  of  complementary  energy  storage  elements, 
will  have  monotonic  response  w'henever  the  energy 
dissipated  per  cycle  becomes  sufficiently  large  com¬ 
pared  to  the  energy  stored  in  each  mode. " 

By  way  of  an  example,  it  is  seen  that  the  simple  first-order 
lag  is  a  monotone  and  is  dynamically  characterized  by  a  single  number, 
the  time  constant.  (See  Fig.  52). 

It  is  possible  to  determine  a  transfer  characteristic  of  the 
form:  (Ref.  34). 


G(D)  = 


_ 1 _ 

1  +  bjD  +  b^D  +  . .  . 


(141) 


which  will  adequately  represent  a  monotone  process  over  a  desired 
range  of  frequencies.  In  order  to  determine  the  coefficients  of  the 
polynomial  it  is  necessary  to  have  a  representative  step  response 
E,  of  the  process  from  which  can  be  calculated  the  time  moments  which 
are  defined  as: 


a 

o 


+  oo 

I 


CD 


dE  (t),  (usually  normalized  to  unity) 


*  +  CO 

ak=  ^  tkdE(t).  (142) 

-  oo 

The  moments  and  the  polynomial  coefficients  are  related  by  the  determinant: 
(Ref.  34). 
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If  the  step  response  is  the  result  of  a  transfer  characteristic 
of  the  form  given  by  Eq.  (141),  then  indeed  one  will  obtain  an  exact 
representation.  However,  if  the  step  response  exhibits  discontinuities 
then  probably  the  process  contains  time  delays  and  consequently  cannot 
be  exactly  represented  by  Eq.  (141).  The  resulting  representation  in 
this  case,  though  not  exact,  should  still  prove  adequate  for  most 
applications.  A  word  of  caution  is  appropriate  at  this  point.  It  may 
seem  that  a  better  approximation  is  always  possible  by  taking  more 
terms  in  the  polynomial.  However,  some  of  the  higher  order  coefficients 
(b^.),  may  be  negative,  resulting  in  an  unstable  representation,  which  is 
generally  useless. 


3.5.2.  Application  to  the  Jet  Equations 
Let  us  consider  the  contributions  to  Y  and  V  caused  by  F 
in  Eqs.  (139)  and  (140): 


Y  _  (1  -  e'AxD) 

G1  ~  ~ 

F  D2 


Axe 


-AxD 


(144) 


D 


and 


„  .  v  (1  -  e-ixD) 
-  - 

F  D 


(145) 


One  can  immediately  w rite  the  solutions  to  these  equations  for  a  step 
change  in  F : 
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F  2  3  2 


=  t  -  U  (t  -  Ax)  (t  -  Ax)  . 
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These  responses  are  shown  in  Fig.  53.  It  is  now  apparent  that  Gj  and 
represent  monotones  with  discontinuities.  Therefore,  it  is  not 
possible  to  find  exact  representations  in  the  form  of  Eq.  (141)  for 
Gj  and  G^. 

The  moments  for  G^  can  be  easily  calculated  with  the  aid  of 
Eq.  (142).  In  this  case  Ax  is  taken  as  unity.  They  are: 


a  =  1,  a,  = 
o  *  1 


»  •  •  *  *  ~ 


k  +  1 


(148) 


From  these  and  Eq.  (143)  it  is  now  possible  to  obtain  the  b^: 


b  = 
o 


bi  -■ 


b,  = 


(149) 


b4  = 


Therefore,  the  best  approximation  to  G-,  is 


G2,= 


1  +  —  D  +  —  D2 


(150) 
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Monotone  Approximations  for  Jet  Equations 


The  step  response  of  -3  compared  to  that  of  in  Fig.  33. 

The  steady-state  value  of  the  unit  step  response  of  is  1/2. 
It  is  usually  more  convenient  to  normalize  this  to  unity  for  the  purpose 
of  calculating  the  b^.  The  steady-state  gain  of  the  resulting  transfer 
characteristic  can  then  be  modified  very  easily. 

Proceeding  as  before,  one  finds  that  for  G-,  the  a,  are: 


a 


o 


2 

k  +  2 


and  that 


b  =  1 
o 


2 
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The  best  approximation  to  is 
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j  - 
1  ' 
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135 


D~ 


6480 


(153) 


where  the  "1/2"  is  the  correction  for  the  steady-state  gain.  The  step 
response  of  Gj/  for  the  2nd,  3rd,  and  4th  order  cases  is  compared  to 
tnat  of  Gj  in  Fig.  53. 
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The  approximations  that  ultimately  will  be  used  for  Gj  and 
depend  on  two  factors: 

1.  The  expected  frequency  content  of  F,  and 

2.  Computer  etorage  and  time. 

It  was  seen  in  Sec.  3.  1.3,  that  the  bubble  region  pressure  looks 
like  the  response  of  a  simple  first-order  system  with  a  time  constant  of 
about  2  units.  The  bubble  region  pressure  is  the  force,  F,  in  Eqs.  (144) 
and  (145).  The  harmonic  content  of  the  bubble  pressure  signal  is  therefore 
small  at  frequencies  above  the  break  frequency,  1  /t  ^  0.5  rad/unit.  The 
natural  frequency  .f  G,1,  Eq.  (150),  for  example,  iso  =  3.46  rad/unit 
which  is  comfortably  above  0.  5.  Since  G^’  is  a  good  approximation  to 
G^  up  to  frequencies  of  3,46  rad/unit  then  it  follows  that  it  also  should  be 
sufficiently  accurate  when  driven  by  the  bubble  pressure.  Similar  con¬ 
siderations  apply  to  Gj'. 

Using  higher-order  approximations  for  G,  and  G^  without 
justification  is  wasteful  of  computer  storage  and  time.  If  a  jet  is 
modeled  by  8  segments,  a  2nd-order  approximation  for  G^,  and  a 
3rd-order  approximation  for  G^,  then  it  becomes  necessary  to  operate 
on  40  x  40  matrices.  A  minimum  number  of  three  such  matrices 
necessitates  4,800  w'ords  of  computer  storage. 

3.6.  A  Possible  Model  for  the  Fluid  Jet  Modulator 

From  the  results  of  the  preceding  sections  it  is  now  possible 
to  propose  a  simplified  model  for  a  single  control  port  fluid  jet  modulator. 
A  non -functional  block  diagram  is  shown  in  Fig.  54. 

The  input  to  the  system  is  the  control  pressure  applied  upstream 
of  the  control  port  restriction.  The  resulting  control  flow  combines  with 
the  flow  being  demanded  of  the  bubble  region  by  the  jet.  The  excess  flow, 
resulting  from  this  combination,  then  causes  the  bubble  region  pressure 
to  "charge-up"  to  a  new  level  because  of  the  compliance  (compressibility) 
of  the  fluid  in  the  bubble. 

The  bubble  pressure  then  causes  the  jet  to  accelerate  to  a  new 
position.  The  change  in  position  of  the  jet  at  the  knife-edge  affects  the 
demand  characteristics.  The  characteristics  are  also  influenced  by  the 
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bubble  pressure.  The  resulting  change  in  the  demand  flow  then  affects 
the  bubble  region  pressure  thus  closing  the  loop  in  the  diagram. 

Figure  55  shows  a  functional  diagram  of  the  same  system.  The 
checked  enclosures  indicate  the  correspondence  with  Fig.  54. 

The  difference  between  the  control  and  bubble  pressures  is 
multiplied  by  the  control  port  conductance,  G,  to  produce  the  control 
flow,  W  .  This  flow  is  being  supplied  to  the  bubble  region.  The  flow 
being  demai?ded  of  the  bubble  region  is  then  subtracted  off.  When  the 
system  is  in  steady  state  this  difference  is  zero.  An  excess  of  one  over 
the  other  causes  a  rate  of  change  of  bubble  pressure,  and  consequently, 
a  transient.  This  pressure  rate  is  integrated  to  produce  the  bubble 
pressure.  The  bubble  pressure  is  fed  back  thus  closing  a  loop. 

The  net  pressure  across  the  jet,  (P  -  P  )  is  applied  to  a  transfer 

2L 

characteristic  similar  to  Eq.  (153)  which  determines  the  position  of  the  jet 
in  the  region  of  the  knife  edge. 

The  atmospheric  flow  is  computed  as  a  product  of  two  distinct 
functions.  It  is  -postulated  that  this  flow  passes  through  some  effective 
area  which  is  a  function  of  Y  and  the  lateral  setback  of  the  knife  edge  and 
is  described  by  a  static  function,  0 The  effective  pressure  driving  this 
flow  is  a  function  of  (P  -  F  ),  or  its  negative,  and  is  described  by  0, .  The 
simplest  possible  functions  that  seem  appropriate  are: 

*  <Yo  -  Ye  -  ~  > 

where  Yq  i;»  the  lateral  setback  of  the  knife  edge  and  b  is  jet  "width  at  the 
knife  edge.  The  atmospheric  flow  would  be 

Wa  *  'Yo  -  Ye  '  J>f<Pa  ‘  P> 

which  is  recognized  as  the  form  of  the  simple  orifice  equation. 
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Fig.  55.  Functional  Block  Diagram  of  a  Pressure  Controlled  Fluid  Jet  Modulator. 


The  return  flow  it*  computed  as  a  function  of  the  jet  position  only. 
The  simplest  possible  reticulation  of  0^  would  be 

W  =  0;  Y+-+6^Y 

r  e  2  o 

W  d(Y  +  ^  +  6-  Y);  Y  +  —  +  6  -  Y 

r  e  2  o  e  ^  o 

where  6  is  a  constant  which  allows  the  =  0  point  to  fall  within  or 
without  the  non-zero  range  of 

The  entrained  flow,  Wg,  is  assumed  constant  although  this  is 
probably  not  strictly  true. 

The  combination  of  these  three  flows  then  combine  to  produce 
the  demand  flow  thus  closing  the  loop. 

The  above  model  assumes  that  the  jet  equations  derived  in 
Sec.  3.  5.2.  are  valid  (at  least  dynamically)  and  that  the  demand 
characteristics  can  be  reticulated  as  indicated  in  this  section.  Both 
assumptions  must  be  verified  by  experiment. 
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4.  THE  DYNAMICS  OF  SYSTEMS 


4*  1.  The  Transmission  of  Waves  in  T  nes 

Viscosity  and  heat  transfer  disperse  the  wave -fronts  in  fluid 
lmes.  LaPlace  operators  for  the  propagation  and  characteristic  im¬ 
pedance  functions  have  been  derived  by  the  author  (Ref.  24)  for  liquid 
and  perfect-gas  filled  lines,  and  the  frequency  response  and  parts  of 
the  step  responses  were  calculated.  Following  the  suggestions  of  the 
author  (Ref.  35),  S.  E.  Nelson  found  the  complete  families  of  the  step 
responses  of  liquid-filled  lines.  This  work  is  presented  in  a  thesis  (Ref. 

36);  a  professional  paper  is  planned.  All  of  the  above  work  assumed 
small  disturbances  in  rigid  cylindrical  lines  with  laminar  flow. 

D.  M.  Auslander  is  presently  finishing  his  S.M.  thesis  in 
which  he  makes  theoretical  predictions  for  the  frequency  responses  of 
liquid  flow  in  cylindrical  rigid  lines  with  a  gross  turbulent  flow.  The 
The  predicted  attenaution  is  much  greater  than  for  laminar  flow.  Both 
smooth -walled  and  rough-walled  pipes  are  considered. 

None  of  the  above  work  was  supported  by  the  present  contract. 

A  question  of  frequent  importance  is  the  effect  of  bends  on  the 
transmission  of  waves  in  lines.  A  rather  simple  answer  can  be  provided 
when  the  wave-lengths  of  interest  are  much  longer  than  the  dimensions  of 
the  bend,  and  the  disturbance  is  small.  Under  these  conditions  the  steady- 
flow  characteristics  of  the  bend  apply,  producing  a  pressure  drop  propor¬ 
tional  to  the  square  of  the  flow.  The  proportionality  factor,  which  depends  on 
the  sharpness  of  the  bend,  is  given  for  several  geometries  in  even  elementary 
fluid-mechanics  texts.  The  small  disturbance  assumption  allows  ar.  incre¬ 
mental  linear  resistance  to  be  found.  This  resistance  is  an  increasing 
function  of  the  gross  flow;  it  is  zero  for  zero  gross  flow. 

Thus  the  bend  acts  like  a  friction  joint  in  a  straight  pipe.  An 
incident  wave  is  partly  reflected  and  partly  transmitted;  ii  R  is  the  linearized 
resistance  for  small  disturbances  about  a  gross  flow  in  a  line  with  charac¬ 
teristic  impedance  Zq  (roughly  equal  to  the  speed  of  sound  divided  by  the 
area  of  the  line),  the  reflection  and  refraction  operators  (coefficients  in 
simple  cases)  are  respectively, 
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reflection  operator: 


refraction  operator : 
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2  Z  +  R 
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(154) 


(155) 


4.  2.  The  Stability  of  Fluid  Systems 

Fluid  systems,  especially  those  containing  proportional  fluid 
amplifiers  with  no  moving  parts,  often  display  little  understood  bistable 
or  cyclic  instabilities.  The  author  has  prepared  a  paper  on  this  subject 
(Ref.  37)  which  will  be  widely  available  before  the  present  report;  the 
present  contract  is  therein  credited  with  the  support  for  the  work.  Therefore, 
the  substance  of  this  work  omitted  herein. 

Two  general  criteria  for  small -disturbance  stability  are  presented, 
with  simple  illustrations  involving  amplifiers,  lines,  and  volumes.  Applica¬ 
tion  of  the  criteria  to  real  systems  is  based  on  separate  measurements  of 
certain  static  and  dynamic  properties  of  the  more  complicated  elements 
and  analysis  of  the  separate  properties  of  the  other  elements. 
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S.  STALL  Ci^LL  SlNSOk 


5.1.  Introduction 

The  concept  and  test  results  on  the  first  model  of  a  stall-cel’ 
sensor  are  described  in  Sec.  3.  2  of  the  first  annual  report.  Analytical 
work  of  both  a  theoretical  and  experimental  nature  has  continued.  A 
static  stability  analysis  has  been  completed  which  should  result  in  useful 
design  formulae.  The  experimental  study  was  continued,  utilizing  a  new 
test  model,  and  although  results  are  mostly  of  qualitative  value,  some 
quantitative  data  is  presented. 

5.  Sizable  Jet  Amplification 

Fluid  jet  amplifiers  are  classed  as  analog  (proportional)  and 
digital  (bi-and  tri -stable)  depending  on  the  stable  operating  positions 
of  the  jet.  A  bistable  jet  is  unstable  in  the  center  position  and  stable 
only  in  its  two  extreme  positions.  The  mechanism  of  bistable  jet 
amplification  is  explained  as  follows. 

While  a  free  jet  is  inherently  statically  stable,  a  jet  which  is 
not  permitted  to  blow  freely  into  the  environment  may  be  rendered 
statically  unstable  through  proper  design  of  the  discharge  region.  Con¬ 
sider  a  jet  discharging  to  a  region  such  as  that  shown  in  Fig.  56.  In 
passing  through  the  discharge  region,  the  jet  entrains  flow,  we»  The  flow 
entrained  is  relatively  unaffected  by  control  port  pressure  differences  and 
jet  curvature  and  hence  remains  constant  on  each  side.  The  entrained 
flow  removed  fromthe  control  port  is  equal  to  the  flow  into  this  region 
from  the  environment,  w^.  Then,  on  the  side  in  the  direction  of  jet 
curvature,  the  velocity  of  the  flow  passing  through  the  area  between  the 
projection  and  the  edge  of  the  jet  increases.  Due  to  increased  frictional 
losses,  the  pressure  on  that  side  (Pj)  decreases.  The  pressure  on  the 
opposite  side  (P2)  increases.  A  pressure  force  is  then  seen  to  act  on  the 
jet  in  the  direction  of  curvature.  If  the  pressure  force  accompanying  a 
''small"  deflection  of  the  jet  exceeds  the  force  due  to  the  accompanying 
Y -component  of  momentum  reaction,  static  instability  results  -  the  jet 
curvature  will  then  increase  until  the  jet  attaches  to  the  projecting  edge. 

In  the  case  of  a  two-sided  amplifier,  such  as  the  one  shown,  this  is  re¬ 
ferred  to  as  bistable  operation. 
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Fig.  56.  Discharge  Region  Giving  Static  Instability. 


W.  Y  =  lateral  setback 

J  ° 


R  =  control  line  resistance 

b  -  nozzle  width 
o 


Discharge  Region  With  Control  Lines  Added. 


Curvature  of  a  bistable  jet  can  be  effected  by  adding  control  lines 
to  each  side  of  the  discharge  region.  This  is  shown  in  Fig.  57  along  with 
nomenclature  that  will  oe  of  use.  By  increasing  the  control  flow  on  the 
side  of  attachment,  decreasing  control  flow  on  the  unattached  side,  or 
both,  the  jet  can  be  flipped  from  one  side  to  the  other.  The  ability  of 
the  power  jet  to  flip  quickly  in  response  to  an  alternating  signal  trans¬ 
mitted  to  the  control  port  led  to  the  proposed  application.  A  schematic 
of  the  stall  sensing  arrangement  is  shown  in  Fig.  58  (see  also  Fig.  20 
of  Ref.  38).  In  Fig.  59c  a  possible  packaging  arrangement  is  shown. 

5.  3.  Theoretical  Study 

5.3.1.  Static  Stability  Analysis 

Brown  has  presented  a  criterion  for  static  and  surge  stability  of 
the  center  position  of  the  power  jet,  based  on  experimental  measures  of 
the  characteristics  of  the  device  as  measured  at  the  control  port-;.  Van 
Koevering  has  generalized  this  to  include  non-center  jet  positions  and  con¬ 
figurational  asymmetries.  Brown  and  Simson  have  reported  criteria  for 
the  wave  (organ  pipe)  stability  using  some  of  the  same  results,  and  in 
Sec.  4.  2.  of  this  report.  Brown  discusses  a  generalized  approach  to  the 
stability  and  response  of  linear  systems. 

The  stability  analysis  for  the  center  position  of  the  jet  is  herein 
interpreted  analytically  in  terms  of  the  jet  model  given  by  Simson  (Ref.  39) 
plus  one  added  assumption,  reducing  greatly  the  requirement  fo-  e.cperi- 
mental  data.  The  analysis  is  for  a  symmetrical  amplifier,  assumes  no 
feedback  effects  of  the  receiver  ports,  and  neglects  the  dynamic  head  of 
the  control  flow  which  enters  through  a  very  large  passage.  In  addition, 
the  "return  flow,  w^",  is  assumed  zero  for  the  center  position  of  a  jet, 
a  condition  which  is  approximately  necessary  i  >r  bistability.  The  flow 
entrained  into  the  jet,  wg,  is  assumed  to  be  unaffected  by  small  deflections 
of  the  jet.  Thu3  the  changes  in  the  control  flow  wg  are  equal  to  minus  the 
changes  in  the  flow  wa  which  enters  the  control  region  from  the  environment 
by  passing  between  the  knife  edge  and  the  jet.  Brown's  criterion  for 
center -position  stability  then  becomes  (See  Ref.  40) 
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in  which  the  subscripts  1  and  2  refer  to  the  two  control  regions,  is 
the  deflection  of  the  jet  as  shown  in  Fig.  60,  and  R.  is  the  resistance  of 
the  upstream  restriction.  The  resistance  R  can  be  assigned  either  steady 
state  or  surge  values. 


The  deflection  of  the  jet,  for  small  pressure  gradients,  is  very 


nearly 


Yi=  - 


X  “  (P  -  P.) 
o  a  1 ' 


(1  57) 


2  p  b  V. 
°  J 


an  assumption  based  on  the  uniformity  of  the  pressure  gradient  across 

the  jet  and  the  conservation  of  the  momentum  of  the  jet.  The  nozzle  width 

is  b  ,  and  the  jet  mean  velocitv  V.. 
o’  J  '  J 

The  most  critical  assumption  is  the  modelling  of  the  flow  w 


=  Cd  YallV 


(158) 


where  t  is  the  depth  of  the  model.  This  is,  of  course,  the  classical  adaption 
of  Bernoulli's  equation  to  flow  through  an  orifice.  The  orifice  here  is  as¬ 
sumed  bounded  by  the  "edge"  of  the  jet.  The  atmospheric  flow  is  assumed 
to  be  into  the  control  port,  a  condition  necessary  for  bistability.  For  positive 
control  flow,  this  imposes  the  requirement  w  <  w  . 

C  c 

The  substitution  of  Eqs.  (157)  and  (158)  into  Eq.  (156),  noting  that 
dYa^=  -  dYj,  gives  the  condition  for  stability 


>  -  1  + 
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r  '  d  ao 
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o  d  ao  j 
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The  remaining  steps  are  to  put  w  and  Y  into  more  convenient  forma. 

a  o  au 

Using  Simson's  model  for  a  jet,  and  assuming  that  the  knife  edge  is  within 
a  distance  of  5.2  nozzle  widths  from  the  nozzle, 


Y  =  Y  -  0.  169  X  -  0.  5  b 
ao  o  o  o 


(160) 


The  flow  w  is  the  difference  between  the  entrained  and  control  flows; 
ao 

using  Simson's  model  for  the  former  and  for  the  latter. 


w  ~  0.035  pX  t  V.  -  w 
ao  r  o  j  co 


(161' 


For  the  general  case  of  finite  resistance  and  Eq.  (159) 
with  the  condition  of  Eq.  (161)  results  in  the  following  criteria  for  bistable 
operation: 
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for  square-law  resistance. 


Yj=  0  Z(Pc  -  P) 


where  Pc  is  the  control  line  supply  pressure,  P  is  the  average  control 
port  pressure  (  =  P.  +  P-)/2),  and  Y  is  given  by  Eq.  (J60).  These 
results  enable  a  theoretical  prediction,  for  example,  of  the  optimum  control 
line  resistance  for  a  particular  configuration. 

For  the  special  case  of  blocked  control  passages  (R  =  oo),  Eqs. 
(160),  (16Z),  and  (163)  reduce  to 


Y  <  0.  169  X  +  0.  5b  + 
o  o  o 


(0.035)2  Xq4 

C  2  b 
d  o 


(164) 


For  the  case  of  finite  R,  this  relation  represents  a  necessary  (although  not 
sufficient)  condition  for  bistability  on  which  dimensional  choices  can  be 


based. 


The  upper  limit  imposed  on  Yq  by  Eq,  (164)  represents  the 


transition  from  bistable  to  tristable  operation.  By  reducing  the  control 
line  resistance  sufficiently,  the  operation  can  then  become  continuous.  The 
lower  bistability  limit  is  approximately  the  point  where  the  return  flow 
becomes  finite  for  a  center -positioned  jet.  This  limit  represents  a  transition 
from  bistable  to  continuous  operation,  and  may  be  estimated  from  Eq.  (160), 
by  setting  Yao  =  0.  The  dimensional  limits  for  bistability  given  by  Eqs.  (160) 
and  (164)  are  plotted  in  Fig.  61. 

5.3.2.  Atmospheric  Flow  Coefficient 

Continuing  the  analogy  of  orifice -type  flow  for  the  atmospheric  flow 
component,  it  is  shown  by  dimensional  analysis  that  the  loss  coefficient  C^ 
has  the  following  functional  relationship  (neglecting  the  effect  of  variations 
in  control  port  shape) 
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w  =  0.  035  X  t  1  2  p  (P  -  P)  -  w  (1 61) 

so  that  the  desired  functional  relation  for  is  determined  experimentally 

by  measuring  P  ,  w  and  P. 

3  6  s’  CO 

Da.a  of  this  type  has  been  extracted  from  a  load-flow  diagram 

presented  by  VanKoevering  (Ref.  41,  Pig.  5.8+)  and  is  presented 

graphically  in  Fig.  62.  This  plot  is  for  fixed  Y  /t;  is  seen  to  vary 

moderately  with  Reynolds  number  over  this  range.  The  latest  test  model 

(described  in  Sec.  5.4.  1)  was  tested  with  Yao/t  ranging  from  0.  39  to  0.  58 

and  Reynold's  number  on  the  order  of  4000.  had  values  ranging  from 

0.  2  to  0.3  and  was  a  moderate  function  of  Y  /l  and  weak  function  of 

ao' 

Reynold's  number.  Leakage  resulting  in  a  large  control  port  pressure 
asymmetry  limited  the  quantitative  value  of  the  data. 

5.  3.  3.  General  Application  of  Theory 

The  procedure  for  applying  the  proposed  theory  in  designing  for 
optimum  bistable  operation  is  outlined  below.  This  is  for  the  most  general 
situation,  where  is  a  function  of  Reynold's  number  in  the  range  of  interest. 

+  The  value  of  R'  on  VanKoevering 's  Fig.  5.8  should  be  6.75,  The  corrected 
value  is  computed  from  p.  98  of  his  thesis. 

-  152  - 
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Fig.  62.  Atmospheric  Flow  Coefficient  from  VanKoevering 
Data. 
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1.  Determine  Y^/t.  Approximating  P  =  WCQ  =  0, 
determine  Re  from  Eqs.  (166)  and  (161),  the  control 

3. 

line  resistance  being  known. 

2.  From  experimental  data  of  the  type  presented  in 
Fig.  62,  determine  C^. 

3.  Solve  Eqs.  (162)  or  (163)  for  the  control  line 
resistance,  still  approximating  P  =  0. 

4.  From  Eqs.  (1  58),  (160),  and  (161)  solve  for  wco 
and  P. 

5.  Using  these  values  of  w  and  P  determine  the  new  Re 

&  co  a 

and  corresponding  C^.  Solve  Eq.  (162)  or  (163)  again,  etc. 
until  the  process  converges.  Better  initial  assumptions 
should  lead  to  faster  convergence. 

5.  3.4.  Application  of  Theory  to  Stall  Sensing  Amplifier 
For  the  application  at  hand,  with  square -law’  control  line  resistance 
and  equal  supply  and  control  port  pressures,  Eq.  (162)  simplifies  to 


0.  035  X  t  (b  Y  +  2X  ) 
o  o  ao  o 

c  TL  (b  y  +  x  2) 

o  .  o  ao  o 
4 


i —  _  — r 

t2  b  ci  Y  -  (0.035)2  X  4 
ro  d  ao  '  o  __ 


C  I  WbY  +  X2 

o  4  I  I  o  ao  o 


(168) 

where  d  and  C  are  the  control  line  orifice  diameter  and  flow  coefficient, 
o  o 

respectively.  For  control  flow  Reynold's  number 

2vd 

Rc  =  - —  >  260 

v 

where  v  is  the  average  velocity  through  a  sharp-edged  orifice,  C  can 
be  assumed  constant  at  about  0.  625  (Ref.  4  2,  p.  183).  This  Reynold's  is 
easily  exceeded  in  the  present  application. 


Brown's  graphical  source -load  matching  technique  was  applied 
to  the  load  flow'  curve  corresponding  to  Fig.  62.  A  value  for  optimum 
control  line  orifice  diameter  of  0.031V  inch  was  predicted  at  en  operating 
point  for  which  the  atmospheric  flow  Reynolds  number  is  determined  to 
be  approximately  105.  From  Fig.  62  this  corresponds  to  a  flow  coefficient 
of  0.66.  The  appropriate  dimensions  and  dQ  =  0.0319  inch  in  Eq.  (168) 
give  =  0.  566  in  fair  agreement.  The  variation  in  d^  for  these  two 
values  of  is  about  30  per  cent. 

Using  the  dimensions  of  the  new  test  model;  Xq  =  0.  253  inch, 

bQ  =  0.  058  inch,  t  =  0.  1 75  inch,  Eq.  (168)  was  solved  for  dQ  in  terms  of 

Y  '  Y  /b  ).  Two  positive  real  roots  were  alwa.ys  obtained  but  a  con- 
o  v  o'  o  r  7 

sideration  of  the  direction  of  inequality  ruled  out  the  larger  (as  they 
indicated  that  orifices  larger  than  the  optimum  favor  bistability).  The 
results  are  plotted  in  Fig.  63  (flft  =  Aspect  Ratio  =  b/t). 

5.4.  Experimental  Study 

5.4.  1.  Test  Apparatus 

The  test  amplifier,  shown  in  Fig.  64,  is  a  modification  of 
several  earlier  models  which  presented  difficulties  due  to  (1)  Control 
port  resonance,  (2)  Knife-edge  reattachment,  (3)  Pressure  asymmetry. 

The  first  problem  arises  from  wave  phenomena  in  the  control 
regions.  It  is  in  part  dv  to  low  frequency  standing  waves  and  the  effect 
is  best  minimized  by  making  the  control  ports  as  small  as  possible. 

Small  control  ports  also  favor  dynamic  response. 

In  order  to  minimize  hysteresis,  sharp  knife  edges  are  employed. 
This  presents  a  problem  due  to  the  formation  of  a  separation  pocket 
downstream  of  the  knife  edges  (see  Fig.  65).  A  deficiency  of  flow  from 
the  environment  into  this  region  can  cause  the  jet  to  reattach  downstream 
of  the  edge  of  the  knife  resulting  in  a  large  increase  in  hysteresis.  Data 
summerized  by  Simson  (Ref.  39,  p.  131)  predicted  reattachment  at  larger 
lateral  setbacks  (hence  small  (3  -  see  Fig.  65)  with  the  present  test  model 
for  all  but  impractical  knife-edge  lengths.  An  alternative  solution  involved 
placing  large  holes  in  the  face  plates  in  the  region  of  the  separation  pockets. 
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Optimum  orifice  diameter 


Fig.  64.  Test  Model  With  Receiver  lubes 
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from  « rn nil  dimensions!  m is s  1 1  g e n t s  s nci  lesk-B.ge.  XKc  extreme  sens! — 
tivity  of  setback  tolerance  on  op  .  ating  characteristics  led  to  a  micrometer 
adjustment  permitting  easy  changes  to  0.001  inch  accuracy.  The  problem 
of  leakage  was  particularly  perplexing  because  the  desired  free  movement 
of  .he  knife  edges  prevented  adequate  sealing. 


The  dynamic  test  arrangement  (Figs.  66  an  1  67)  is  similar  to 
that  of  the  first  annual  report.  To  eliminate  control  line  dynamics,  t^e 
control  lines  are  almost  negligibly  short,  although  they  were  initially  long 
enough  to  insert  piezoelectric  pressure  transducers.  During  static  testing 
and  while  adjusting  the  signal  line  supply  pressure  for  dynamic  testing,  a 
monometer  is  added  to  the  signal  control  line.  The  test  model  has  dimensions 
given  in  Sec.  5,  3.4. 

5.4.2.  Static  Test  Results 

A  large  amount  of  experimental  data  has  been  taken  using  the 
present  and  earlier  test  models.  Unfortunately,  dimensional  misalignment 
and  leakage  problems  limited  the  quantitative  value  of  most  of  the  results. 


The  upper  instability  limit  (for  R  ~  oo)  as  given  by  Eq.  (164) 
is  plotted  in  Fig.  61.  Setbacks  are  nondimensionalized  with  respect  to 

X 

nozzle  width.  The  results  of  two  test  runs  at  a  supply  pressure  of  3.45  psia ' 
and  with  Xo*  =  3.89  are  shown  in  Fig.  61.  The  smallest  possible  setback 
used  in  practice  to  decrease  the  hysteresis  bond  of  the  control,  increasing 
the  gains. 

From  Eq.  (168),  dQ  is  plotted  versus  Y  1  and  in  Fig.  63,  for 
the  dimensions  of  the  present  test  model.  The  experimentally  determined 
optimum  orifice  diameters  for  Y0'  =  2.07  are  plotted.  Agreement  's  fair 
within  the  range  of  validity  of  the  theory  (i.e.,  incompressible  flc  .  -  P  <  2  in,  Hg). 
Data  points  further  up  in  the  compressible  range  (6  to  10  in.  Hg)  are  off  the  graph, 
A  sudden  jump  between  4  and  5  in  Hg  was  observed. 


3.45  psi  is  the  average  stall-free  pressure  ratio  across  a  single  compressor 
stage  as  determined  in  the  first  annual  report. 
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nic  Test  Set  Up 


Fig.  67.  Dynamic  Test  Apparatus. 
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One  of  the  principal  objectives  of  the  experimental  study  was  to 
determine  if  the  hysteresis  could  be  reduced  sufficiently  to  permit  bi¬ 
stable  operation  with  the  available  pressure  swing  (PIvIAX  -  PMIN  <0.1  pai). 
This  was  done  successfully  and  reported  in  the  First  Annual  Report 
(see  Fig,  21).  With  the  new  test  model,  operation  with  very  small 
hysteresis  was  difficult  due  to  control  port  noise  of  the  same  order  of 
magnitude  as  the  hysteresis.  The  signal  (slightly  filtered)  from  a 
piezoelectric  in  one  control  port  is  shown  in  Fig.  68.  Control  port 
resonance  caused  random  flipping  of  the  jet,  the  large  jump  indicating 
jet  movement.  The  increase  in  noise  is  possibly  due  to  the  sharp  knife 
edge.  Tilting  the  knives  at  an  angle  less  than  90°  would  probably  reduce 
the  noise. 

Defining  Gain  as  Supply  Pressure  times  Recovery  Ratio  divided 
by  Hysteresis;  with  a  recovery  of  70  per  cent,  hysteresis  within  0.  1  psi, 
supply  pressure  of  3.4  psi,  the  gain  would  be  slightly  greater  than  24. 

5.4.3.  Dynamic  Test  Results 

A  primary  purpose  of  the  experimental  study  was  to  determine 
the  frequency  response  characteristics  of  the  amplifier.  In  Sec.  3.2.1 
of  the  First  Annual  Report,  it  was  established  that  the  amplifier  must 
have  a  frequency  response  better  than  400  cps  at  5  per  cent  modulation, 
corresponding  to  a  rise  time  of  0.  125  millisecond.  The  test  results 
reported  in  the  first  report  indicated  a  maximum  frequency  response 
of  260  cps  at  37  per  cent  modulation  and  only  60  cps  at  8  per  cent 
modulation.  This  corresponds  to  a  rise  time  of  about  1.4  milliseconds. 

For  the  nozzle  supply  pressure  of  2.94  (pressure  ratio  =  1.2)  the  nozzle 
flow  velocity  is  598  ft/sec.  Then  with  the  lateral  setback  -  0.06  inch, 
the  transport  time  defined  as  2X  /V  is  0.0167  millesecc  d.  Thus  the 
response  time  is  greater  than  80  transport  times  which  Is  much  larger 
than  the  results  of  Sec.  2.  1. 5.  in  the  first  report  would  indicate.  The 
discrepancy  was  attributed  to  dynamics  associated  with  the  long  control 
line. 

In  the  new'  test  model  the  control  line  was  virtually  eliminated. 

The  control  ports  however,  were  larger  than  before  to  permit  instrumenting 
each  control  port  with  a  piezoelectric  pressure  transducer. 
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Fig.  68.  Random  Flipping  of  Jet  Due  to  Control  Port 
Noise. 
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The  stall  cell  signals  are  shown  in  Fig.  20b  of  the  first  report. 
Xhc  of  rTiCdulstiiis  ths  d^d  no^  f'6 rmit  ftx3ir^  s irnul— 

of  the  stall  signal .  the  minimum  nreasure  beine  lower  than  in  the  actual 

C?  •  4  V 

case.  The  experimental  Input  signal  is  shown  in  Fig.  69.  The  rise  time 
associated  with  the  signal  increase  should  closely  simulate  the  rise  time 
with  the  actual  signal. 


Dynamic  test  results  are  shown  in  Fig.  70.  The  lower  trace  is 
from  the  piezoelectric  in  the  signal  control  port  and  the  upper  trace 
from  a  piezoelectric  placed  downstream  of  the  knife-edge  as  in  the  first 
report.  The  modulation  is  approximately  25  per  cent. 

The  asymmetry  in  the  input  signal  is  evident  in  the  lower  trace, 
the  rise  time  in  one  direction  being  much  greater  than  in  the  other.  The 
effect  on  the  upper  trace  is  to  lengthen  the  output  pulse  width.  At  455  cps, 
the  control  port  transducer  indicates  sharp  response  to  the  larger  de¬ 
creasing  pressure  signal. 

5. 5,  Conclusions 

Experimental  data  enabling  an  assessment  cf  the  theoretical 
stability  criteria  is  incomplete.  To  avoid  the  dimensional  asymmetry 
and  leakage  problems  a  large  scale  test  model  should  be  constructed. 

With  such  a  model  these  extraneous  effects  would  have  less  consequence. 

The  dynamic  test  results  demonstrate  the  strong  dependence  cf 
response  on  pressure  swing.  A  reduction  in  size  of  the  control  ports 
could  increase  the  response  significantly.  With  the  larger  pressure 
swing  it  appears  that  the  response  is  adequate.  By  staging  amplifiers, 
using  proportional  devices  in  all  but  the  last  stage,  the  stall  cell  signal 
could  be  amplified  to  the  necessary  pressure  swing. 
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Fig.  69.  Pulse -Width-Modulator  Input  Signal. 


6„  DESIGN  OF  A  TEMPERATURE -INSENSITIVE  PNEUMATIC 

OSCILLATOR 


6.1.  Introduction 

A  constant -frequency  pneumatic  oscillator  is  a  necessary  element 
in  pneumatic  logic  systems  and  could  form  part  of  the  basis  of  temperature- 
measuring  and  speed-measuring  devices.  The  frequency  of  a  really  useful 
oscillator  should  be  insensitive  to  the  temperature  and  pressure  of  the  fluid. 
The  objective  of  this  report  has  been  the  development  of  a  fixed-frequency 
pneumatic  oscillator  which  requires  n-  temperature  control  and  only  the 
pressure  control  of  a  common  regulator. 

A  very  simple  oscillator  consists  of  a  pneumatic  jet  amplifier 
with  transmission  lines  connected  to  its  control  ports  in  such  a  way  as 
to  cause  a  dynamic  instability.  See  Fig.  71.  Since  the  velocity  of  sound 
varies  as  the  square  root  of  the  absolute  temperature,  the  lengths  of  the 
transmission  lines  would  have  to  be  automatically  adjusted  over  a  large 
range,  and  the  adjustment  would  have  to  be  very  precise.  The  lengths 
wo\?ld  also  have  to  be  adjusted  for  pressure  changes,  since  the  entrain¬ 
ment  of  flow  and  the  velocity  of  th<;  jet  significantly  affect  the  frequency 
of  oscillation. 

Another  oscillator  is  one  in  which  a  linear  mechanical  spring -mass 
system  (essentially  a  tuning  fork)  ccntrols  the  frequency.  Since  the  damping 
is  small,  the  frequency  is  limited  tc  the  vicinity  of  the  natural  frequency 
of  the.  spring-mass  system;  compensation  for  changes  in  pressure  and 
temperature  are  smaller  and  need  not  be  so  precise.  For  this  reason., 
only  oscillators  controlled  by  spring-mass  systems  have  been  considered, 
although  an  oscillator  with  a  pneumatic  resonator  might  be  satisfactory 
for  seme  purposes. 

To  meet  the  requirement  of  pressure  insensitivity,  the  goal  has 
been  driving  pressure  forces  which  are  as  nearly  as  possible  in  phase 
with  the  velocity  of  the  mass  so  that  they  subtract  from  the  damping  forces 
but  neither  add  to  nor  subtract  from  the  spring  forces.  Then  the  phasing  of 
the  pressure  forces  would  be  independent  of  both  temperature  and  pressure 
except  to  the  extent  that  the  modulus  of  elasticity  of  the  spring  and  the 
dimensions  of  the  vibrating  system  vary  with  temperature.  The  amplitude 
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Receivers 


Fig.  71.  Jet  Amplifier-Oscillator, 
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Fig.  72.  Vibrating  Shaft  and  Disk. 


Cutoff  region 

Fig.  73.  A  Suggested  Driving -Force  Versus  Radius  Curve. 
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of  motion  would  change  with  pressure,  however,  requiring  some  non- 
critical  control  to  increase  the  range  of  pressures  over  which  the  device 
would  operate  usefully. 

6.2.  Circular -Motion  Oscillator 

A  disk  on  a  shaft  of  negligible  mass  can  be  made  to  vibrate  (without 
rotation)  in  a  circular  path.  This  motion  can  be  considered  as  the  super¬ 
position  of  two  mutually  perpendicular  lateral  vibrations  with  equal  ampli¬ 
tudes  and  a  ninety  degree  phase  difference.  See  Fig.  72. 

From  steady-state  circular  motion,  the  conditions  for  equilibrium 

are 

mrco^  =  kr  4-  G  (169) 

in  radial  direction 


Brco  =  F  (170) 

in  the  tangential  direction  where 
m  =  mass  of  the  disk 
co  =  the  angular  frequency  of  vibration 
r  =  the  displacement  from  the  center  position 
k  =  the  lateral  stiffness  of  the  shaft 
G  and  F  =  forcing  functions 

Notice  that  if  G  -  0,  then  =~'jk /m,  F  might  be  a  function  of  r,  F(r). 

Thus,  if  a  configuration  can  be  found  which  will  provide  a  driving  force 
which  is  perpendicu1  ar  to  the  displacement  at  all  points  of  the  cycle  and 
of  constant  magnitude  for  any  given  amplitude,  then  the  frequency  of 
oscillation  will  be  the  undamped  natural  frequency  of  the  disk-shaft  system. 
F(r)  might  be  a  function  of  pressure  and  temperature,  but  only  the  amplitude 
of  oscillation  is  affected  by  changes  in  pressure  and  temperature. 

The  relationship  F(r)  might  have  a  sharp  cutoff  at  a  certain  amplitude, 
r,  as  shown  in  Fig.  73.  The  equilibrium  amplitude  is  determined  by  the 
intersection  of  the  damping  force  and  the  F(r)  curves. 


6.3.  Problems  Affecting  the  Configuration 

Two  problems  greatly  affect  the  configuration  and  limit  the  usefulness 
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and  pressure  forces  which  affect  the  effective  spring  constant.  This  will  be 
discussed  in  detail  iater. 


The  second  problem  is  that  oi  the  stability  of  the  circular -motion 
equilibrium.  It,  may  be  desirable  to  substitute  for  the  ideal  F(r)  another 
F(r,G),  where  0  is  the  angle  or  ^  uxlt.  for  example,  in  the  shaded  portion 
of  Fig.  74a,  F{r,  0)  might  be  zero.  Now  the  intersection  of  the  damping 
force  and  the  driving  force  curves  has  no  meaning  except  in  a  qualitative 
way,  and  equilibrium  motion  is  not  strictly  harmonic,  although  perhaps 
very  nearly  so. 

Assuming  harmonic  motion,  it  can  be  seen  that  the  circular  equili¬ 
brium  motion  might  be  unstable.  A  sl;ght  increase  in  the  amplitude  in  the 
X-direction  would,  because  of  the  slope  of  the  F(r,  0;  curve,  result  in  a 
decrease  in  the  Y  direction  driving  force  and,  consequently,  the  Y-direction 
amplitude.  By  similar  argument  the  X-direction  amplitude  would  increase 
some  more,  and  the  final  steady-state  motion  would  be  elliptical. 

For  certain  applications,  an  elliptical  motion  would  be  acceptable, 
but  more  often,  circular  motion  is  desired.  If  the  number  of  regions  of 
positive  driving  force  (shaded  areas  in  Fig.  74b)  is  increased,  then  F(r,Q) 
approaches  .-a  ideal,  co  uinuou3  F(r).  For  this  case,  a  more  careful  anal /sis 
is  required  to  determine  the  stability  of  the  circular  motion.  Such  an  analysis 
was  motivated  by  initial  experimental  observations  of  sharply  elliptical  rather 
than  the  desired  circular  motion. 


6.4.  Stability  Analysis  for  Continuous  F(r). 

In  this  analysis  the  motion  of  the  end  support  is  taken  into  account. 

It  is  assumed  that  the  movable  end  support  of  Fig.  75  moves  up  and  down 
without  sidewise  or  rotational  motion  and  that  its  motion  is  directly  propor¬ 
tional  to  the  motion  of  the  vibrating  disk.  (We  assume  that  the  wire  is  stiff 
axially).  The  result  is  that  ^he  effective  mass,  spring,  and  damping  constants 
for  vibration  in  the  radial  direction  are  different,  generally,  from  those  for 
vibration  in  the  tangential  direction.  (For  purely  circular  motion  of  the  disk, 
:he  end  support  does  not  move  at  all). 
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Fig. 74.  Stability  of  an  Oscillator  with  a  Discontinuous  Driving  Fore 


Movable  end  support 
Spring  wire 


L. 


Vibrating  disk 


7777777 


Fig.  75.  General  OsciUator  Configuration  for  the  Stability  Analysis, 


The  kinematic  equations  are: 


Displacement 


r  u 

r 


\i  t :  ) 


Velocity  =  ru  r  rOu„ 
7  r  0 


(172) 


Acc  deration 


fv  -  rtT  i 

U  T 

«  «  • 

2r0  +  r  0 

1 _  _J 

r 

_  _ 

where  is  a  unit  radial  vector,  u^  is  a  unit  tangential  vector,  and  r 
and  0  are  polar  coordinates  in  the  system  whose  center  is  the  neutral 
position  of  the  disk. 


(173) 


Next,  using  only  the  radial  acceleration,  velocity  and  displacement 
and  assuming  linear  damping,  the  dynamic  equation  for  the  radial  direction 
is  written. 


m 

r 


+  B  f  +  K  r 
r  r 


0 


(174) 


where  r-subscripts  refer  to  the  radial  direction. 
Let 

-y  K  B 

2  =  r  .  ~r  ~  r 

wr  "  - -  ’  2^rWr - : 

m  m 

r  r 


(175) 


then 


fr‘  -  rO2)  -f  2C,  o  r  4-  co  2  r  =  0 
r  r  r 


0  76) 


The  dynamic  equation  for  the  tangential  direction  is 


m 


o  !_ 


rO  +  2r0 


+  Bft  rO  -  F(r)  =  0 


(177) 


where  ©-subscripts  refer  to  the  tangential  direction. 
Define 


2  _  K0  „ 

“0  = - 2hl“8 

m0 


m. 


,  and  f  (r)  = 


i  _  F(r) 


m. 


078) 
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where  F(r)  is  the  continuous  driving  force; 
then 

»  »  ♦  « 

\r©  +  2r0)  f  2£a  r0  .  f  (r)  ^  0-  (1 79) 

In  analyzing  the  stability  of  circular  motion,  we  can  linearize 
equations  (176)  and  (177)  and  investigate  the  effect  of  a  small  departxire 
from  equilibrium.  Equation  (176)  becomes 

aY  -  2rOAG-Q^Ar+2£,  u>^Ar4-  A  r  =  0  (180) 

Now 

Af(r)  =  —  Ar  =  f'Ar,  (181) 

dr 

so  Eq.  (179)  becomes 

*  •  •  •  «  •  *  * 

0  Ar  +  rAC  +  ZhAO  +  2©Al  +  2£qU)q©  Ar  +  2 w^rA©  -  f'Ar 

For  the  equilibrium  condition  under  investigation, 

r  is  constant 
V  =  r  =  0 

0  =  (183) 

•  • 

0=0 

so  that  Fqs.  (18^)  and  (182)  become 

AY  +  21,  uir  2  -«n2)  Ar  -  2r  <*>  Aft  =  0  (184) 

r  r  r  ©  0 

•  •  • 

r  A  0  +  2  <o^a  r  +  2£q Ar  +  2^  ^  rAO  -  f'Ar  =  0  (185) 


=  0 

(182) 
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When  the  Eqs.  (184)  and  (185)  are  combined  and  the  aQ  i.erms  eliminated 
(one  differentiation  involved),  the  result  is 


r 

aV‘+  2«  (t  +  CQ)  A*r  b  |u>  Z{1 
r  r  0  j _ r 


* 3  ,ji o 


n  r 


H#"r  (»/  +  “02>  -  2  “#f 


A  r  =  0 


(186) 


The  stability  of  several  situations  might  be  investigated.  The 
simplest  one  is  that  for  which 


“r  =  “0  *  ^0 


(187) 


(end  supports  rigid,  supporting  wire  stretchable)  and 


iri6<'  i  (188> 

(a  very  reasonable  assumption  in  any  case  of  interest  here).  We  now  get 

(D3  +  I;2  .«•  4ooq2D  !  4^  ug)  Ar  =  2  u>gf'Ar  (189) 

where 

D  =  —  (190) 

dr 

Figure  This  a  polar  plot  of  the  left-hand  side  of  Eq.  (189)  with  D  replaced 
by  ju>.  The  Nyquist  stability  criterion  for  this  inverse  Nyquist  plot  re¬ 
quires,  for  a  stable  circular  motion,  that  2  have  a  value  of  the  real 
axis  within  the  curve;  for  stability, 

-  6;0wc2<  f'  <2  ;u>02  (191) 

This  result  is  expressed  graphically  in  Fig,  77.  The  lower  limit  for  the 

2 

cutoff  slope,  -6£q  Wq  ,  is  the  one  of  interest.  Physically,  an  instability 
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Fig.  77.  Interpretation  of  Stability  Analysis. 

I 


Fig.  78.  Inverse  Nyquist  Plot  for  System  With  Lag. 
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at  this  slope  is  not  very  spectacular;  the  path  of  motion  inerely  begins  to 
become  elliptical.  Note,  from  the  polar  plot,  that  the  frequency  for 
oscillations  of  Ar  about  the  equilibrium  radius  is 

The  range  of  permissible  cutoff  slopes  can  be  increased  by  intro¬ 
ducing  a  lag  into  the  driving  function.  At  the  operating  radius,  we  use 


Af  _  g1 

Ar  (TD  +  1 ) 


(192) 


to  replace  the  f'  previously  used.  Physically,  this  means  that  changes  in 
the  driving  force  are  related  to  changes  in  the  radius  at  an  earlier  time, 
i.e.,  at  an  earlier  angle,  0,  in  the  motion.  If  the  lag  angle  is  almost 
ninety  degrees,  then  an  increase  on  the  amplitude  in,  say,  the  X-direction 
will  result  in  a  decrease  in  the  driving  force  in  the  same  (X)  direction, 
but  ninety  degrees  later. 

To  find  the  effect  of  this  lag  upon  the  stability,  the  left-hand  side 
of  Eq.  (189)  and  the  polar  plot  of  Fig.  74  are  multiplied  by  (jToj  +  1). 

If  Tu  »  1  at  u  =  coq,  the  w  =  point  will  be  rotated  count-  rclockwise 
almost  ninety  degrees,  and  the  magnitude  at  co  =  will  be  multiplied  by 
T  w  q.  The  new  polar  plot  is  shown  in  Fig.  78.  We  find  that,  for  stability. 


and 


Of  course. 
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(193) 


(194) 


(195) 


(196) 
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The  optimum  lag  is  actually  slightly  less  fhan  ninety  degrees  and  rotates 
the  point  of  maximum  amplitude  at  u,  -  1.15  to  the  negative  real  axis. 

If  the  g '  cutoff  slope  is  too  steep,  the  resulting  ''instability’'  has  an 
oscillation  frequency  o  t  This  means  that  the  motion  is  not  elliptical 

but  epicyclic. 

Although  the  use  of  the  lag  provides  for  a  much  steeper  allowable 
cutoff  slope  and  consequently  a  greater  pressure  range  for  satisfactory 
operation,  its  use  increases  the  complexity  of  the  oscillator.  For  this 
reason,  the  oscillators  built  and  tested  all  lacked  the  lag. 

A  lead  could  also  be  used  and  would  be  even  more  successful  than 
a  lag,  although  it  would  make  a  still  more  complex  system.  In  this  case 

—  =  h’  (TD  +  1}  (197) 

A  r 

would  be  used  to  replace  £;  in  Eq.  (  39).  The  polar  plot  would  be  divided 
by  (jTu>  +  1).  The  result  (See  Fig.  79)  would  be  a  stable  system  for  all 
cutoff  slopes  and  all  pressures. 

A  physical  interpretation  cf  these  results  is  now  given  with  the  aid 
of  Fig.  80,  which  shows  a  hypothetical  disk  motion  in  which  r^  is  t  large. 
The  simplest  system  senses  the  error  in  r&  and  corrects  r^  by  correcting 
Fa«  If  is  reduced,  rc  in  turn  becomes  elongated,  rc  reduced,  ax<d  the 
motion  becomes  elliptical.  The  system  with  a  ninety  degree  lag  corrects 
r  „  If  r  is  reduced,  r  is  consequently  increased  more,  and  the  motion 

C  C  d 

becomes  eccentric.  In  the  system  with  a  leac,  r  itself,  rathe  r  than  Its 

a 

innocent  complex  ents,  is  corrected. 

Now  suppose  that  f  :o^,  but  that  ^  L ,  .  Assume  a  system 

like  that  of  Fig.  81.  If  the  amplitude  of  radial  vibration  of  the  disk  is  6, 

then  the  amplitude  of  vertical  vibration  of  the  top  support  is  approximately 

2 

6  /L.  The  effective  mass  of  the  support  as  seen  at  the  disk  is 


'ig.  81.  Notation  for  Analysis  of  Oscillator  with  Movable  End 


inc  the  natural  radial  frequency  is 


U> 


r 


f 


M 


(199) 


.t  the  disk  is  vibrating  in  a  circular  path,  the  top  support  will  not  move 
it  all,  so  the  tangential  natural  frequency  is 


(200) 

Fouth's  stability  '*  pterion  applied  to  Eq,  (175)  yields,  for  the  lower 
limit  of  stability, 


to  = 
0 


T 


K 


rr. 


to  CO 

~  t,  (CO  2  X  to  2) - £_ 

to  0  r  0 

0  0 


(Cr  +  C0) 


CO  2  n  4-  4C  CA)  x  3  to  2 
j^r  ^r  0  0 


(201) 


Assuming  «1,  this  becomes 


£'  =  —  <-irurZ  -  Hr  "82  -  259  “62)  (202) 

W  - 

9  $ 

A 

Now,  with  hindsight  not  yet  available  to  the  reader,  we  assume  that 
'/L  =  0,  02  and  M  =  20  m.  sis  a  result,  =  co^/1.4. 

An  application  of  Routh’s  criterion,  Eq.  (202),  show's  that  for 
lability,  the  steepest  allowable  cutoff  slope  is  -3.944gto^  compared  to 
6£gtOg2  for  rigid  end  supports  and  stretchable  spring  wires,  or,  equi- 
alently,  zero  support  mass.  Th3.s  indicates  that  it  is  necessary  to 
linimize  the  movable  support  mass. 

Extra  damning  might  be  given  the  vertical  motion  of  the  support. 

2 

:  we  assume  that  =  10  Eq.  (202)  shows  that  f’m^n  =  -  26.  5  u>  q  . 
he  additional  damping  seems  to  be  desirable. 
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6,  5.  Temperature  Compensation 

The  frequency  of  the  ideal  circular-motion  oscillator  is  not 
affected  b”  changes  in  the  physical  properties  of  the  gas  flowing  through 
it.  Changes  in  temperature  will,  however,  change  in  the  size  and  the 
stiffness  of  the  mechanical  parts. 

Changes  in  stiffness  dominate  changes  in  size.  According  to 
simple  beam  theory,  the  lateral  stiffness,  K  of  a  shaft  is  prop  >rtional 
to  El/L  where  I  is  the  second  area  moment  of  the  cross-section,  E 
is  the  modulus  of  elasticity  in  tension,  and  L  is  the  length.  1  is  pro- 

4 

portional  to  L  at  all  temperatures,  so  K  is  proportional  to  El.  Figure 
82  (Ref.  43)  shows  the  variation  in  E  with  temperature  for  some  steels. 
Along  the  linear  portion  of  the  curves, 


2.  13  x  10 


(203) 


EAT 


From  any  materials  handbook,  a  typical  value  for  the  coefficient  of  linear 
expansion  of  steel  is 


6.7  x  10 


(204) 


LAT 


For  purposes  of  analysis,  the  change  in  size  will  be  neglected.  Fine  ad¬ 
justment  of  an  actual  oscillator  can  overcome  this  neglect. 

Now  we  have 


2x10 


(205) 


KAT 


eat 


If  AT  -  1000°F,  then  AK/K  ~  20  per  cent.  This  number  will  be  kept  in 
mind  henceforth.  Notice  that  linearity  is  here  assumed. 

Of  several  possible  mechanisms  investigated  for  temperature 
compensation,  the  only  one  that  seems  really  satisfactory  for  such  a 
large-scale  compensation  is  the  varying  of  the  axial  compression  or 
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tension  of  the  supporting  wires  by  means  of  a  sealed  bellows.  VTth 
this  mechanism,  one  end  support  can  be  made  non-rigid,  as  required 
for  linearity  of  the  vibrating  spring,  and  the  actual  compensating  motion 
can  be  kept  very  small,  although  the  force  changes  by  a  substantial 
amount. 

Figure  83  shows  schematically  the  oscillator  with  the  bellows 
temperature  compensator.  The  bellows  are  designed  to  be  flexible  in 
the  axial  direction  and  stiff  in  the  lateral  direction.  The  vibrating  disk 
is  supported  for  by  four  wire  springs.  Because  of  the  parallelogram 
action  of  the  four  wires,  the  net  moment  felt  by  the  movable  end  support 
is  small,  and  the  only  force  to  be  resisted  is  the  lateral  force. 

Sealed  chambers  A  and  B  encircle  the  axis.  The  gas  in  chamber  B 
is  at  a  higher  pressure  than  in  chamber  A.  The  pressure  difference  applied 
to  the  movable  support  is  proportional  to  the  absolute  temperature  of  the 
gas  in  chambers  A  and  B.  With  the  spring  of  the  bellows  or  with  a  separate 
spring,  the  wires  can  be  given  an  initial  compression.  As  temperature  in¬ 
creases,  the  stress  in  the  wires  changes  from  compressive  to  zero  to 
tensile.  Since  the  fatigue  resistance  of  the  metal  becomes  less  at  high 
temperature,  the  adjustment  should  be  such  that  the  compressive  stresses 
become  zero  at  the  highest  temperature. 

The  use  of  two  chambers,  A  and  B,  eliminates  the  influence  of 
changes  in  the  environmental  pressure.  In  some  applications,  one  chamber 
might  be  omitted.  If  the  oscillator  is  subjected  to  acceleration,  the  counter¬ 
balancing  mechanism  shown  in  Fig.  81  Is  necessary  in  order  that  the  os¬ 
cillator  be  insensitive  to  acceleration  in  the  axial  direction.  The  oscillator 
will  still  be  sensitive  to  accelerations  and  vibrations  in  the  lateral  direction 
to  a  degree  depending  upon  its  mass  ana  spring  constants.  Because  of  the 
sensitivity  to  forces  in  the  axial  direction,  the  oscillator  driving  ar.d  porting 
mechanism  at  the  vibrating  disk  must  be  designed  so  that  the  axial  pressure 
forces,  as  well  as  the  radial  flow  forces,  are  kept  very  small. 

6.  6.  Lateral  Stiffness  of  the  Springs  as  a  Function  of  the  Compression 

For  the  cantilever  beam  shown  in  Fig.  84  the  maximum  deflection  is 
given  by  (Ref.  44). 
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.  p*  t  , 
tan  i -  -  i 


P1  I  P 


This  can  be  modified  for  the  case  of  Fig.  85: 


6=  M  LSL  tan,  121 
P'J  P'  4  El 


Equation  (207)  can  be  put  into  the  form  (Ref.  45) 


(206) 


(207) 


4  El 


tan  a  -  a 


\  a' 


(208) 


where 


P1  L 
4  El 


(209) 


For  P1  =  0,  (K^at  )p,_  q  is  12EI/L  from  simple  beam  theory  or  by  reducing 
of  Eq.  (208).  Then 


^fatj  p= 


tan  a  -  a 


Den  Hartog  then  expands  this  in  an  infinite  series: 


(210) 
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2  \2 


(211) 


=  1-0.  987 


.  0.012  ^  -  0. 
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2  2 

where  P'cr^  =77*  EI/L  either  from  the  Euler  buckling  formula  or  the 

reduction  of  Eq.  (210)  for  K„  =  0. 

f  at 

^  ^fat^^fat)  P=  0  *s  P^0^e<^  as  a  function  of  P*/P'  ^  the  function 

is  nearly  a  straight  line.  Assuming  that  it  is  linear,  and  assuming  a  linear 
relation  between  change  in  temperature  and  change  in  the  modulus  of 
elasticity. 


K  =  Kq(  1  -  aAT) 


P 


P  ..  (1  -  qAT) 
o  cnt  ' 


(212) 


where  P  ,  is 
o  crit 


the  total  buckling  load  for  all  the  wires. 


K  and  P  are  for  T  =  0, 

o  o  crit  * 


AT  =  temperature  change 


(213) 
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o.  =  - 1 

Eat 


and 


o  +  AT  \ 

- - J  +  CqA  f  (1  -  aA  T)  (214) 

o  • 

where  Pq  is  the  pressure  difference  between  chambers  A  and  B  in  Fig.  81  and 
Co  is  the  bellows  spring  constant,  both  for  AT  =  0,  and  Af  is  the  initial 
deflection  of  the  bellows.  A*  is  the  area  of  the  ends  of  the  chambers.  Then, 
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For  temperature  insensitivity,  the  coefficient  of  AT  must  be  zero; 


(215) 
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So  then 


/ 

/ 

C  A  i 
o 

P  crit 

V 0 


T  P 
o 


o  crit 


AP  A 
o 

Po  crit 


c  a  i 

o 

Po  crit 


(217) 


where  K  is  the  lateral  stiffness  of  the  springs  supporting  the  disk.  In  the 
physical  system  APq  and  A  l  are  both  easily  adjusted.  Hence  both  K  (and 
the  frequency)  and  the  temperature  compensation  can  be  adjusted. 

The  two  assumptions  of  linearity  made  here,  namely  a  =  a  constant 
and  =  (^at)p_0  “  ^/^Critjcannot  f°r  really  precise 

vompensation  without  a  much  more  extensive  investigation  of  the  relation 

а,  between  modulus  of  elasticity  and  temperature.  Assume  a  actually  is 
constant.  For  a  1  000°F  temperature  range,  the  maximum  P  in  Eq.  (212) 

is  20  per  cent  of  P  The  deviation  of  Eq.  (212)  from  a  straight  line  over 

this  range  is  about  +_  0.006  per  cent.  Remember,  however,  that  a  beam  in 
bending  is  a  linear  spring  only  for  small  deflections.  The  size  of  the  assumed 
linear  range  is  related  to  the  precision  required. 

б. 7.  Configurations  to  Provide  Driving  Force  and  Output 

The  first  circular-motion  oscillator  designed  and  built  consisted  of 
a  flat  disk  suspended  by  four  wires  as  in  Fig.  86.  The  circular  path  of  the 
disk  was  in  a  plane  adjacent  to  a  plate  containing  air  ports  which  were  covered 
and  uncovered  by  the  motion.  Pressure  was  supplied  to  the  housing  for  the 
system.  Output  lines  were  connected  to  the  ports  in  the  base  plate.  The  two 
crescent -shaped  ports,  shown  in  Fig.  86b  and  c,  provided  two  outputs  180 
degrees  out  of  phase  with  other.  (A  photograph  is  shown  on  Fig.  87a). 

The  driving  force  for  the  disk  was  provided  by  vanes  set  at  a  spiral 
angle  around  its  periphery.  The  flow  through  the  vanes  was  greatest  at  the 
largest  opening  uncovered  by  the  disk,  and  the  flow  force  caused  by  the 
vanes  was  always  tangential  to  the  direction  of  the  displacement  of  the  disk. 
See  Fig.  86c. 

When  tested,  this  oscillator  showed  some  sensitivity  to  pressure 
because  of  two  undesirable  forces  affecting  the  spring  constant  One  of 
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these  was  the  column  loading  on  the  face  of  the  disk.  The  other  radial 
pressure  force  was  a  result  of  flow  into  the  partially  closed  ports. 

The  flow  force  on  the  edge  of  a  plate  closing  a  hole  was  measured 
roughly  and  found  to  be  generally  as  shown  in  Fig.  87.  Because  of  the 
skew  of  the  curve,  this  radial  force  could  not  be  balanced  adequately  by 
the  corresponding  force  on  the  opposite  edge  of  the  disk. 

Another  design  was  built  and  tested.  The  vibrating  disk  had 
eight  holes  evenly  spaced  in  a  circle.  Each  hole  was  perpendicular  to 
the  radius  at  that  point  but  slanted  with  respect  to  the  axis  of  the  disk. 

Thus,  flow  through  a  hole  was  giver  a  velocity  a.i  d  momentum  in  the 
tangential  direction.  The  base  plate  contained  another  ring  of  eight  holes 
at  a  slightly  greater  radiu3.  See  Figs.  88and89b.  A  displacement  of  the  disk 
from  the  center  position  resulted  in  a  greater  opening  f  >r  flow  on  one 
side  than  the  other,  and  therefore,  a  tangential  flow  force.  As  the 
amplitude  became  greater,  however,  the  disk  and  plate  holes  overlapped, 
less  so  that  the  flow  force  was  reduced.  The  amplitude  was  not  limited, 
though.  After  the  circular  motion  attained  a  certain  amplitude,  further 
increases  in  pressure  caused  the  motion  to  continue  to  grow  with  an 
elliptical  shape.  The  stability  analysis  was  then  made  to  understand  this 
phenomenon. 

It  seems  reasonable  to  expect  that  the  flow  forces  arising  from  flow 
through  a  hole  be  less  than  the  T.OW  forces  on  an  open  edge,  since  the  re¬ 
duction  of  pressure  by  Bernoulli  flow  is  felt  at  least  partially  all  around 
the  interior  of  the  hole.  The  first  of  the  two  oscillators  just  described 
has  a  pressure  sensitivity  of  about  1  per  cent  per  psi  supply  pressure.  The 
second  had  a  sensitivity  of  about  0.  02  per  cent  per  psi.  The  comparison  is 
crude  because  the  oscillators  were  only  roughly  of  the  same  size. 

Both  oscillators  had  the  vibrating  disk  upstream  so  that  the  output 
could  be  taken  directly  from  the  base  plate.  A  simple  experiment,  the 
results  of  which  are  summarized  in  Fig.  89,  emphasized  that  it  would 
have  been  better  to  place  the  disks  downstream. 

An  oscillator  with  a  pneumatic  lag  in  the  amplitude -sensing  mechanism 
was  considered  but  not  built,  even  though  the  stability  analysis  had  indicated 
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Fig.  88.  Second  Test  Oscillator. 


that  a  lag  was  quite  desirable.  A  possible  configuration  is  given  in  Fig.  90. 

The  holes  in  the  upper  disk  whic.»  meter  the  supply  air  are  aligned  with  the 
corresponding  holes  in  the  supply  plate  when  the  disk  is  centered,  but  be¬ 
come  more  and  more  misaligned  (orbiting  around  the  supply  hole)  as  the 
amplitude  of  oscillation  increases.  The  chamber  causes  a  lag  between 
changes  in  flow  through  the  upper  disk  and  the  driving  force  provided  by  the 
vanes  on  the  lower  disk.  The  flow  is  in  the  proper  direction  for  keeping  the 
radial  flow  forces  low,  and  consequently,  the  oscillator  as  shown  has  no  useful 
output;  auxiliary  ports  are  necessary. 

The  third  oscillator  built  was  quite  the  same  in  principle  as  the 
second,  but  the  flow  was  reversed.  The  disk  was  made  very  thin  so  that 
the  sides  of  the  holes  would  not  interfere  with  the  flow  (see  Fig.  91a)  and 
vanes  were  placed  immediately  behind  the  holes  to  provide  the  driving 
force,  as  shown  in  Fig.  91b.  Instead  of  being  in  a  slightly  larger  ring, 
the  holes  in  the  supply  plate  were  drilled  in  line  with  the  holes  in  the  disk 
and  then  half -plugged,  as  shown. 


It  is  interesting  to  examine  in  greater  detail  the  operation  of  this 
oscillator  and  to  examine  especially  the  action  of  the  amplitude  cutoff. 
Consider  air  flow  in  the  configuration  of  Fig.  91b.  The  momentum  of  the 
fluid  flowing  through  the  holes  is  WV,  where  W  is  the  mass  flow  rate  and 
V  is  the  velocity.  If  the  vane  is  set  at  a  45  degiee  angle,  then  the  lateral 
force  on  the  vane  is  approximately  WV/lfzT  This  assumes  that  the  flow 
is  straight  through  the  holes,  which  is  not  strictly  true.  In  fact,  with 
this  assumption  we  overlook  a  small  source  of  radial  force  which  will  be 
discussed  later.  It  is  assumed  also  that  the  velocity  of  the  disk  is  small 
compared  to  the  velocity  of  the  air.  It  is,  in  fact,  30  inches/second  for 


u  j  =  500  cps  at  an  amplitude  of  0.  01  inch. 

The  product  WV/|  Z  is  now  derived  for  both  incompressible  and 
compressible  flow.  For  incompressible  flow. 


(218) 


where  AP  is  the  pressure  drop  across  the  holes  and  p  is  the  fluid  density. 
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Fig.  90.  Oscillator  with  a  Pneumatic  Lag  in  the  Amplitude  -Sensing 
Mechanism. 


Fig.  91.  Third  Test  Oscillator  Driving  Mechanism. 


W  =  P  VA 


(219) 


and 

WV  _  PV2A  AAP 

F=  fz  =  fr  fr 

For  isenttopic  compressible  flow, 


V  = 


(220) 


(221) 


where  2vl  is  the  Mach  number,  and  •y  is  the  ratio  of  specific  heats. 

WV  =  pV2A  =  pAM2yRT 
Jsing  the  equation  of  state 


(222) 


we  get 


P 

p  =  - 

RT 


_  WV  PAMZv 

~ir~  fr 


(223) 


(224) 


where  P  is  the  stagnation  (supply)  pressure,  and  P  is  the  pressure  where 

°  2 
V  is  taken.  From  isentropic  flow  tables(Ref.  46)  P/PqM  y  can  be  calculated, 

and  it  is  plotted  as  a  function  of  P/P  in  Fig.  92.  For  sonic  flow  (P/P  <  0.  528), 

O  O 

P/P  M  y  remains  constant  at  0.7  39. 


In  each  of  these  cases  the  driving  force,  F,  is  proportional  to  an 
effective  area,  A.  Figure  93  represents  the  disk  given  an  X-direction  dis¬ 
placement  from  its  neutral  position  with  respect  to  the  supply  plate.  On 
each  vane  the  momentum  force,  which  is  proportional  to  the  common  area 
of  the  supply  plate  hole  and  the  disk  hole  at  that  location,  is  indicated  by  an 
adjacent  vector  in  the  figure.  Ideally  the  net  X-direction  force  on  the  disk 
is  zero.  This  is  not  strictly  true  because  the  flow  is  given  a  slight  deflection 
in  going  through  the  holes  and  impinges  upon  each  vane  with  a  slightly  different 
initial  angle.  Setting  each  vane  at  a  slight  angle  with  the  radius  might  correct 
this. 
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Fig,  92.  P/P  Versus  P/P 
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Fig.  93.  Third  Oscillator  Driving  Forces  tor  a  Given  Displacement. 
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Fig.  94.  Effective  Area  Versus  Radius. 


The  net  Y-direction  force  is  the  desired  tangential  driving  force. 

This  net  force  is  proportional  to  the  effective  area,  A,  which  is  the  sum 
of  all  the  open  areas  multiplied  by  the  proper  angle  cosines.  This  effective 
area  was  obtained  graphically  as  a  function  of  the  displacement  and  is 
presented  as  Fig.  94. 

In  the  stability  analysis  it  was  found  that  for  =  oj^/1.4  the  steepest 
permissible  cutoff  slope  is  about  twice  the  damping -force  slope.  In  Fig.  93 
a  straight  line  has  been  placed  so  that,  at  its  intersection  with  the  area  curve, 
this  relation  between  slopes  is  met.  The  area  curve  is  proportional  to  the 
driving  force,  so  the  straight  line  is  proportional  to  the  damping  force. 

This  intersection  gives  tne  maximum  amplitude  of  oscillation  with  circular 
motion  as  0.66  times  the  diameter  of  the  holes  in  the  disk.  To  find  the 
pressure  range  for  the  oscillator,  a  minimum  useful  amplitude  must  be 
chosen.  The  ratio  of  the  maximum  to  the  minimum  driving  force  is  deter¬ 
mined  by  the  multiplier  needed  to  make  the  damping  and  area  curves  inter¬ 
sect  at  the  minimum  amplitude.  If  the  minimum  amplitude  is  0.  35  times 
the  diameter  of  the  disk  holes,  then  the  ratio  of  maximum  to  minimum 
driving  force  is  found  to  be  2.  35.  If  the  driving  force  is  proportional  to 
the  pressure  (it  is  for  critical  flow  and  for  incompressible  flow),  then  the 
ratio  of  maximum  to  minimum  pressure  is  also  2.  35. 

Different  types  of  damping  curve  are  shown  in  Fig.  95.  The  square- 
law  damping  is  somewhat  more  realistic  than  the  linear  damping.  The  addi¬ 
tion  of  Coulomb  damping  introduces  a  stable  200 -amplitude  equilibrium  point, 
which  can  keep  the  oscillator  from  starting  itself. 

Equating  the  driving  force  to  the  damping  force,  we  get 


P*A  2 

■s2^e  rm 


YT 


(225) 


where  P'  is  a  pressure  function  and  m  is  the  mass  of  the  oscillating  disk. 

Since  A  is  a  function  of  the  radius,  the  preceding  equation  can  be 
rewritten 

P’d 


r 

“H" 


(226) 


2;q  u>0  m 


-  194  - 


0.5  1.0 


r/d 

Fig.  96.  (r/d)/(A/d“)  and  (A/d^)/{r/d)  Versus  r/d. 
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where  (r/dJ/A/d-),  a  function  of  r/d  alone,  has  been  plotted  in  Fig.  96.  The 
term  on  the  left  in  Eq.  (226)  is  independent  of  the  size  of  the  oscillator  for 
geometrically  similar  oscillators.  To  show  this,  we  write 


d 


This  incidentally,  is  itself  of  interest,  although  not  surprising. 
Then 


d 


2^e  "o '  m 


CL 


The  damping  force  is 


2  2 

2;0  om  a  d 


The  air  damping  force  is  proportional  to  the  projected  area  of  the  moving 

2 

parts  and  to  their  velocity,  or  for  geometrically  similar  oscillators,  d  . 


Therefore, 


CQ  *  1 


The  internal  damping  moment  is  given  by  the  formula 


Mdc 


=  cr 


(227) 


where 


cr  ,  a  c  cr  6u>_  a  1 
a  w 


Hence  the  internal  damping  force,  proportional  to  M^/L,  is  proportional 
to 


T  7 

J-  j  *- 

-  a  a 


cL 


and  again 


cl 
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Equation  (226)  would  be  used  in  determining  the  size  of  the  oscillator  for 

а.  specified  frequency  and  supply  pressure  level.  Some  information  about 

would  be  required. 

б. 8.  Design  of  the  Springs 

The  next  task  is  the  design  of  the  spring  support  system  for  the 
disk.  A  minimal  size  is  desired.  For  each  wire,  considered  as  a  beam, 
the  bending  formula  is 

QL3 

6  =  -id—  (228) 

12EI 

and  the  spring  constant  is 

J9l  =  — 2£L  (229) 

6  L3 

For  the  disk-spring  system,  the  spring  constant  is  affected  by  n,  the 
number  of  wires: 


K  = 


n  UEI 


(230) 


where  K  is  the  spring  constant,  n  is  the  number  of  half-wires,  L  is  the 
wire  length  on  each  side  of  the  disk.  For  a  round  wire,  the  moment  of 
inertia,  I,  islTc^/4,  where  c  is  the  wire  radius.  The  maximum  moment 
on  a  wire  is  the  product  of  the  spring  constant,  the  deflection,  and  the 
length: 

M  =  K6L  (231) 

max 


The  maximum  stress  is  given  by  the  relation 


M  c 
n'.arc 

I 


ab 


max 


(232) 
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Equations  (230),  (231),  and  (232)  combined  yield  the  two  equations. 


.  2 

L>  cr, 
o 


12  E6 


and 


Equation  (234)  rearranged  gives 


In  operation. 


K 

3En7 T 


so  that 


K  =  m  Wg 


2 


(233) 


(234) 


(235) 


(236) 


(237) 


Equation  (235)  gives  the  minimum  length  for  a  specified  disk  configuration, 
natural  frequency,  and  allowable  stress. 

6.  9.  Construction  of  the  Third  Oscillator 

The  output  of  the  third  oscillator  was  a  flow  of  air  modulated  by 
the  vibration  of  the  disk.  The  output  air  and  the  driving  force  air  were 
supplied  separately.  The  modulation  of  the  output  flow  was  done  by  four 
ports  in  the  disk  arranged  in  a  ring  at  a  larger  radius  than  the  driving  holes 
and  vanes.  There  were  corresponding  ports  in  the  adjacent  stationary  port 
plates,  and  the  result  was  a  set  of  four  valves  similar  to  conventional  plates 
valves.  In  making  four  valves  it  was  intended  that  different  combinations 
of  the  outputs  could  give  a  single-phase  output  with  frequency  u>  q,  a  push- 
pull  output  with  a  frequency  a  push-pull  output  with  frequency  2 
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and  a  s ingle  output  with  frequency  Details  of  the  porting  are  shown 

in  Fig.  97.  The  vibrating  disk  is  detailed  in  Fig  98;  Figs.  99  and  99a 
show  the  whole  oscillator.  Detail  drawings  of  all  the  parts  are  in  the 
Appendix. 

The  output  valves  were  designed  to  balance  the  pressure  forces  on 
the  faces  of  the  disk.  Baffles  were  put  in  to  keep  the  jets  at  the  valve 
openings  from  hitting  the  opposite  sides  of  the  holes  and  creating  radial 
forces.  Leakage  flow  channels,  placed  around  each  of  the  driving  flow 
ports  in  the  port  plate  adjacent  to  the  disk  led  into  an  exhaust  passage 
in  the  plate.  By  placing  a  slight  restriction  in  the  exhaust  line  downstream 
from  the  vanes  and  keeping  the  leakage  flow  exhaust  line  unrestricted,  the 
axial  pressure  forces  in  the  center  of  the  disk  could  be  balanced  for  a  given 
pressure  and  approximately  for  neighboring  pressures.  This  was  not  done 
in  the  tests  of  the  oscillator,  however,  because  the  pressure  sensitivity 
was  already  quite  low. 

The  starting  point  in  the  design  was  the  choice  of  a  frequency,  a 
convenient  yet  small  amplitude,  and  a  maximum  allowable  wire  stress: 

co  ^  =  500  cps  -  7T x  10  rad/sec 
d  =  0.  010"{the  hole  size  is  roughly  the  amplitude) 

<r^  =  60,  000  psi  with  E  =  30  x  10^  psi 

The  layout  of  the  disk  came  next,  and  aa  a  result, 

c  2 

m  =  0.  258  x  10  lb.  sec  /in. 
n  =  8 

The  value  of  K  used  was  120  per  cent  of  the  desired  K  ;o  allow  for  an 
initial  20  per  cent  compression  of  the  wires  at  ambient  temperature. 
Equations  (233)  and  (235)  yielded  a  wire  size  of 

L  =  0.706  in. 
c  =  0.00831  in. 
diameter  =  0.0166  in. 

Area  =  2.  17  x  10**^  in.  ^ 

I  =  0.374  x  10"8  in.4 
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Fig.  97.  Temperature-Insensitive  Pneumatic  Oscillator 
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Investigation  showed  that  for  a  wire  compression  of  20  per  cent  of 
the  buckling  lead*  the  maximum  bending  moment  is  98  per  cent  of  that  of 
the  uncompressed  wire.  The  lesson  is  that  for  all  practical  purposes,  the 
maximum  stress  in  the  wire  can  be  calculated  as  the  sum  of  the  compressive 
stress  and  the  maximum  no-compression  bending  stress: 

cr  =  <r,  +  or 
max  b  c 

where  cr^  is  the  fluctuating  bending  stress  and  <rc  is  the  constant  compressive 
stress. 

In  the  actual  oscillator,  the  wire  sizes  were 

L  =  0.72  in. 

c  =  0.  0080  in. 

diameter  =  0.0160  in. 

.  -4  2 

Area  =  2.01  x  10  in. 

I  =  0.322  x  10"8  in.4 


We  recall  that  the  critical  loading  of  a  single  wire  is 


crit 


_  ir  ei 


(238) 


For  all  the  wires  together,  the  critical  load  is 


P  =  4  P' 
crit  crit 


(239) 


Consequently,  it  was  found  that 


P*  =  1.84  lb. 

crit 


0.2  Pl  =  0.  367  lb. 
crit 


0.  2  P« 


(T  = 
C 


=  1830  1b/in.2 


area 


or  -  (1.8  +  62)  x  103  lb /in.  “ 

max  —  ' 
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7.  34  lb. 


and 

p 

crit 

0.  2  P  =1.47  lb. 
c  nt 

An  estimation  of  the  linearity  of  the  3pring  can  be  obtained  using  the 
relation 

K  =  K 

o 


The  axial  stiffness  of  the  wires  is  very  great  in  comparison  with 
the  axial  stiffness  required  in  the  end  support.  Hence,  for  a  displacement, 
6,  of  the  disk,  the  axial  displacement  of  one  end  support,  the  other  being 
fixed,  is  2  6  /L.  The  column  loading  caused  by  this  displacement,  is  the 
product  of  the  axial  stiffness  of  the  support  and  displacement,  2  Ka&  /L. 
The  relative  change  in  the  lateral  stiffness  at  the  disk  is  therefore 


and  is  a  measure  of  the  deviation  from  linearity. 

In  this  test  model,  the  movable  end  support  was  suspended  on  four 
wires  each  0.47  inch  long  and  of  the  same  diameter  as  the  disk  spring  wires. 
For  this  case,  Eq.  (230)  yielded 


12  n  El 


52  lb/in 


The  mass  of  the  end  support  was  about  20  times  the  mass  of  the 
disk.  It  has  been  shown  that  for  this  situation  the  maximum  radius  for 
circular  motion  of  0.  65  times  the  disk  hole  diameter,  the  minimum 
amplitude  for  operation  is  about  0.  35  times  the  hole  diameter.  Then  for 
the  third  oscillator,  the  deviation  from  linearity  over  the  range  of  operation 


2K  {62  -  62  ) 

a  max  min 


L  P 


—  =  0.  05  per  cent 
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The  materials  used  were  chosen  for  availability  and  machinability 
and  would  not  be  suitable  for  high -temperature  applications.  The  disk, 
port  plates,  and  end  supports  were  brass,  the  springs  were  steel  music 
wire,  and  the  housing  was  structural  steel.  The  effective  diameter  of  the 
spring  wire  wa3  made  larger  where  it  passed  through  the  disk  by  forcing  a 
short  piece  of  hypodermic  tubing  over  a  slightly  flattened  portion  of  the  wire; 
the  tubing  was  then  pressed  into  the  disk.  Getting  the  wire  mountings  tight 
and  the  internal  damping  to  a  minimum  required  installing  the  wires  several 
times.  Even  then  the  wires  had  sightly  different  damping  in  different 
direction. 

6.  10.  Testing  of  the  Third  Oscillator 

Figure  100  is  a  representation  of  the  test  setup.  A  more  elaborate 
setup  would  have  been  desirable  but  would  have  taken  much  more  time  than 
was  available.  For  convenience,  frequency  change  was  measured  as  a 
function  not  of  compression  but  of  tension  in  the  spring  wires.  All  the  air 
connections  were  made  with  1/16  in.  I.  D.  plastic  tubing.  Properly,  the 
flows  as  well  as  the  pressures  should  have  been  measured  and  the  pressure- 
flow  characteristics  of  the  tubing  obtained  so  that  the  actual  pressures  at 
the  oscillator  could  be  determined.  Time  did  not  permit  this. 

Each  output  line  was  equipped  with  a  valve  in  series  and  a  bleed 
valve.  1  he  valves  in  aeries  allowed  the  output  lines  to  be  opened  to  the 
pressure  transducer  alternatively  or  in  combinations.  The  bleed  valves 
adjusted  the  pressure  and  flow  levels  of  the  output.  The  valves  were  small 
and  simple,  made  with  6-32  screws  in  a  plastic  block. 

There  were  three  provisions  for  adjusting  the  position  of  the  disk 
relative  to  that  of  the  facing  port  plates.  The  spacing  between  them  was 
adjusted  by  a  screw  through  the  base  spring  support  and  bearing  upon  the 
bottom  port  plate,  controlling  its  vertical  position.  The  clearance  between 
the  port  plates  was  fixed  by  spacers.  The  port  plates  could  be  rotated  in 
the  housing  to  adjust  angular  position.  The  lateral  position  of  the  disk  could 
be  adjusted  by  adjusting  the  position  of  the  top  support  for  the  springs.  This 
adjustment  was  extremely  difficult;  it  would  have  been  better  to  make  port 
plates  adjustable  laterally  so  that  their  position  in  all  directions  could  be 
easily  adjusted  by  screws  while  the  oscillator  was  running. 

-  205  - 


5CS" 


The  largest  pressure  signal  at  the  output  was  obtained  with  the 
output  line  bleed  valves  completely  closed.  The  leakage  flow  in  the 
oscillator  was  more  than  the  flow  required  for  maximum  output  pressure. 
The  only  output  measurements  taken  were  pressure  measurements  with 
the  bleed  valves  closed. 

The  output  measurements  were  very  disappointing.  Figure  101 
a  and  b,  show  five  results  for  two  different  input  conditions,  stated  there. 
These  were  taken  with  the  oscillator  in  the  same  condition  of  adjustment 
which  was  the  best  adjustment  obtained.  Lack  of  good  adjustment  allied 
with  inaccuracies  in  the  oscillator  parts  were  probably  responsible  for 
the  unequalness  of  the  different  outputs.  Notice  that  the  driving  force 
supply  pressure  for  both  the  output  measurements  is  quite  a  lot  larger 
than  the  output  supply  pressure.  With  these  two  pressures  equal,  the 
oscillator  would  not  run.  Considerable  damping  was  added  by  leakage 
flow  in  the  clearance  spaces.  This  damping  needs  to  be  taken  into 
account  in  the  design  of  an  oscillator  of  this  type. 

The  result  of  the  frequency  versus  wire  tension  test  is  given  in 
Fig.  i02.  These  data  were  taken  with  the  output  supply  shut  off.  The 
curving  of  the  line  could  well  be  the  result  of  error  in  measurement. 

Using  equation  (212)  and  the  value  of  P  =  34  lb.  obtained  for  this 

oscillator,  we  find  that  for  a  loading  of  10  oz. ,  K  should  increase  by  a 

i - 

factor  of  1.085  and  that  should  increase  byljl.085  =  1.042.  For  no 
tension  the  frequency  was  486  cps.  This  times  1.042  is  506  cps.  The 
measured  frequency  for  10  oz.  tension  was,  from  the  curve,  507  cps,  and 
from  the  data  point  itself,  108.5  cps. 

These  frequency  measurements  were  made  by  counting  the  number 
of  cycles  in  oscilloscope  photographs  of  leakage  flow  in  the  output  lines. 
The  sweep  time  was  2  milliseconds  for  19  cm  sweep.  Better  apparatus, 
such  as  a  precise  electronic  osciLlator,  could  beat  against  the  pneumatic 
oscillator,  indicating  the  short-time  stability  of  the  pneumatic  oscillator 
frequency  as  well  as  the  long-term  changes. 

With  the  output  air  supply  turned  off,  the  frequency  as  a  function 
of  driving  force  supply  pressure  was  measured.  As  mentioned  earlier 
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Fig.  101a.  Third  Oscillator  Output  (Tracings  from  Photographs). 
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-balancing  restriction  in  the  exhaust  line.  The  results 


are  gi  /en  in  the  following  t  .e , 


Driving  Force 
Supply  Pre33ure,  psig 
22 
26 
30 
34 
38 


Frequency,  cps 

484.  0 

485.  5 

486.  0 
485.  0 
484.  5 


The  frequency  figures  are  precise  to  not  more  than  +  0.  1  per  cent, 
or  0.  5  cps.  The  Largest  frequency  difference  is  2  cps  for  a  pressure 
difference  of  8  psi,  or  0.05  per  cent  per  psi. 

6.11.  Conclusions 

Because  of  the  output  of  the  third  oscillator  was  so  small,  one 
might  conclude  that  it  would  be  useful  only  for  driving  an  amplifier, 
probably  cf  the  fluid-jet  type.  A  diffeient  approach  would,  ther.  ore, 
be  to  assume  that  anything  vibrating  would  be  useful  if  amplifiers  were 
properly  applied  to  it,  and  then  to  design  something  to  vibrate  with  a 
constant  frequency.  An  oscillator  driven  as  the  second  and  third  os¬ 
cillators  were  driven  would  merely  use  the  exhaust  behind  the  vanes  to 
drive  one,  two,  or  more  amplifiers. 

A  possible  design  of  a  very  precise  oscillator,  and  also  one  with 
relatively  large  clearances,  is  given  in  Fig.  103.  The  amplitude  would 
not  be  large  enough  to  pull  the  inner  ends  of  the  slots  out  of  the  air  jets 
or  the  outer  ends  into  the  jet.  Pressure  forces  on  all  surfaces  per¬ 
pendicular  to  radii  would  be  nearly  the  same  all  around  the  disk.  The 
driving -force  versus  radius  curve  would  not  turn  downwards  in  the  range 
of  amplitudes  of  interest,  so  to  make  the  oscillator  operate  at  all  properly, 
a  jet  amplifier  control  could  be  necessary. 

Pitot  lubes  set  beneath  the  slots  at  the  proper  radius  and  connected 
together  would  sense  the  amplitude  and  cause  the  jet  supplying  the  oscillator 
to  flip  away  from  its  receiver  port  wnen  the  oscillation  amplitude  became 
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Fig.  103.  Jet  Amplifier-Controlled  Oscillator. 
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too  large.  A  pneumatic  lead  network  could  here  be  employed,  furnishing 
almost  ideal  amplitude  control  if  the  lag  inherent  in  the  amplifier  and 
lines  was  not  too  large.  As  shown,  output  pickups  would  mirror  the 
amplitude  “Sensing  pickups.  Axial  forces  would  be  balanced  by  virtue 
of  the  twin-disk,  design. 


Two  subjects  which  might  be  investigated  are  the  effect  of  external 
random  vibration  on  the  frequency  of  a  circular -motion  oscillator  and  the 
effects  of  slight  variations  of  spring  and  damping  constants  with  angle 
around  the  oscillator  axis. 


6.12.  Other  Oscillators 

Quite  different  oscillator  designs  were  briefly  considered  also. 

A  pneumatic  lag  causes  oscillation  in  the  device  shown  in  Fig.  104,  but 
the  driving  force  has  a  large  component  in  phase  with  the  flapper  dis¬ 
placement.  Since  the  lag  parameters  are  pressure  and  temperature 
dependent,  this  oscillator  has  an  additional  temperature  sensitivity  and 
would  need  a  good  temperature -insensitive  pressure  regulator.  With  a 
bellows -column  temperature  compensator,  it  might  be  used  if  the  re¬ 
quirement  for  precision  was  not  too  great.  It  has,  over  the  circular- 
motion  oscillator,  the  advantage  of  one -dimensional  motion.  Oscillation 
might  be  angular,  rather  than  translational,  so  that  a  properly  balanced 
oscillator  would  be  insensitive  to  linear  acceleration  in  any  direction. 

W*hen  the  flow  of  a  fluid  is  reduced,  pressure  forces  are  needed 
to  change  the  momentum;  these  pressure  forces  are  in  phase  with  the 
velocity  of  the  valve  causing  the  change  in  flow.  For  the  oscillator  shown 
in  Fig.  105,  the  motion  of  the  oscillating  member  toward  one  side  increases 
the  restriction  at  a  on  one  side  and  decreases  it  on  the  other  side.  The 
fluid  momentum  in  region  c  increases  on  one  side  and  decreases  on  the 
otheT*.  The  pressure  changes  necessary  to  produce  the  momentum  changes 
have  opposite  signs  on  the  two  sides  of  the  oscillating  member,  and  a  net 
lateral  force  results.  It  is  in  phase  with  the  rate  of  restriction  change  at 
a  and  the  velocity  of  the  oscillating  member.  By  adjusting  the  restriction 
at  b,  the  net  steady-state  pressure  moments  active  on  both  sides  of  the 
oscillating  member,  including  regions  c  and  d,  can  be  reduced  to  zero. 
Because  of  this  sensitive  adjustment,  this  oscillator  would  be  quite  sensi¬ 
tive  to  changes  in  the  downstream  impedance  and  would  have  to  be  used  in 
connection  with  some  form  of  amplifier.  Furthermore,  no  convenient  cut¬ 
off  technique  for  high  pressure  is  apparent. 
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Fig.  104,  Another  Oscillator. 
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7.  FLUID  DIODE  SENSOR 


7.1.  Introduction  and  Summary 

The  concept  of  fluid  diodes  used  tc  sense  the  width  of  the  wakes 
behind  compressor  blades,  along  with  initial  results  of  a  study  conducted 
by  W.  B.  Bails,  were  given  in  the  first  annual  report.  This  section  is  a 
continuation  which  assumes  the  reader  is  familiar  with  that  report. 

W.  B.  Bails  presented  his  thesis,  "Fluid  Diode  Sensor",  in 
September  1963  ,  and  performed  a  litule  analysis  thereafter.  Of  the 
three  ’large-scale"  diode  types  tested,  only  one,  the  "vortex  diode", 
gave  a  resistance  ratio  approaching  the  desired  value.  The  flow  was 
‘•urbulent.  For  corresponding  small-scale  models  the  flow  was  baacially 
lamine  r ;  all  the  resistance  ratios  were  reduced,  with  even  that  for  the 
vortex  diode  being  unsatisfactory.  Nevertheless,  dynamic  tests  were  run, 
with  the  realization  that  better  diodes  would  eventually  be  needed.  The 
vortex  diode  appeared  to  nave  a  low  cut-off  frequency,  rendering  it  useless 
for  the  high  frequency  response  required  for  blade -wake  sensing. 

Subsequently,  A.  A.  Khaliq  built  a  model  of  the  Tesla  diode  which 
was  brifely  tested  statically  by  F,  Brown.  The  static  operation  was  supsrior 
to  that  of  the  other  diodes,  although  a  turbulent  flow  is  still  essentially 
necessary.  Further  work  should  be  done  on  these  diodes,  including  dynamic 
testing;  they  seem  to  be  a  practical  basis  for  a  diode  sensor. 

7.2.  Dynamic  Analys:s 

The  lineax-ized  dynamic  analysis  was  presented  in  the  first  annual 
report.  The  nonlinear  square -law  analysis  was  carried  so  far  as  the  step 
responses  of  the  resistanc e -volume  combination,  tut  not  to  the  steady-state 
solution  tc  the  complete  square -wave  input.  This  latter  is  outlined  briefly 
below;  complete  details  have  beer,  given  by  Bails  in  his  thesis. 

Three  solution  domains  exist:  unsaturated,  partially  saturated,  and 
saturated.  In  the  unsaturated  case  the  individual  step  responses  never  reach 


+  Available  at  cost  from  the  Document  Room,  3-156,  Engineering 
Projects  Laboratory,  M.  I.  T.  ,  Cambridge,  Mass.  02139. 
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the  limiting  values.  In  the  partially  saturated  case  the  limiting  value  is 
reached  for  steps  in  one  direction  only;  it  Is  readied  for  steps  in  both 
directions  in  the  fully  saturated  case. 


The  unsaturated  case  is  shown  in  Fig.  112.  Matching  the  response 
segments  head-to-tail  gives,  in  the  steady  state,  the  deviation  between 
average  and  the  middle  value  of  pressure  of 
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(244) 


For  Dartial  saturation  in  the  forward  direction,  when 


<  ■)  ,  >  ->ft.  a  <  1 

a  -  2;  y  -  Zp;  —  -  -- 


Y  P 


the  solution  is 


P  -  P 


a  4  y  i  6 


2  ,2 

+  a  -  a  -  v  +  -  P 

3  \ 


t-  2a 


(245) 


In  the  negative  direction,  when 


a  *  1;  y  ~  2p;  -  *  - 

Y  P 


(246) 


the  solution  becomes 


P  -  P 


5 

Y 


a  +  -y 


7  /  .2 
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3  V  2P; 


Finally,  for  the  fully  saturated  case,  when 

a  ~  2 ;  y  ?  ?.p 


(247) 


(248) 


the  solution  is  simply 


P  -  P 


P,  4 

— -  (p  -  1)  +  a  -  y 


(249/ 


a+Y  L3  J 

1  /  1 

The  solutions  over-all  domains  for  a  particular  value  j  A-r/Tavj  =0.7 
are  plotted  in  Figs.  112  -  115.  In  Figs.  116  -  118  linear  and  nonlinear  analyses 
are  compared,  with  some  expeilmental  data  lor  a  small  value  of 
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7.3.  Static  Testing  -  Small-Scale  Models 

Three  models  of  the  nozzle  diode  were  made,  employing  two  basic 
de  ns.  The  two-dimensional  diode,  shown  in  Fig.  119b,  was  modeled 
directly  from  the  large-scale  water-table  models  using  three  cascaded 
nozzles.  The  othe  r  two  diodes,  shown  in  Fig.  119a,  were  three-dimensional, 
and  contained  six  nozzles  each.  The  vortex  diode  testes  is  shown  in 
Fig.  1  ZOa  and  the  diffuser  diode  in  Fig.  120b. 

The  flow  through  these  diodes  was  nominally  laminar,  as  shown  by 
Table  2  which  reports  the  maximum  Reynolds  numbers.  Since  the  flow  in 
the  large-scale  models  was  turbulent,  quite  different  performance  would  be 
expected,  and  in  fact  was  observed.  Experimental  results  are  shown  in 
Figs.  121  -  125.  As  mentioned  earlier,  only  the  vortex  diode  has  even 
marginally  satisfactory  static  characteristics.  Approximate  corres¬ 
ponding  values  of  the  linear  and  nonlinear  parameters  are  given  in  Table  3. 

7.4.  Dynamic  Testing  -  Small  Scale  Models 

A  schematic  of  the  test  setup  is  shown  in  Fig.  126.  The  slotted 
disk  pictured  in  the  first  annual  report  was  replaced  with  the  disks  pictured 
in  Fig.  127,  since  the  first  disk  wa3  deflecting  the  jet  significantly.  Typical 
resulting  dynamic  pressures,  measured  with  a  0,02-inch-diameter  Kistler 
piezoelectric  transducer,  are  shown  in  Fig.  128. 

The  final  experimental  results  °re  given  in  Figs.  116  -  118  and 
129  -  137.  Since  the  static  performances  w^ere  already  known,  the  principal 
interest  here  is  any  upper  frequency  limit  to  the  characteristics.  Such  a 
limit  is  decidedly  apparent  for  the  vortex  diode,  which  unfortunately  was 
the  only  diode  tested  by  Bails  with  even  marginally  satisfactory  static 
performance.  This  limit  apparently  is  inherent  to  the  transient  vortex¬ 
forming  process.  Detailed  discussions  of  the  experimental  procedure  and 
results  are  given  in  Bail's  thesis. 

The  failure  of  the  various  diodes  tested  has  led  to  the  later  work 
on  the  Tesla  diode. 

7.5.  Static  Testing  -  Tesla  Diode 

The  Tesla  diode  tested  is  shown  in  Fig.  138  and  139.  Tests  were 
performed  with  hydraulic  oil  to  achieve  low,  transitional,  and  high  Reynold’s 
numbers.  (The  characteristic  distance  used  tc  compute  the  Reynold's  number 
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Fig.  133.  Frequency  Response,  Vortex  Diode,  Forward. 
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Table  2.  Comparison  of  Maximum  Reynolds  Number 
for  Large  and  SmaLl  ScaLe  Models 


LINEAR  NONLINEAR 


136.  Frequency  Response,  Diffuser ,  Diode,  Forward. 


sure  difference  aP/P 


137.  Pressure  Difference,  Diffuser  Diode. 


Fig.  138.  Photograph  of  Tesla  Diode. 


is  the  narrowest  width  of  the  main  channel,  which  is  equal  to  the  depth 
of  the  model,  namely  one -quarter  inch). 

The  results  are  given  in  Fig.  139-  fn  the  laminar  flow  range 
(Re  <  100)  the  diode  is  almost  totally  ineffective.  Its  effectiveness  in¬ 
creases  with  Reynold's  number  as  expected;  the  limiting  value  of  the  ratio 
of  pressure  drops  at  high  flows  was  not  actually  found,  but  must  be  near 
5.0.  This  is  adequate  for  useful  fluid  diode  sensors,  although  a  somewhat 
higher  value  would  be  preferable.  More  sta£  -s  of  cusps  should  i.elp. 

A  fundamental  cause  of  the  fall-off  in  •  'rmance  of  fluid  diode 
sensors  at  high  frequency  can  be  traced  to  this  Reynold's  number  effect. 
As  frequency  ascends,  the  flow  rate  of  fluid  slushing  back  and  forth  is 
reduced,  and  the  effectiveness  of  the  diode  consequently  is  diminished. 

A  satisfactory  design  rule  would  be  to  expect  significant  attenuation  in 
the  performance  of  diode  sensors  when  the  peak  Reynold's  number  falls 
below  about  2000.  This  may  occur  at  an  insufficiently  high  frequency. 
Increasing  the  size  of  the  device  decreases  this  problem,  but  introduces 
wave -speed  limitations.  A  compromise  optimum  size  must  exist. 


Pressure,  drop,  psi 


Fig.  139'  Results  of  Static  Measurements  of  Tesla  Diode. 
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